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Abstract 



This thesis consists of two parts. In part I we consider a discrepancy in the deriva- 
tion of the electromagnetic self force for a point charge. The self force is given 
by the Abraham-von Laue vector, which consists of the radiation reaction term 
proportional to the 4-acceleration, and the Schott term proportional to the 4-jerk. 
In the point charge framework the self force can be defined as an integral of the 
Lienard-Wiechert stress 3-forms over a suitably defined worldtube. In order to 
define such a worldtube it is necessary to identify a map which associates a unique 
point along the worldline of the source with every field point off the worldline. One 
choice of map is the Dirac time, which gives rise to a spacelike displacement vector 
field and a Dirac tube with spacelike caps. Another choice is the retarded time, 
which gives rise to a null displacement vector field and a Bhabha tube with null 
caps. In previous calculations which use the Dirac time the integration yields the 
complete self force, however in previous calculations which use the retarded time 
the integration produces only the radiation reaction term and the Schott term is 
absent. We show in this thesis that the Schott term may be obtained using a null 
displacement vector providing certain conditions are realized. 

Part II comprises an investigation into a problem in accelerator physics. In a 
high energy accelerator the cross-section of the beampipe is not continuous and 
there exist geometric discontinuities such as collimators and cavities. When a 
relativistic bunch of particles passes such a discontinuity the field generated by a 
leading charge can interact with the wall and consequently affect the motion of 



trailing charges. The fields acting on the trailing charges are known as (geometric) 
wakefields. We model a bunch of particles as a one dimensional continuum of point 
charges and by calculating the accumulated Lienard-Wiechert fields we address the 
possibility of reducing wakefields at a collimator interface by altering the path of 
the beam prior to collimation. This approach is facilitated by the highly relativistic 
regime in which lepton accelerators operate, where the Coulomb field given from 
the Lienard-Wiechert potential is highly collimated in the direction of motion. It 
will be seen that the potential reduction depends upon the ratio of the bunch 
length to the width of the collimator aperture as well as the relativistic factor 
and path of the beam. Given that the aperture of the collimator is generally on 
the order of millimetres we will see that for very short bunches, on the order of 
hundredths of a picosecond, a significant reduction is achieved. 
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Guide to Notation 



All fields will be regarded as sections of tensor bundles over appropriate domains of 
Minkowski space A4. Sections of the tangent bundle over A4 will be denoted FTA^ 
while sections of the bundle of exterior p-forms will be denoted FA^A^. Given a 
single worldline C in free space sections over the whole of spacetime excluding 
the worldline will be written rT(A^\C) and rAP(A^\C). We use the SI unit 
convention. Appendix |A] provides a brief summary of the dimensions of various 
mathematical objects. 



Chapter 1 



Introduction 



In this chapter we introduce the defining characteristics of Minkowski space, namely 
the metric and the affine structure, and the fundamental equations of Maxwell- 
Lorentz electrodynamics. We use the term Maxwell- Lorentz electrodynamics to de- 
note the microscopic vacuum Maxwell equations, first derived by Lorentz from the 
macroscopic Maxwell equations (see [2111]) and often called the Maxwell-Lorentz 
equations, together with the Lorentz force equation. We introduce the general 
form of the electromagnetic stress 3-forms and show that they give rise to a set 
of conservation laws. A brief introduction to the necessary mathematics can be 
found in Appendix [Bj 

1.1 Minkowski space 

Definition 1.1.1. Minkowski space is the pseudo-Euclidean space defined by 
the pair {Ai,g), where Ai is the four dimensional real vector space M^ and g 
is the Minkowski metric. With respect to a global Lorentzian coordinate basis 
{y^ , y^ , y"^ , y^) on Ai the Minkowski metric is defined by 

g = -dy° (g) dy° + dy^ (g) dy'^ + dy"^ ® dy^ + dy^ (g) dy^. (1.1) 

Lemma 1.1.2. Given a new set of coordinates {z^, z'^, z"^, z^) on Ai we write g in 



Minkowski space 



terms of the new basis using the transformations (B.19), 



Qya Qyb 






g = gabdy" ® dy' = gab^^dz' A dz" = g^^'dz' A dz"" (1.2) 



where 



^c. -^a.^^. (1.3) 



Lemma 1.1.3. The coordinate basis {y^ , y^ , y"^ , y^) on Ai naturally gives rise to 
the basis {dy^ , dy^ , dy"^ , dy^) of 1— forms on T*A^. This forms a g-orthonormal 
basis and from definition \B.2.5 it follows 



-kl = dy" A dy^ A dy^ A dy"^ (1.4) 

is the volume form on M.. In terms of a different coordinate basis {z'^ , z^ , z'^ , z^) 
it is given by 

7^1 = ^/ \det {g')\dz^ A dz^ A dz^ A dz^ (1.5) 



(^) 



where g^ is the matrix of metric components g^l defined by (1.3). 



Definition 1.1.4. Let ^ be a vector field over set of points M. If tlie map 

MxM — >V, {x,y) — >x-y, (1.6) 

exists and satisfies 

1. For all X e M, for all v eV 

there exists y G M such that y — x = v, 

2. For all x,y,zEM, {x — y) + {z — x) = z — y, (1.7) 

then M is an afiine space. 



Maxwell-Lorentz Equations 



For any integer n the space M" is affine. It follows that Minkowski space is an 
affine space. In some calculations it will be necessary to endow Minkowski space 
with an origin, thus transforming it into a vector space, however the results of 
such calculations will not depend on the vector space structure but only the affine 
structure. 

1.2 Maxwell-Lorentz Equations 

The equations which describe the interaction between matter and the electromag- 
netic field were first formulated by Maxwell in 1865 [5]. Maxwell's equations form 
a continuum theory of electrodynamics due to their origins in macroscopic experi- 
ment. In this thesis we are interested in the interaction of point charges and their 
fields, therefore we need equations which are valid on the microscopic scale. 

Definition 1.2.1. The Maxwell-Lorentz equations, or the microscopic vacuum 
Maxwell equations are given by 

dJ^ = (1.8) 

eod*J^ = J (1.9) 

where J-" G FA^A^ is the electromagnetic 2-form, J' G FA^A^ is the current 3-form 
and eo is the permittivity of free space. If we introduce the 1-form potential 

^gFA^A^ such that J^ = dA, (1.10) 



then (1.8) and (1.9) reduce to the single equation 

eodi<dA = J, (1.11) 



where (1.8) is satisfied automatically because the double action of the exterior 



derivative is zero. 



Maxwell-Lorentz Equations 



Definition 1.2.2. In terms of distributional forms (see C.1.2) 



A"" G TdA^M, 7"" = dA"", 



(1.12) 



Maxwell's second equation is given by 



eod^dA''[ip] = J''[ip]= I ipAj. 

'm 



where cp G TqA^M is any test 1-form (see C.1.1). 



:i.i3) 



3 + 1 decomposition 

Definition 1.2.3. Given any velocity vector field U G FTA^ satisfying 
g{U, U) = —1, the electromagnetic 2-form J-" may be written 



T = g AU + cB, 



:i.i4) 



where £ G FA^A^ and B G FA^A^ are the electric 1-form and magnetic 2-form 
associated with U and J-", and satisfy 



iu£ = "it/B = 0. 



:i.i5) 



Here is the metric dual operator defined by (B.29) and c is the speed of 

light in a vacuum. 

Lemma 1.2.4. According to observers whose worldlines coincide with integral 
curves of U, the electric field S G FTA^ is given by 



E =iuJ^, 



;i.i6) 



Proof of 1.2.4- Follows trivially from (1.14) 
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Maxwell-Lorentz Equations 



Definition 1.2.5. We may use the vector field U to write tlie Minkowski metric 
g in terms of a metric g on tlie instantaneous 3-spaces 

g = -U(S)U + g. (1.17) 

Let # be the Hodge map associated with the instantaneous 3-space such that for 
a e TAPM 

# : TA^M -^ TA^-PM, a^i^a = {-lY+Hu * a (1.18) 

The Minkowski Hodge dual is then given by 

^a = (-l)Pf7A#a. (1.19) 

Lemma 1.2.6. The Hodge dual oj T is given by 

^jr=^^_c#BA?7 (1.20) 

Proof of 1.2.6. 



i.J^ = ^{S AU) + c*B, (1.21) 



It follows from (1.18) that if « G A^A^ then #« = iu i^ a = i^{a A U). Thus 



^(^Af/) = #^. Similarly it follows from (1.19) that if /3 G A^M then */3 = f/A#/3 



thus *B = U A #B. D 

Lemma 1.2.7. Let B = — #B where B G FTA^ is the magnetic field, then ac- 
cording to observers whose worldlines coincide with integral curves of U it is given 
by 



B = -U^Zf, . (1.22) 

c 



Maxwell-Lorentz Equations 



Proof of 1.2.7. Consider (1.20). Since iuifct = iuifP = it follows that iu -^^ J^ 

-c#B. 

Lemma 1.2.8. In terms of S and B the 2-forms T and -kj-" are given by 



D 



i,jr =^s^cBA U = i.{SAU) + cBA U. 



:i.23) 
:i.24) 



Proof of 1.2.8. Since i3 = -#B and ##B = B it follows that #i3 



B. 



Substituting this into ( |1.14[ ) yields (1.23). Similarly substituting the first relation 

D 



into (1.20) yields (1.24). 



The Lorentz Force 

Definition 1.2.9. Let C : / C M — t- A^ be the proper time parameterized inex- 
tendible worldline of a point particle with observed rest mass m and charge q. For 
re/ 

C = C,{d/dT), C = VcC, and c = VcVcC (1.25) 

are the velocity, acceleration and jerk of the particle respectively. Here the push- 



forward map =„ is defined by (B.79)-(B.83) and V is the Levi-Civita connection (see 



B.2). In this introductory chapter and in Part II we assign the dimension of time 



to proper time r such that 



g{C,C) = -c' 



1.26) 



However the reader should note that Part I we will find it convenient to assign the 
dimension of length to proper time so that 



g{C,C) = -l, 



:i.27) 



Conservation Laws 



For further details about dimensions see appendix \K\ 



Lemma 1.2.10. 



g{C,C) = 0, 
and giC,C) = -giC,C). 



1.28) 
1.29) 



Proof of 1.2.10. Equation (|1.28|) follows by differentiating (|1.26) with respect to 

D 



r. Similarly, equation (1.29) follows by differentiating (1.28). 



Definition 1.2.11. The force on a point particle with worldline C{t) due to an 
external field Fext £ FA^A^ is given by the Lorentz force /l, where 



/l e FTA^, /l 



Q. 



"^C ext- 



:i.3o) 



In 1916 Lorentz writes [6] 

Like our former equations [Maxwell's equations] , it is got by generaliz- 
ing the results of electromagnetic experiments 

1.3 Conservation Laws 

Definition 1.3.1. A vector field \^ is a Killing field if it satisfies 



Cvg = 0. 



(1.31) 



In terms of coordinate basis {y^} the metric may be written g = gab{y^)dy°' ® dy'', 

a 



thus for vector field V = J^ the left hand side of (|1.31|) yields 



^9k9 = -Q^dy" (g) dy\ 



1.32) 



where Ok = -J^- In Minkowski space g^Q = — 1 and gab = S^ for a = 1, 2, 3. Thus 



dy^ 



for the four translational vectors g^, tJ^, ^, ^ (|l-31|) is trivially satisfied. In 



Conservation Laws 



fact there are 10 killing vector fields on Minkowski space. 

Let V^ be a Killing vector, then another property of Killing fields we shall use 
is 

Cyi^ = i,Cy^ (1.33) 

Definition 1.3.2. The electromagnetic stress 3-fornis iSk G FA^A^ are given by 



S^=^^^{l9^J^^^7-^9^^7^7) (1.34) 



where c?k = qt^ are the four translational Killing vectors. These 3-forms can be 
obtained from the Lagrangian density for the electromagnetic field using Noether's 
theorem, see [^ for a detailed exposition. The stress forms are related to the 
symmetric stress-energy-momentum tensor T G V^^ ' 'A^ by 

r°" = !^*5K, 5k = *((r(d/. -))-) (1.35) 



Lemma 1.3.3. The stress forms satisfy 



dSK = —td^J^Aj, (1.36) 



c 



and thus for any source free region N d M. 

dSK = 0. (1.37) 



Conservation Laws 



Proof of 1.3.3. 



dSK =^d{iQ^F /\ -kT — iQ^-kT l\T^ 



"2c 



[diQ^T A kT — iQ^T f\d^ T — diQ^ ^ T f\ T ^ iq^^^ T f\ rfj-") 



eo 



[diQ^T A kT — diQ^ -k J^ /\ T^ — —idK-F Ad-k J^. 



:i.38) 



From (1.8) and (B. 73) it follows that 



9k -^ uZ^y". 



;i.39) 



Using (1.39), (B.2.11) and (1.33) respectively yields 



diQ^J^ A kj^ = J^ A kdiQ^J^ = J-" A kCg^J^ = J-" A Cq^^ k J^ = J^ a dlQ^ k J^ + J^ A iQ^d k T 

(1.40) 



Substituting (1.40) into (1.38) yields 



dSvi = 7^ (-^ A d%Q^ k T + T AiQ^dk T — d%Q^ kT AT) — —idi^T A dk J^ 



eo 
"2c 



J-" A iQ^d k J^ — iQ^J-' Adk J^) 



:i.4i) 



Since J-" is a 2-forni and d-k J^ is a 3-forni it follows that 



ig^i^F A dk J^) = iQj^J^ Adk J^ + J^ A iQ^dk J^ = 0. 



:i.42) 



Thus substituting J^ AiQ^dkJ^ = —iQ^J^AdkJ^ and (1.9) into (1.41) yields result. 

D 



Lemma 1.3.4. For any source free region N G Ai 



Sk = dSx = 0. 

dAf JM 



:i.43) 



Proof of 1.3.4- Follows trivially from (1.37) and Stokes' theorem (B.103). D 
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The source J for a point charge 



Lemma 1.3.5. If U is a timelike Killing vector then applying the 3+1 decompo- 
sition yields 



Su = eoSABAU+ ^{S A #f + c^B A #S), (1.44) 



where cqS A B is the Poynting 2-form, and |^(<^ A ^£ + c^B A #i3) the energy 
density 3-form. 



Proof of 1.3.5. 



Using definition 1.3.2 



Su =y{iduJ^ A-kJ^ -id^j-kj^ AF) 



Substituting (1.24) and (1.23) and using the relations (1.16) and (1.22) yields 



eo 



Su =^{£ A {i^£ + cB AU) ~ cB A {£ AU - c^^B)) , 



eo 
"2c 



(^ A #^ + c^B A #i3) +e^E ABAU. 



D 



1.4 The source J^ for a point charge 

We now consider the particular form of the current 3-forni J G A^A^ for a point 
charge. We use notation J = jTpoint charge in order to emphasize that J is a particular 
choice for J . The source is located only on the worldline of the particle therefore 
we expect the source distribution J^ G F/jA^A^ to have the form of a Dirac delta 
distribution. 

Definition 1.4.1. Given the four 0-form distributions 2°" G F/jA^A^, where for 
X ^ M 



3 \^) = q 



C"^(r)5(^)(a;-C(r))t/r, 



:i.45) 



11 



The source J for a point charge 



we define the distributional current vector field by 

J=/(x)|„ (1.46) 

The distributions j°'{x) are non-zero only when x = C{t). The 3-form J G A^A^ 
is given by 

J = ^j. (1.47) 

Lemma 1.4.2. The current 3-form distribution J^ G F/jA'^A^ is given by 

J^M=g /CV (1.48) 

for any test 1-form ip G TqA^M. 
Proof of 1.4.2. 



I 



From (1.47) the distribution J^ is given by 



Jm 

ij -kl A ip, 



M 

fia^-kl, 
M ^" 

J Va * 1. 
M 
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Worldline geometry 



Substitution of (1.45) yields 



iDf 



J M Jt 

/Jfia 
Va{C{T))—dT 

qfMC{T))dC'' 
q Ua{C{r))C\dy^) 



q I C*^ 



where C*(y") =y''oC = 0". 



D 



1.5 Worldline geometry 



Given the proper time parameterized inextendible worldhne 



C -.1 cR^M, t^C{t), 



:i.49) 



we required a way to locally map each point x G {Ai\C) to a unique point C{t'{x)) 
along the worldline. Consider the region N = N\C where A^ C A^ is a local 
neighborhood of the worldline. The afiine structure of A^ permits the construction 
of a unique displacement vector Z\^ defined as the difference between the two points 



(see figure 1.1), 



Z\, = x-C{r'{x)). 



:i.5o) 



Note that the definition of Z only requires the affine structure of A^. It does not 
require Ai to be converted into a vector space by assigning an origin. 
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Worldline geometry 




Figure 1.1: Displacement vector Z\ 



We may construct a local vector field Z G FTA^ such that 






where 



x^-Cir'ix)). 



:i.51) 



for all X E N. Here y"'{x) = x"". 

Since Z is defined for every x E N the only requirement needed to define Z 
completely is to fix r'(x). We are free to choose t'{x) in any way we like however 
particular choices are beneficial for certain problems. In one choice the vector Z 



lies in the plane perpendicular to C{t'{x)) (see figure 1.2). In this case we use 
the notation t' = td, where r/5 is the Dirac time. The Dirac time associates each 
point X E N with the time td{x) given by the solution to 



g{x-C{Tn{x)),C{TD{x)))=0. 



:i.52) 
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Worldline geometry 




D\x 



Let the norm ||.|| be defined by \\Z\\ = \/ g{Z^ Z^. 
Then ||y]||| =0, and ||F|| = \\Y^\\ = g{Y,Y). 



Figure 1.2: Displacement vector Y\x- 



In this case Z\x = Z± = x — C{t£){x)). We use the special notation 



Y = x-C{td{x)). 



:i.53) 



The map r^ : A^ — )■ C is not unique for every a; G A^, for example in figure 1.3 



we see that a single point can be mapped to multiple points along the worldline. 
However for a sufficiently small neighborhood N C {Ai\C) uniqueness can be 



ensured. In appendix D we explore this geometry further. 



In another choice the vector Z lies on the null cone (see figure 1.4). In this 
case we use the notation r' = r,., where r.^ is the retarded time. The retarded time 
associates a point x G {Ai\C) with the time Tr{x) given by the solution to 



g{x- C{Tr{x)), X - C{Tr{x))) =0, x'^ > C°(r,(x)) 



:i.54) 
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Worldline geometry 



plane perpendicular,,,— 

to C\r'p 

C 



C^v 



,6tV' 



etv' 



.^;^c^ 



A&t 



TO 




Figure 1.3: Globally the map td is non-unique. 



In this case Z\x = x — C{Tr{x)). We use the special notation 



X = X-C{Tr.{x)). 



:i.55) 



There is another possible choice in which Z is a vector in the advanced null 



cone at x (see figure 1.4). In this case we use the notation t' = Ta, where r^ is the 
advanced time. The advanced time associates a point x G {M.\C) with the time 
Ta{x) given by the solution to 



g{x- C(r„(x)), X - C(r,(x))) =0, x° < CO(r,(x)) (1.56) 



In this case Z\x = x — C{Ta{x)). We use the special notation 



W = X-CiTa{x)). 



1.57) 



The maps r,.(x) and Ta{x) are not necessarily defined for all X G A^. Figure 1.6 
shows the path of a curve undergoing constant acceleration. The backwards light 

16 



Worldline geometry 







Let the norm ||.|| be defined by \\Z\\ = ^/g{Z, Z)'^. 
Then ||X|| = and ||X|||| = \\X^\\ = g(X,CLMf. 

II II II ||ll II -Lll J \ 7 \ I ry-f-- ) / 



Figure 1.4: Displacement vector X|a.. 




Figure 1.5: Displacement vector W\ 
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C(t) 




# 



Figure 1.6: Globally the retarded and advanced times are not necessarily defined 
for all X G iM\C). 

cone from an arbitrary point in quadrant A or B intersects the worldline once in 
quadrant B, hence the retarded map Tr{x) is well defined in A and B. However 
the backwards light cone from any point in quadrants C or D will never intersect 
the worldline, therefore the map Tr{x) is not defined for x G C or x G D. We can 
ensure the existence and uniqueness of the maps r^ and r^ by working exclusively 
in a sufficiently small (and appropriately chosen) neighbourhood N C {Ai\C). 



Null geometry 

In this section we explore further the consequences of choosing t' = Tr. This map 
is particularly suited to electromagnetic phenomena which propagate on the light 
cone. We call the resulting geometry null geometry. We begin by consolidating 



equations (1.49), (1.51), (1.54) and (1.55). 
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Definition 1.5.1. Given the one-parameter curve C(r) which traces the path of 
a point charge in spacetime, then for every field point x G {Ai\C) there is at most 
one point Tr{x) at which the worldline crosses the retarded light-cone with apex 
at X . 

C:R^M, t^C{t) (1.58) 

Tr-. M^R, X ^^ Tr{x) (1.59) 

Definition 1.5.2. The null vector X G rT(A^\C) is given by the difference 
between the field point x and the worldline point C{Tr{x)) 

X\, = X~C{Tr{x)), (1.60) 



where 



g{X,X) = g{x - C{Trix)),x - CiTrix))) =0. (1.61) 

Definition 1.5.3. The vector fields V,A,Ae rT(A^\C) are defined as 

VU = C^{Tr{x)) — , AU = C^{Tr{x))— and A\, = &{Tr{x)) — , (1.62) 
oy^ oy^ oy^ 

hence from lemma 11.2.101 it follows 

g{y,V) = -c\ g{A,V)=Q, and g{A,V) = -g{A, A). (1.63) 

Definition 1.5.4. We define the normalized null vector field by 



V 

n=—, where R=-g{X,V) (1.64) 

R 
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The normalized vector satisfies 



g{n,n) = and g{n,V) = —1. 



(1.65) 



Lemma 1.5.5. The exterior derivative of the retarded proper time Tr is given by 



dTr 



X 



9{x,vy 



;i.66) 



Proof of 1.5. 5. 



Definition 1.5.2 requires only the affine structure of A^. For the following proof we 
demand the stronger requirement that the points x E M. and C{T.r{x)) G C{t) C 
Ai are attributed with a vector structure on A^, such that 



xgTTTW, x|^=a;" 



dy" 



and C G TTM, C^ = C"^(r,(x)) 



d 

(1.67) 



however the result (1.66) requires only the afiine structure. 



We begin with the light cone condition 



= giX,X), 
= (7(x-C,x-C), 
= ^(x-C,x)-(7(x-C,C), 
= ^(x,x)-2^(C,x) + (7(C,C) 



[1.68) 



Therefore 



= d 



^(x,x)-2(7(C,x) + (7(C,C) 



dg{^, x) - 2dg{C, x) + dg{C, C) 



:i.69) 
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Now the first term in (1.69) yields 



(i5f(x,x) = d{gabx'^x^), 

= gab{dx'')x^ + gabx'^idx''), 

Note that x'* = y°'{x) thus dx°' = dy"- and therefore 



(i5f(x, x) = Xady'' + Xady", 
= '2xady'', 
= 2X. 



:i.7o) 



Similarly the second term yields 



dg{C,^) = d{gabx''C\Tr.)), 

= gMx'')C\rr)+gabx''d{C\T,)\ 



where rf(C"(r,)) = ^(C"(r,))rfr, = VdTr, therefore 

dTr 



:i.7i) 



dc/(C,x) = C + XaV^''rfr„ 
= C + 5f(x, V)dTr 



:i.72) 



The third term gives 



dg{C,C) = d{gabC\r,)C\Tr)), 

= {dC\Tr))gabC\Tr) + {dC\Tr))gabC\Tr 
= 2CaiTr)V''dTr, 

= 2g{C,V)dTr. 



1.73) 
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Substituting (1.70), (1.72) and (1.73) into (1.69) yields 



therefore 



= 2x - 2(^0 + ^(x, V)dTr) + 2g{C, V)dTr 



St-C = g{^-C,V)dTr, 



and upon rearrangement yields 



dTr 



x-C 



X 



g{X,V) g{X,Vy 



:i.74) 



Lemma 1.5.6. 



D 



dV = dTr A A 



:i.75) 



Proof of 1.5. 6. 



Corollary 1.5.7. 



dV =^dTrA dy" 

dTr 

= AadTr A dy" 

= dTr A A 



D 



dV 



X AA 
9{X,V) 



(1.76) 



Proof of 1.5.7. Follows directly from (1.74) and (1.75) 



D 
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Lemma 1.5.8. 



^i*V) ^ fll^ * 1 (l.TT, 



Proof of 1.5.8. 



dya 



A°-dTr Ai_a_ •! 



dTr A i^A 

X Ai<A 
9{X,V) 
9{X,A) 
9{X,V) 



^1 



Lemma 1.5.9. 
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D 



,M.?An^™£Ai)-3*V- (1.78) 

9[-^, V ) 



Worldline geometry 



Proof of 1.5.9. 

di<{X AV) = di. {Xady" A V^dy^) 



diYbXa^idy'^Ady')) 

d{VbXa) A ^(rfy" A dy^) + V^X^rf^ {dy'' A dy^) 

dVb A X„ ^ {dy'' A dy^) + H^^^a A *(rfi/" A d?/'') 

rf\4 A ^(X A rfy^) + dXa A *(rfi/" A y) 

"^^"^ rfr, A ^(X A dy') - ^((^"'X,) A 29 „ ^ K 



Aforfr^ A 7^(X A rfy^) - rfX"" A ia^„ • K 

rfr, A ^(X A I) - dy'' A ig^^ ^V + dC'A Iq^^ ^ V 

dr., A 7^(X A 1) - rfy" A ia^^ -kV + dTr A -k{V A V) 
X A ^(X A I) 



Need also to show that 



dy'' A iQ. -k V 



dy''Ma^a^V = ^^V (1.79) 

Let V = Vady"", then 

^V = VaQ'^Ho , * 1 = V'ta , ^ 1 

a° y° 

Substituting -kl = dy^ A dy^ A rfy^ A dy^ and contracting yields 

^\/ = Vody^ A rf|/2 A dy^ - Vidy^ A rfy^ A dy^ + Varfy" A dy^ A rfy^ - Kgrfy^ A dy^ A rfy^ 

(1.80) 
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Thus 



dy" Aidya^V = - Vidy° A dy'^ A dy^ + V2dy^ A dy^ A dy^ - Vsdy^ A dy^ A dy^ 

- Vody^ A dy"^ A dy^ - V2dy^ A dy° A dy^ + Vsdy^ A (iy° A rfy^ 
+ Vody^ A rfy^ A dy^ + Vidy"^ A dy° A dy^ - V^dy^ A d?/° A dy^ 

- V^oC?2/^ A dy^ A dy^ - Vidy^ A rfy° A dy^ + V2dy'^ A (iy° A dy^ 

(1.811 



Collecting terms in (1.81) and comparing with (1.80) yields (1.79) 



Lemma 1.5.10. 



^.(XaI)^^A^E4^ + ^M44^-3.A 



9{X,V) 



9{X,V) 



Proof of 1.5.10. 



d * (X A I) = d{XaAb * (dy" A c^y'')) 

= d{XaAb) A 7^(rfy" A c?y^) + XaAbd i. {dy'' A dy^) 
= {dXa) AAb-k {dy" A dy^) + X^rfA^ A -kidy" A t/y'') 
= t/X^ A -k{dy'' A 1) + rfAfc A ^^(X A dy^) 
= -d^g^^Xa) Mo ,i^A + AhdTr A .^(X A dy^) 

X A ^(X A A) 



-dX" A za,a * A + 



^(^,V^) 



-dy'* A ioya kA + dC A io^, *A + 
-3i<A + VdTr Aia ,i^A + 



X A *(X A A) 



9{X,V) 

XA^(XAi) 



-3i.A 



9iX,V)^ 

XA^AAV) XA*(XAi) 



9{X,V) 



9{X,V) 



D 



1.82) 



D 
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Lemma 1.5.11. 



dg{A, A) = 2g{A, A)dTr 



:i.83) 



^^^-'-^-^H'^^M^'^- 



:i.84) 



Proof of 1.5.11. 



Proof of (1.83) 



^^^^'^)-^~+(g^'^ 



:i.85) 



dgiA,A)=digatC''iTr)C\Tr)) 
= 2gabdC''C' 

= 2gabC"'CbdTr 

= 2g{A, A)dTr 



Proof of (1.84) 



dg{X,V) = dgix-CiTr),V) 



d 



g{x,V)-g{C{rr),V) 
dg{x,V)-dg{C{Tr),V) 
d{gabX^V')-d{g,bC''{Tr)V') 

gabV'dy'' + gabx'^dV' - gabV'dC^Tr) - gabC%rr)dV' 
V + g{x, A)dTr - g{V, V)dr,. - g{C{Tr), A)dTr 



V + 

V + 



g{x-C{Tr),A)-g{V,V) 



dTr 



giX,A)-giV,V) 



dTr 
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Substituting (1.26) yields result 



Proof of (1.85) 



dg{X, A) = dg{x - C{Tr), A) 



d 



g{x,A)-g{C{Tr),A) 
dg{x,A)-dg{C{Tr),A) 
digabx'^A') - diga,C''iTr)A') 

gabA'dy'^ + gabX^dA' - g^tA'dC^iTr.) - gatC^{Tr)dA' 
A + g{x, A)dTr - g{A, V)dTr - g{C{Tr), A)dTr 



A + 
A + 



gix~CiTr),A)-giA,V) 



dTr 



giX,A)-giA,V) 



dTr 



Substituting (1.28) yields result. 



D 



1.6 Newman- Unti coordinates {t, R,6, 4>) 

We introduce a system of coordinates adapted to the null worldline geometry. The 
coordinates were first introduced in a general form for arbitrary manifolds by Tem- 
ple in 1938 [8j, where they are referred to as optical coordinates. In 1963 Newman 
and Unti [S] claim to introduce a new coordinate system "intrinsically attached to 
an arbitrary timelike worldline", however the coordinate system they investigate 
is none other than the specialization of Temple's coordinates to Minkowski space. 
Since in this thesis we work explicitly with Minkowski space we have chosen to re- 
fer to the coordinates as Newman-Unti (N-U) coordinates in the spirit of Galt'sov 
and Spirin [TU], however the general class of coordinates should be attributed to 



Temple. Similar coordinates were used by Trautman and Robinson [IT] in their 
work on gravitational waves, and in the 1980's Ellis [I2] and others use similar 
coordinates in problems in relativistic cosmology where they are called Observa- 
tional coordinates. Other variations on the name include retarded coordinates, null 
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Newnian-Unti coordinates (t, R, 6, 



geodesic coordinates and lightcone coordinates. 



We recall from (1.60) that 



X = x-C{i 



d 



d 



d 






a 



:i.86) 



Definition 1.6.1. Given the global Lorentzian frame {y^ , y^ , y^ , y^) on Ai, the 
Newman-Unti coordinates (r, R,6,(f)) are defined by the coordinate transformation, 



y = Ci(r)--sin(^)cos(0), 
a 

y^ = C\T)--sm(e)sm((j)), 
a 



and i/=^ = C=^(r)--cos(^), 



(1.87) 



where a G FA^A^ is defined by 



a(r, 6', I 



g{X,dyo) 
- C\t) + C\t) sin(e) cos((/)) + C\t) sin(e) sin(0) + ^^(r) cos(^). 

(1.88) 



From (1.87) and (1.88) it follows 



R=-g{X,C{T)) and r = r^(x(r, i?, 0, 0)). 



1.89) 



The spherical coordinates 6 and are given naturally from the global Lorentzian 
frame {y'^ , y^ , y^ , y^) . 

Lemma 1.6.2. In Newman-Unti coordinates the vector fields X and V are given 
by 



X = R 



d_ 



:i.9o) 
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Newnian-Unti coordinates (t, R, 6, 



and 



OT a 



(1.91) 



Proof of 1.6.2. 



Proof of (1.90) Differentiating the coordinate transformation (1.87) with respect 



to R yields 



_d_ 
OR 



dy^ d dy^ d dy"^ d dy^ d 



OR dy^ 



dRdy^ dRdy^ dRdy^ 



a 



Therefore 



X = x- C{t) 



i? 



ai? 



Proof of (1.91) 



d dy^ d dy^ d dy"^ d dy^ d 
dr dr dy^ dr dy^ dr dy"^ dry dy^ 
d 



C\r) - R^/-)) A + [c\r) - RsiniO) cos(0)^rl^^ ^ 



+ CHr) 



QyO ' y^ "- ' "" '-" ^ "-y-rjg^ ^ J dy^ 
i.s.n(.)sm(0)^(l))|,+ (c3(.)-i.cos(.)^(l))A 






--V-X 



a 



a 



OT a 



D 
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Newnian-Unti coordinates (t, R, 6, 



Lemma 1.6.3. In Newman-Unti coordinate the Minkowski metric g E 0' ' ^ M 
is given by 



g =(-c^ + 2R-)dT ®dT - (dr ®dR + dR® dr) 
a 

T-)2 d2 

+ -^dO (^de + ^^ sm{e)'^d(f) (g) d(f), 



a^ 



a^ 



1.92) 



and inverse metric g ^ E <^ ' M is given by 



[V,V] 



-1 



-c'^a + 2Ra d 



d a^ d 
+ 



d 



+ 



a 



a 

2 



dR dR R^ de de 

d d , d d d 

+ 



d 



i?2 sin(^)2 dcj) dcj) ^dr dR dR dr'' 



;i.93) 



Let z^ = T, z^ = R, z^ = ^, z^ = 0, then the matrices G = Gab = g{dz'^, d^t) 
and G^^ = G'^^ = g^^i^dz"", dz'^) are given by 

\ 



and 



/ -c^a + 2Ra 
a 


-1 








-1 












G 








R^ 



a^ 







i?2 sin(^)^ 



a^ 



:i.94) 



/ 



G-^ 






-1 

-c^a + 2Ra 








-i - 


a 



rv2 





u 





i?2 




a" 








- 





R^sm{ey I 



:i.95) 
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Newnian-Unti coordinates (t, R, 6, 



Proof of 1.6.3. Differentiation of the coordinate transformation (1.87) gives 



d \ s\n{9) cos[(f)) 



dy^ = (c^{t) + Rsin{9) cos(0)- 
V ( 



dr 



dR 



— [agsm{9) cos(0) — a cos{9) cos{(j))) d9 H (^a^ sin{9) cos{(f>) 



dy^ = (c'^{t) + Rsm{9) sin(0) 



a_^_ M9)M^) ^^ 



i{9)cos{(j)))d(t), 



— {^ag sm{9) sin{(j)) ~ acos{9) sin(0))(i0 H (^a^ sin{9) sin{(t)) — a sm{9) cos{(f>)'^ c 



a 



cos{9) 



R 



R 



dy^ = (c^{T) + Rcos{9)^)dT-^^°^^dR+^(aeCos{9). + asm{9))d9+^cos{9)a^ 



(1.96) 



where 



a 



da 
97' 



ae 



9q; 



da 



:i.97) 



Substitution of (1.96) into (1.1) yields (1.92). The dual metric (1.93) follows from 

D 



(1.95). 



Lemma 1.6.4. The 1-forms X,V E rA-'^(A^\C) are given by 



X = -RdT 
~ Ra — c^a 



dr — dR 



a 



:i.98) 

^1.99) 



Proof of 1.6.4- 



= -RdT 



V 



. d ^^a 
g{—+X-, 
OT a 

Ra — c^a 



dr — dR 



a 



D 
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Newnian-Unti coordinates (t, R, 6, 



Lemma 1.6.5. 



A ={g{A, V) + R^dr dR + R — —^ -d9 + R — ^^ ^( 



(1.100) 
(1.101) 



Proof of 1.6. 5. 



dr 



A 



-C\T)dy'' + C\T)dy^ + C^{T)dy^ + C\T)dy' 



dr 



dy" 



-C\r)dy^ + C\T)dy^ + C\r)dy^ + c\T)dy^ 



result follows on substitution of (1.96). 



D 



Lemma 1.6.6. 



' — " d ' — " d (yd 

dT = -— and dR = -— + {c'-2R-)— (1.102) 

oR or ^ a' oR 



Proof of 1.6.6. Follows from (B.29) and (1.93) 



D 



Lemma 1.6.7. 



A={R--c 

a 



ad ad aae — adg d da^ — ad^ d 



+ -^ + 



+ 



adR adr R dO Rsm^{e) 



(1.103) 



Proof of 1.6.7. follows from (B.30), (1.93) and (1.101) 



D 



Lemma 1.6.8. 



giX,A) = -R 



a 
a 



(1.104) 



Proof of 1.6.8. Follows by substitution of (1.90) and (1.103) into (1.92). D 
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The Lienard-Wiechert field 



Lemma 1.6.9. 



*1 = ;r^dT A dR A de A d(f) (1.105) 



Proof of 1.6.9. Follows from (1.5) and (1.87). D 



In appendix [E] we present a different coordinate system which we have called 
adapted N-U coordinates. They will be used in Part II of this thesis. 

1.7 The Lienard-Wiechert field 

The Lienard-Wiechert potential is the solution to the Maxwell-Lorentz equations 
when the source J G FA'^A^ is given by J, the current 3-form for a point charge 



moving arbitrarily in free space (1.47). In this section the Lienard-Wiechert poten- 
tial and associated fields are given in term of the null geometry formalism developed 
in the proceeding section. We use the notation A for the Lienard-Wiechert 1-form 



potential as a special case for A G FA^A^. It is the solution to (1.10) given the 
source J. 

Definition 1.7.1. The Lienard-Wiechert Potential of the point charge at point 
X E {Ai\C) is given by 

A|, . mM\C), AU = ^^. (1.106) 

We associate with A the 1-form distribution A-^ defined by its action on test 3-form 
^ e FqA^A^ by 

A^eTdA^M, A^[ip]= f ipAA. (1.107) 

Jm 

Lemma 1.7.2. The electromagnetic 2-form attained by substituting A = A in 
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The Lienard-Wiechert field 



(1.10) will he called the Lienard-Wiechert 2-form and is given by 



F e TA\M\C), 



q g{X,V)X ^A- g{X,A)X ^V - c^X ^V 



47reo 



g{x,vf 



[1.108) 



Proof of 1.7.2. 



F = dA 



-d 



V 



47reo \giy,X) 



-d 



47reo \g{X,V) 



V + 



q 



Aner 



AneogiX.V) 
1 



9ix,vy 



dgiX,V) + 



9iX,V) 



d{V) 
d{V) 



;i.l09) 



Substituting (1.84) and (1.76) into (1.109) yields 



q g{X,V)X AA~g{X,A)X AV-c^X AV 



47reo 



9ix,vy 



(1.110) 



D 



Definition 1.7.3. We associate with F the regular 2-forni distribution F^ defined 
by its action on test 2-forni (f G FqA^A^ by 



F^ e TdA^M F^[ 



ipAF. 



(1.111) 



M 



Readers familiar with the 3- vector notation for the Electric and Magnetic Lienard- 
Wiechert fields can look at lemma [6. 1.71 to see how these relate to F. 

Definition 1.7.4. We split the Lienard-Wiechert 2-form F into two terms where 



'R 



q g{X,V)X AA-g{X,A)X AV 



4vreo 



9{x,vy 



:i.ll2) 
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The Lienard-Wiechert field 



will be referred to as the radiation term, and 



q -c^XAV 
47reo g{X,Vf 



(1.113) 



will be referred to as the Coulomb term. 



Lemma 1.7.5. The Lienard-Wiechert potential (1.106) satisfies the Lorentz gauge 
condition 



d^A = 0. 



(1.114) 



Proof of 1.7.5. 
Let 



K, 



47reo' 



(1.115) 



then 



-d-kA = d( 

K \g{X,V) 

1 



giv,x) 



^dg{V,X)Ai<V 



giv,x) 



d{i<V) 



Substituting (1.84) and (1.77) yields 



1 



.— .7^- *^ - (f#y^)- *- .-a^<i^ -' 



and using lemma B.2.7 gives 



-d'kA = 

1^ g{y,x)^ 



.l-i'-^^^^^).l + m^.l^O (1.116) 



gix^vy 



gix,vy 



D 
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The Lienard-Wiechert field 



Lemma 1.7.6. Off the worldline the Lienard-Wiechert potential satisfies 



d-kdA = di^F = 0. 



Thus given arbitrary region N C {Ai\C) it follows that 



(^Adi<F = 



N 



(1.117) 



for any test 1-form if G T^h^Ai. 



Proof of 1.7.6. 



From (1.108) 



1 ^ g{X,V)^{XAA)~{g{X,A) + c^)^{XAV) 



*F 



K 



g{x, VY 



Thus 



-.'''^='K-^{XA^)-'y w^n' 



d 



rf * (X A A) 

" g{x,vY ' ^\g{x,vy 
'g{X,A) + c' 



A.iXAA)-^-^^^^^d.iXAV) 



9ix,vy 



d 



9{x,vy 



^{X AV) 



;i.ii8) 



Using the chain rule for differentiation yields 



d 



1 



and d 



2dgiX,V) 
9{X,V)^ 
dg{X,A) 3igiX,A) + c'] 



gx,v^ 

9{X,A) + c^ 
g{X,Vf J g{X,Vf g{X,VY 



dgiX,V) (1.119) 



Substituting l\L78^, (|L82]), ([L84j), ([L85]) and ( |LT19| ) into ( |LT18| ) and using lemma 

D 



B.2.9 yields result. 



Lemma 1.7.7. In Newman- Unti coordinates the Lienard-Wiechert potential A G 
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The Lienard-Wiechert field 



rA-'^(A^\C) and the electromagnetic 2- form F G rA^(A^\C) are given by 



and 



A = 

Fr 

Fc 



q 



c 
R 



4vreo 
q fadg — aag 



)dT+ -sdR), 



Aire 







q c^ 



a^ 



dr A dR. 



drAdO + ^^^^^^drA^ 



(1.120) 
(1.121) 
(1.122) 



It follows from theorems C.2.1 and C.2.2 that the distributional l-form A^, E 
r^iA-'^A^ and distributional 2-form F^ G TdA'^M. are well defined. 

Proof of 1.7.7. 



Equations (1.120) and (1.122) follow by substitution of (1.89) and (1.99) into 



(1.106) and (1.113). Equation (1.121) follows by substitution of (1.89), (1.99), 



(1.101) and (1.6.8) into (1.112). 



D 



Lemma 1.7.8. 



^ Fk = ^ r^4^^rfi? A rf^ - ^^^^^^^^^4^^^rfi? A 



and 



• Fr 



Aneo V a^ sin(6') 
q c^ sin(6') 



a" 



Aireci o? 



dO A d(j). 



(1.123) 



Proof of 1.7. 8. Follows from definition B.2.6 , lemma 1.6.6 and the equations 

D 



(1.105), (1.121) and (1.122). 



Lemma 1.7.9. The distributional Lienard-Wiechert field F G F/jA A^ satisfies 



e^d-kY 



D\ 



iDr 



(1.124) 



for any test l-form if G FqA^A^. 

Proof of 1.7.9. First we consider the form of ip close to the worldline. A general 
test l-form if G FqA^A^ is given in Minkowski coordinates by 



v^ = Vii.y ,y ^y ,y )dy 



(1.125) 
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Now making transformation (1.6.1) to Newman- Unti coordinates, such that 



ip = 0r{T, R, 6, (p)dT + 0r{t, R, 6, (j))dR + (peir, R, 6, (fydO + (^^(r, i?, 6, 



[1.126) 



yields 



0R=YlliT,9,<P), 

ipg=Ye\r,9,<P)R, 

ip^=YhT,9,<P)R, 



(1.127) 



where Y^ is a bounded function of r and the rest of the F/'s are bounded functions 
of r, 6 and d). Thus to zero order in R 



f = Y^dT + YJ^dR + 0{R), 



[1.128) 



and 



dip = — ^c^r Ad6 Tr^dr A 

o9 od) 



^yo QY^ dY'^ 

+ -^dr AdR- -^dR Add- -^dR A rf0 + OiR). 

OT OO 0(j) 

(1.129) 



Now definition CI. 4 yields 



D\ 



rf^F^[^] =*F^[d^] 



A^F 



(1.130) 



M 



We split the integral over M. into a region away from the worldline and a region 
containing the worldline. Let the four dimensional region B C A^ be defined in 
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N-U coordinates by 



{ T,R,9,(I) I re I, 0<R<k, 0<9<n, < 4> < 2n } (1.131) 



where / is the domain of C. The boundary d{Ai\C) = dM is given by 



dM = { T,R,9,(^ I re I, R = k, 0<9<7t, < (f) < 2n }. (1.132) 



We calculate d * J-"^ [(f)] with the assumption that fc — )■ so that the approximation 



(1.128) remains valid 



d(hAi^F 



M 



dcpAi<F+ / dipA^F 



/ d{if Ai.F)+ ipAd^F+ difA^F (1.133) 

Jm\«> Jm\m Jm 



The second term in (1.133) vanishes due to lemma (1.7.6). Consider the third 
term. 



I dip A-kF = I dip A -kFc + / 

Jb Jv, Jm 



Using (1.123) and (1.129) yields 



dip A*F 



q f f dYJlc^sm{9) 



4vreo V Jj^ dr o? 



dr AdRAd9 A 



+ 



+ 



dY^ a^a — ad^ 



_A4 i90 a^ sin(6') 

dY^ sin(6')(ad6) — a^d 



M 
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dr A dR A d9 A d(j) 

dr AdRAd9 A 



a^ 



(1.134) 



;i.l35) 



All three terms vanish under integration with respect to R when A; — > 0, therefore 



the third term in (1.133) vanishes. Finally we consider the first term. We note 
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that R is constant on the boundary and therefore dR = 0. By Stokes' Theorem 



d{(p A T^F) = / (/? A T^F 



(^A*Fc+ / <^A*Fr 



(1.136) 



The second term vanishes because dR = 0. We are left with the first term, 



V9 A*Fc 



4vreo 



r"fr 



c^ sin(6') 



dr Ad9 A 



a^ 



Q 



4vreo J^ 



K°M 



TT />27r ^2 



c sin(6') 



J<l>=0 



a^ 



d9d(p ) dr 



(1.137) 
(1.138) 



"IT l'2lV 2 a\y.(fi\ 

Let 1=1 / — dOdcj), then substituting (|1.88|) we obtain 

=0 J^=o ^ 



■K c2-K 



(? sin(6') 



e=o ^0=0 (-CO + Ci sin(e) cos(0) + C^ sin(^) sin(0) + C^ cos(^))2 



(1.139) 



Let z = e^^ such that 



sin(0) = —(2; - 2; ^), cos(0) = -{z + z ^), 



(i^; 



«2; 



(1.140) 



Substitution yields 



1 = 0^ 



42sin(^)2; 



^(0,1) ((-zC2 + Ci) sin(^);z2 + (2C3 cos(^) - 2C0)2 + {iC^ + C^) sin(e))2 

(1.141) 



-dz A d^ 



where /i(0, 1) represents circle of radius 1 centred at the origin. The quadratic: 



iC + C^) sin{e)z' + {2C^ cos(^) - 2C")^ + {iC + C^) sin(0) = (1.142) 
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has roots at 



sin(^)(-iC2 + Ci) 



C^ cos{9) + C° 



± y'CO - 2C0'C3 cos(^) + C^' cos^{9) - C^' - C^' + cos-'{9)C^' + cos{efC^ 

(1.143) 

Denoting the roots by a;(+),/3(— ) yields, 

f 4zzsin(^) , ,, , , 

|a;| > 1 therefore it hes outside the contour. The residue of / at 2; = /3 is given by 

^ 4.sin(g)(a + /3) 

(-/3 + a)3 ' ^'-'^^^ 

Therefore by the residue theorem (see for example [13]), 

io ((C3cos(^)-C0)2 + ((:72^+Ci')sin2(^))i 
and integration using a computer gives 

47r^ (1.147) 



hence 



<^ A *Fc = - / Y^{T)dT = - / '^.(r)dr = ^ /c*0, (1.14J 

Co Jt=~oo ^0 Jt=-oo ^0 J7 



and comparison with lemma 1.4.2 yields 



eorf*F^[0] =g /"c*0 = J^[0] (1.149) 



D 
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PART I 



The self force and the Schott 



term discrepancy 
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Chapter 2 



Introduction 



In chapter [I] we assign the dimension of time to proper time r so that (1.26) is 



satisfied. We will return to this convention in Part II where we derive the electric 
and magnetic fields in the standard 3-vector notation. In Part I it is convenient 



to assign the dimension of length to r so that (1.27) is satisfied. For details see 
appendix |A| 

2.1 The self force, mass renormalization and the 
equation of motion 

It is a consequence of Maxwell-Lorentz electrodynamics that any source of an elec- 
tromagnetic field will be subject to interaction with that field. The resulting force 
on the source is known as the self force. For a charged particle undergoing inertial 
motion the self force is zero, however for an accelerating charge the force is non- 
zero and tends to act as a damping term [H] . It is well known that an accelerating 
charge loses energy due to the emission of radiation, where the instantaneous loss 
of momentum due to radiation, Prad ^ FTA^, is given by the L armor- Abraham 
formula [13] 

pRAD = T^9iC,C)C. (2.1) 
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The negative of this force is the radiation reaction force, which must be a con- 
tribution to the self force. This has led many authors to use the term radiation 
reaction synonymously with self force, however in the fully relativistic case there is 



an extra term in the self force in addition to the negative of (2.1 ). This additional 
term is known as the Schott term, and has lead to some controversy. The fully 
relativistic self force is given by the Abraham-von Laue vector [16] 



fseu=7^(C-9{C,C)C), (2.2) 



where the Schott term is third order with respect to the worldline. 

The zeroth component of the radiation reaction force is Larmor's equation 
for the rate of radiation. The spatial part is proportional to the negative of the 
Newtonian velocity and may be interpreted as the radiation reaction force of the 
particle. The physical nature of the Schott term has been a topic for debate. Its 
presence leads to two interesting results: i) the self force can vanish even when 
the radiation rate is non-zero, for example in the case of uniform circular motion, 
and ii) the self force can be non-zero even when there is momentarily no radiation 



being emitted. Thus the identification of the whole of (2.2) as a radiation reaction 
force would be misleading. The Schott term is a total derivative, so it does not 
correspond to an irreversible loss of momentum by the particle, but plays an 
important role in the momentum balance between the radiation and the particle 



With the self force given by (2.2 ) the resulting equation of motion for a charged 
particle undergoing arbitrary motion is given by the Abraham-Lorentz-Dirac (ALD) 
equation 

mV^C^ = /L + /self+/ext, (2.3) 

where m is the observed rest mass of the particle, /l G FTA^ is the Lorentz force 
due to the external field, /ext ^ FTA^ is the force due to non-electromagnetic 
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effects R and fscu G FTA^ is given by (2.2). All three forces on the right of 
( |2.3[ ) are vector fields with support on finite closed regions of the worldlinaj thus 
Rohrlich's dynamic asymptotic condition [5], 



\imVcC = 0, (2.4) 



is satisfied. The third order nature of the Schott term has instilled doubts about 
the validity of the ALD equation, since it leads to particular classes of solution 
which are foreign to classical physics. These solutions include preacceleration, 
where a particle may begin to accelerate before a force has been applied, and 
runaway solutions, where a particle may continue to accelerate exponentially even 
for a static force (see |H 1201 121] )• Further doubts about the validity of the ALD 



equation are raised by the fact that there remains to this day no derivation of (2.3) 
which is completely free from ambiguity. The most widely known difficulty is that 
of mass renorniahzation. 

The origin of mass renorniahzation can be found in the early attempts to 
calculate the self force based on extended models for the electron. At the dawn of 
the twentieth century the limitations imposed by quantum physics were unknown 
and it was widely believed the dynamics of an electron could be established by 
supposing a classical model for the particle. The model was based on the idea of a 
macroscopic charged object reduced to the microscopic scale. There is an inherent 
problem with this approach because macroscopic charged objects are stable only 
because of the intermolecular forces binding them together. As an elementary 
particle the electron is necessarily devoid of these forces, thus within such a model 
the particle would have a tendency to blow itself apart due to the mutual repulsion 
of its volume elements. The solution of this difficulty, proposed by Poincare, was to 
postulate the existence of an additional cohesive force which would exactly cancel 



^In general these are not known but could include effects due to gravity or collision with 
neutral particles. It is common to assume /cxt = 0. 

^When looking at the solution of this equation it is often useful to consider the external force 
(EM or non-EM) to be a pulse [T71 HH US] , however this is a mathematical idealization. 



45 



The self force, mass renorniahzation and the equation of motion 



the repulsion. This cohesive force would enable the electron to remain stable, 
however it would by definition have no effect on the motion of the particle and its 
physical nature would remain unknown. 

If we accept the Poincare's hypothesis and assume an extended model for the 
electron, then the self force may be calculated using the Lorentz force law. It is 
possible to calculate the Lorentz force acting on a particular volume element due 
to the rest of the charge distribution. The self force is then given by net force on 
the particle due to the respective Lorentz forces on each of the volume elements. 
In order to calculate this force it is necessary to postulate an additional condition 
on the model, that of rigidity. The most common notion of rigidity is that of Born 
rigidity, where the particle is rigid in its rest frame. In the early 1900s Lorentz 
[22] and Schott[23], amongst others, were able to calculate the resulting force for 
a number of different charge distributions. Non-relativistically, for a Born rigid , 
spherically symmetric charge distribution instantaneously at rest, the calculation 
yields [1] 

n=2 

where /self is a 3-vector, x is the 3-acceleration of the charge and the dot denotes 



differentiation with respect to time. The constant k is defined by (L115) and r 
denotes the radius of the distribution. The constant U is given by, 

U = 1 1 ^M^d'^d',' (2.6) 

where n{x)/q is the normalized charge distribution. In the limit r — )■ 0, i.e. the 
point charge limit, the terms in the summation vanish. The resulting equation of 
motion for r — )■ is given by 



2 2 

'moX^-—qKUx + —qKX_ + ji + j^, (2.7) 
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where mo is the bare mass and /l and /ext are 3-vectors. We notice the first 
term on the right hand side is proportional to the acceleration of the electron. 
This led to the identification of the coefficient rrie = -^qKU as an electromagnetic 
contribution to the observed rest mass of the particle. This enables the term to 



be shifted to the left hand side of (2.7), resulting in the equation of motion 



2 

mx = :^qK'x_ + A + ^ct, (2.8) 



where m = m-o + ^^e is the observed rest mass given by the sum of the electro- 
magnetic mass and bare mass. This is known as the Lorentz-Abraham equation, 



and is the non-relativistic limit of (2.3). The process of shifting the term meX to 
the left hand side is known as mass renorm,alization. In the point charge approach 
mass renorniahzation is still required, however the electromagnetic mass of the 
point particle is found to be infinite. This means the bare mass must be assumed 
to be negatively infinite in order to leave a finite observed mass and a meaningful 
equation of motion. This process of adding two infinite quantities to give a finite 
mass is undesirable and brings into question the validity of the resulting equation 
of motion. 

With the advance of physics since the early twentieth century it is now clear 
that any notion of rigidity is incompatible with special relativity. It is also known 
that electrons and other charged elementary particles exhibit wave particle duality 
and other quantum behavior. This has lead to almost complete abandonment of 
the macroscopic model in favor of other models which do not cling to the idea of 
miniature classical distributions of charge. The simplest such model is that of a 
point charge. However if we adopt the point charge model from the outset it is not 
obvious how to define the self force because the Lienard-Wiechert field is singular 
at the position of the particle. In 1938 Dirac proposed a method by which the self 



•^There have been attempts to eradicate mass norm.alization, for example see |24j . where 
different mathematical techniques are used to cancel the singular terms, however there remains 
no physical justification. The inability of classical physics to consistently treat field divergences 
has lead to further needs for renormalization in quantum field theory. 
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force arises as an integral of the stress-energy-momentum tensor associated with 
the Lienard-Wiechert field. In 1973 Rohrlich writes 



Whatever one may think today of Dirac's reasons in developing a clas- 
sical theory of a point electron, it is by many contemporary views (and 
I completely concur), the correct thing to do: if one does not wish to 
exceed the applicability limits of classical (i.e., non-quantum) physics 
one cannot explore the electron down to distances so short that its 
structure (whatever it might be) would become apparent. Thus for 
the classical physicist the electron is a point charge within his limits of 
observation. 

2.2 The point charge approach and the Schott 
term discrepancy 

Within the point model framework the components of the instantaneous change 
in electromagnetic 4-momentum Pem G FTA^ arise as integrals of the Lienard- 
Wiechert stress 3-forms over a suitable three dimensional domain of spacetime. 
This instantaneous change in 4-momentum is identified as the negative of the 
self force but with an additional singular term which can be discarded by mass 
renormalization. In Dirac's calculation the domain is the side S^ of a narrow tube. 



of spatial radius Rbo, enclosing a section of the worldline C. See FIG. 2.1 The 
displacement vector Y defining the Dirac tube is spacelike, therefore the Lienard- 
Wiechert potential is not naturally given in terms of the Dirac time T£){x). However 
in appendix Id] we show that for small Rd = g{Y, Y), and hence small td — Tr due 



to (D.20), it can be expressed as the series 



+ (vD{-g{AD, Ad) - -g{nD, Ad^ - -g{nD, Ad)) - -Ad - -g{nD, Ad)Ad^ Rd 
+ 0{Rl), ( [d:29| 
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C{r) 




Figure 2.1: The Dirac tube (left) and Bhabha tube (right) 



where the vector fields Vd, ^d, ^d G rT(A^\C) are defined as 



C\ CS f) 

VDU = C'{rDix))—-, AdU = &{td{x))—- and AdU = C\rD{x))—-, 



(D.4) 



When using a Dirac tube the integration of the stress 3-forms gives for the 
instantaneous EM 4-nionientuni pTl |2S1 [IHl 121] 



EM 



2.. 



q^[i;(C-9{C,C)C) 



lim - 

Kdo-5>o 2Rdq 



-C 



(2.9) 



This is the Abraham-von Laue vector with the additional singular term which 
depends on the shrinking of the Dirac tube onto the worldline. 

An alternative approach, first used by Bhabha|,27j in 1939 , is to integrate the 
Lienard-Wiechert stress forms over the side St of the Bhabha tube with spatial 
radius Rq. The principal advantage of this approach is that the displacement 
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vector X is lightlike and as a result the Lienard-Wiechert potential can written 
explicitly, 



1 V , . 

-Al, = ——- ( P06| 
K 9{V,X) ^ ^ 



It follows that the corresponding stress 3-fornis can also be written explicitly. 
However previous articles which use a Bhabha tube to evaluate the instantaneous 
EM 4-nionientuni give the following expression ^M ETJ EU [281 ED] 

-P^^ = -qjL(C,C)C + hm -^C) = -P^^-lqKC. (2.10) 

Vd Ro-*0 ZKq / O 

This is the radiation reaction force with the additional singular term which depends 
on the shrinking of the Bhabha tube onto the worldline. The Schott term 2qn^ 
is missing from the approaches employing the Bhabha tube. In 2006 Gal'tsov 
and Spirin [10] draw attention to this discrepancy. They claim the Schott term 
should arise directly from the electromagnetic stress-energy-momentum tensor and 
provide a derivation using Dirac geometry in order to show this. However they 



propose the missing term in (2.10) is a consequence of the null geometry used to 
define the Bhabha tube. We show in this thesis that the term may be obtained 
using null geometry providing certain conditions are realized. 

2.3 Regaining the Schott term 

Addition to non-EM momentum 

The standard approach which has been used in articles [211 [271 EB [28] , is to simply 
add the term to the non-EM momentum of the particle. This method will give the 
correct form for the ALD equation, however it is not physically justified since the 
self force is by nature an electromagnetic effect. 
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We suppose a balance of momentum 



PpART + PeM — /( 



cxt 



(2.11) 



where total momentum has been separated into electromagnetic contribution Pem 
and non-electromagnetic contribution Ppart, and P = ^qP- All the external 
forces acting on the particle are denoted by /ext- A suitable choice for the non- 
electromagnetic momentum Ppart has to be made. Most external forces /ext, 
including the Lorentz force, are orthogonal to C: 



g{U,C) = Q. 



(2.12) 



For such an external force, if (2.9) is obtained then a natural choice for Ppart is 



-PpART = moC. 



(2.13) 



This is the correct term for the 4-momentum of a particle if its spin has been 



neglected. Combining (2.9), (2.11) and (2.13) gives 



TTIqC = /, 



ext Pem 



/ext - qti { MC, C)C + C- lim 



^Do-^o 2Pp)o 



-C 



(2.14) 



Thus assuming the observed rest mass m to be given by 



m = mo + lim 



qK 



■Rdo-^o 2Rj^q 



C 



(2.15) 



we satisfy the orthogonality condition (1.28). By contrast, if (2.10) is obtained 
one cannot set 



moC = /ext - P 



EM and m = mo + lim t— -C 

Ro^o 2Po 



(2.16) 



and satisfy (1.28). This has lead some authors [211 [271 ED EH] to add an ad hoc 
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term to the non-electromagnetic contribution to the force. 

PpART = moC + lqKC. (2.17) 

This ad hoc term will ensure the orthogonality condition is satisfied and hence 
compensate for the missing Schott term. 

Regaining the term by careful analysis of limits 

We will show that the calculation of the self force using null geometry requires 



three limits to be taken (see figure 2.2), the shrinking of the Bhabha tube Sx onto 
the worldline C i.e. Rq — )■ 0, and the bringing together of the lightlike caps Si 
and S2 onto the lightlike cone with vertex C{tq) i.e. ti — ;■ tq r2 — )■ tq, where tq 



is the proper time at which we wish to evaluate the self force (see FIG, 2.1). We 
therefore have the freedom to choose the order of these limits. We choose to let 
the three limits take place simultaneously, subject to the constraint that 



lim 



T2— S>TO 



n +T2- 2To 

AR 



A = .%-t°( " ;; " ) (2.18) 



where A G M is finite. This gives the self force as 

Ui=-q^^(lg{C,C)C + \C+ lim ^C) (2.19) 

V i?o->-0 ZRq / 

which is in agreement with f^^^i ^^ ^ ~ ~\-i hence the Schott term arises by direct 
integration of the stress forms using null geometry. 
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n = T2 





TO = n= T-2 




i?n = 



/ 



Figure 2.2: Three independent limits are required. The converging of the caps 
Ti — )■ r2, the movement of the apex of the squashed tube to the point where the 
self force will be evaluated r2 — )■ tq, and the shrinking of the radius Rq ^ 0. 
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Defining the self force for a point 
charge 



In this chapter we formally define the Dirac and Bhabha tubes. We also present 
the definition of the self force based on conservation of four momentum within 
a Bhabha tube. For this definition we take the limit as the tube approaches an 
arbitrary point on the worldline. 

3.1 The Dirac and Bhabha tubes 

Definition 3.1.1. Consider the region A^ = N\C where A^ C A^ is a local neigh- 
borhood of the worldline. Suppose the two continuous maps 

t' : N -^R, x^ r'{x) (3.1) 

i?' : iV ^ M+, x^ R'{x), (3.2) 

are well defined for all x & N. Here IR+ denotes the positive real numbers and we 
shall call R'{x) the displacement of x from C{t'{x)). Furthermore for A G M"^ let 

r'fxfx - C{t\x))) + C{t'{x))) = r\x), (3.3) 

and R'(x(x - C{t'{x))) + C{t'{x))) = XR'{x). 
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These relations ensure that for an arbitrary point tq on the worldhne with C(ro) G 

N 

hm t'{x) = To and hm R'{x) = (3-4) 

x-^C{to) x->-C{to) 

Since N is open there exist values Tmin, Tmax, -Rmax G N such that the 4-region 

S = I X Tmin < t'{x) < Tmax, < R' (x) < /?[, }, (3.5) 

where R'q < -Rmax, is well defined . 

Definition 3.1.2. The 3-boundary of this region T = (9S is a known as a world- 
tube and is defined by T = S'^ U S2 U S'rp where for i = 1,2 

E^ = I X t'{x) = Ti, < R'{x) < R'q }, (3.6) 

where R'q and Tj are constants and Ti = t'{x) for all x E T,'^. 

S't = I X T2<t' < Ti, R' = R'q }. (3.7) 

We call S'j the caps of the worldtube T and they are surfaces of constant r' whose 
boundaries are topological 2-sheres. We call S'rp the side of T and it is a timelike 
surface of constant R' = R'q topologically equivalent to a cylinder. We call r' the 
worldline map associated with T' and R' the displacement map. 
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Figure 3.1: The Dirac Tube 



Lemma 3.1.3. When using Dirac geometry t' = Td and R' = Ro = ^yg(Y,Y). 
The surfaces S? are subregions of the planes of simultaneity according to an ob- 
server comoving at C{Ti). The Dirac tube T/j is defined by S^^ U S2 U Srp where 
for i = 1,2 



S? = {c(r,) + r I giY,C) = 0, giY,Y) < {Rj^of} 
S? = {c(r) + F I giY,C) = g{Y,Y) = (R^o) 



Tl<T<T2 



}■ 



The parameter Rdq > is a measure of the cross-sectional radius of the Dirac 



tube, see figure 3.1. 
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Figure 3.2: The Bhabha Tube 

Lemma 3.1.4. When using null geometry t' = t^ and R' = R = —g{X, V). The 
surfaces Sj are subregions of the forward null cones at C{Ti). The Bhabha tube 
T^ is given by Si U S2 U St where for i = 1,2 



S. = {c(r,)+X I g{X,X) = 0, -g{X,C) < Ro} , 
^t = {c{t)+X I g{X,X)=0, -g{X,C) = Ro, n < r < r.y (3i 



The parameter Rq > is a measure of the cross-sectional radius of the Bhabha 



tube, see figure 3.2. 
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3.2 Conservation of 4- momentum 



Lemma 3.2.1. Consider figure 3.3 The two Bhabha tubes T™ and T°"* given by 



V" =^f u ^'^ u s^, 

and T°^* =Sf * U S°"* U S°"*, 



(3.9) 



have different radii i?"^ and i?°"*. The surfaces T,f^^ and T,2^^ are the differences 
between the caps of the two tubes, 



yidiff y>out\ y^in 

and J:f^ =S°"*\S^ 



(3.10) 



Let the ^-region S enclosed by the two tubes be finite and source free, with boundary 



^(n v-\difE yidiflf |^ y^out yin 



(3.11) 



then the following relation is true 



sfff 



Sk — / Sk 



Vdiff 



Sk — 



yout 



(3.12) 



Proof of 3.2.1. Using Stokes Theorem (B.103) and (1.37) it follows that 



uSk — / Sk 



N 



ON 



Sk — / Sk 



Vdift' 
^2 



yiQut 



Sk — / Sk, 



(3.13) 



D 
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Figure 3.3: Stokes theorem applied to two worldtubes 
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3.3 The instantaneous force at an arbitrary point 
on the worldhne 

Definition 3.3.1. The k-component of 4-niomentuni flux Pj^ G FA^M through 
an arbitrary source-free 3-surface S C A^ is defined by 

Pi?^ = / Sk, (3.14) 



We cannot use definition |3.3.1| to give the flux of momentum through the caps 
Sj of the Bhabha tube because they are not source free; there is a singularity 
at the point intersected by the worldhne. Instead in the following we employ 
Stokes theorem in order to heuristically justify defining the difference in momentum 
between the two caps as an integral over the side St. 

Definition 3.3.2. Consider the Bhabha tube T = S1US2USX with R = Rq where 
Ro is constant. Let ri = Trijli) and T2 = Tr(T,2) with T2 > ti. The instantaneous 
change in 4-momentum at arbitrary proper time tq is defined by 

PK(ro)=.^'^°f^— / SkY (3.15) 



This definition is justified heuristically as follows. Inspired by definition 3.3.1 
we wish to write 

APk = Pif ^^ - Pf ^^ = f Sk (3.16) 

However the integrals Pj^ and Pj^ ^ are both infinite rl Ignoring this, we assert 

PK(ro) = l-Slif ^— (P?^ - Pif^O) (3-^^) 



Inserting (3.16) into (3.17) yields (3.15) 



^There are methods by which these infinities can be avoided, for example Rowe |29j uses 
distribution theory in order to eradicate singularity, and Norton |24j uses a method where the 
zero limit is not used. 
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Definition 3.3.3. The vector Pem(to) G Tq^to)-^ is defined by 



Pem(t-o) = PK(ro)/' 



_d_ 



(3.18) 



where g^'' = g ^{dy^, dy'') and g ^ is the inverse metric on Ai. Since tq is arbitrary 
there is an induced vector field Pem on the curve C. 

Lemma 3.3.4. Let Lienard-Wiechert stress 3-forms Sk G FA'^A^ be given by 



Q _ qR I qC I qR c 
OK — oj^ "t" Oj^ -f- Oj^ 



(3.19) 



where 



Si =1 (^^kFr a ^Fr - Z9^ * Fr A Fr) 
Sk =1 (^SkFc a *Fc - 29k * Fc a Fc) 
^T =1 (^SkFc a *Fr - t9^ * Fc A Fr + ta^Fn A *Fc - ta^ * Fr A Fc) (3.20) 



Then 



R ^ g{A,A)-g{n,Ay 



^ 167r2eoi?4 



qRC 



q^ 



Sir^eoR^ 



(nK -kV + VK-*^n — nii-ku + gKa * dx^) , 



(Ak -kn + uk-^ a + g{A, n) (Vk ^^71 + n^^-kV — 2nK -kn)). 



(3.21) 



A'"ote that the factor of c in (1.34) ^■5 absent from (3.20). This is due to our decision 



to assign the dimension of length to proper time. 



Proof of 3.3.4. 



We will show only for for S{|, the other results follow similarly. From (1.112) we 
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have 



^''^^''-g{X,V) 



-,^^AXAA)-^^^^^a^{XAV) 



9iX,Vr 



^ (XkA - AkX) - ^/^'tII (XkV - VkX) 



9iX,VY 



where k is defined by (1.115). Thus 



1 



K,^ 



idK^R A ^Fr 



1 



9{X,Vr 
9{X,A) 



+ 



9{X,Vr 

9{X,A) 

9{X,Vf 

9{X,Ar 

9{X,Vr 



XkA a *(X a A)- AkX a *(X a A) 
XkA a *(X AV)- AkX A *(X A V) 
XkV a *(X a I) - VkX a *(X a A) 
XkV a *(X AV)- VkX A *(X A V) 



Now using lemma B.2.9 yields 



(3.22) 



-i^FRAvtFR 



9iX,Vr 
9{X,A) 



+ 



9iX,V)^ 

9{X,A) 

9{X,V)^ 

9iX,Ay 

9{X,Vf 



XK9iA, A)*X- XK9iA, X)^A- Ak9{X, A) ^ X 

- Xk9{A, X)^V- Ak9{V, X)^X 

- Xk9{V, X)^A- Vk9{A, X)^X 
-Xk^X~Xk9{V,X)^V-Vk9{V,X)^X 



Expanding and cancelling like terms yields 



K 



1 ^ ^ / / . .X 9iX,Af\ Xk^X 
^t9^FKA^Fn={9iA,A) ^^ ' ^ » 



giX, Vr / 9{X, Vy 



(3.23) 
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Now we need to calculate the second term in S|J in (3.20). From (1.112) we have 



K- 



;id^ ^ Fr A Fr 



9iX,V) 

9iX,A) /. 



^^i9j^i<{XAA)AXAA 



+ 



9{X,Vy 
9iX,A) 

9{x,vy 

9{X,A)^ 

9{x, vr 



tQj^^iX AA) AX AV 
iQ^i^iX AV) AX AA 
io^^iX AV)AX AV 



(3.24) 



Now using lemma B.2.10 yields 



i9^*{XAA)AXAA ={AKg{X, A) - XKg{A, A)) * X + X^giX, A)^A- g{X, AfgKa * dy" 

id^*{XAA)AXAV =g{A, X)Vk *X + X^giX, V)^A- g{X, A)g{X, V)gKa * dy'' 

id^^{XAV)AXAA =g{A, X)Xk *V - X^giX, A)^X- g{X, A)g{X, V)gKa * dy'' 

ia^*{XAV)AXAV ={VKg{X, V) + Xk)*X + X^giX, V)^V - g{X, VfgKa * dy"^ 

(3.25) 



Substituting (3.25) into (3.24) and expanding yields 



Thus 



-ia^ vc Fr A Fr 



9iX,Af 

9{x,vy 



9{A,A) 



X}^ -k X 

9{x,vy 



Sk ='^'1 (^^kFr a *Fr - tQ^ ^ Fr A Fr) 



Q 



IGvr^eo 



9{A,A) 



9iX,Af \ Xk^X 
giX,Vr giX,V)^ 



(3.26) 



(3.27) 



D 
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Chapter 4 



The resulting expression 



In this chapter we give the result of carrying out the integration in definition 



(3.3.2). We use a computer to carry out the calculation and make use of the 
Newman-Unti coordinate system. The input code for the MAPLE mathematical 
software can be found in appendix [Fj Here we state the result. 



4.1 Arbitrary co-moving frame 

Setting T = tq + 6 we expand Sk in powers of 6. We adapt the global Lorentz 
frame such that 



y^(C(ro)) = for j = 0,1,2,3 
and C{to) = —, C{To) = a—, c{r^) = hi — 



(4.1) 



where a,V eM. are constants given by 

a = ^g{C{T,),C{T,)), V = dy^iCiro)) (4.2) 



and from (1.28) 



b' = a' 
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Thus expanding C and C we have 



0(5 + To) 
C{6 + To) 



1 b\2\ 9 b\^ d 



2 9y2 



a6 + ts')^ + 0{6') 



2 /9y3 



''°4"'''4^'^4^("-^'''^)|5-^°<*^* 



From (1.62) and (1.89) we have 



(4.3) 



V\(s+To,R,e,<t>) = CiS + To) and A\(^s+To,R,e,<t>) = C{5 + tq) (4.4) 

It is useful to express V and A in mixed coordinates, with the basis vectors in 
terms of the global Lorentz coordinates, but the coefficients expressed in terms of 
the Newman-Unti coordinates. 

The result 



Here we outline the steps taken in the MAPLE code. We begin with the expression 



(1.120) for the Lienard-Wiechert potential A in Newman-Unti coordinates. Taking 



the exterior derivative we obtain the field 2— form F and its Hodge dual -kF. We 
obtain expressions for the four translational Killing vectors t^ in Newman-Unti 



coordinates and using (1.34) we obtain expressions for the four electromagnetic 



stress 3— forms Sk. Substituting the expansions (4.3) into these expressions we 



obtain the integrands, and finally using (4.14) we integrate over Et. The result is 



-/ 

1 



1 



qf^ Jet 



So 
Si 

S2 

s. 



1 A2 _ A2 2 

'4 



Rq 3 

+ \b'^-^ + 0{6l) + 0{6l) 

4 /to 



2 



^ui^^^jii _ Za\S2 - Si) - ^ab%6l - 6^) + 0{6l) + 0{6: 



I), 



^-'^'' ^' ' ^^^' ^'+ la\Sl - Si) + 0{5l) + 0{5l) 



^f R. 



+ -a- 
2 Ro 



(4.5) 
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where 



Sl = Ti- To, 



r2 - To 



and K is given by (1.115). 



Combining (4.5) into a single expression and using (3.15) and (3.18) we obtain 
the following expression for P(ro) G Tc{to)-^ 



-Hro) 



3" ay' 



a"7^ + lim 



1 



Ro'^^o2Ro^^y' ' ll^?,\ 4i?o 



■a^ + . 



Vi, + 0{51)+0{51). 



(4.6) 



Hence from the definition of A (2.18) and (4.1) 



-P(ro) = +l9{C{To),C{ro))C{ro) + lim -^C'(ro) + \C{to) + 0{5l) + 0(52^). 

(4.7) 



The first term in (4.7) is the standard radiation reaction term and the second term 



is the singular term to be renormalized. The third term is proportional to C*(to) 
and therefore may be recognised as the Schott term providing the coefficient is 
well defined in the limit. 

If A is chosen to be finite it follows immediately that all higher order terms in the 
series vanish. This is because Rq^ is the most divergent power of Rq appearing in 
the series. Mathematically we are free to choose A to diverge, in which case higher 
order terms could be made finite. However this would require extra renormalization 
in order to accommodate the A terms and the resulting equation of motion would 
not resemble the Lorentz-Abraham-Dirac equation. 

Choosing A to be finite yields for Pem(t) G '^c(t)M. 



qn 



EM 



g(C,C)C + \C+^m--C 
-R.o-i>0 ZRq 



(4. 
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The value of A may be fixed by satisfying the orthogonahty condition (1.28), 



= -^(Pem, C) = -lg{C, C) + Xg{C, C) = -{1 + X)g{C, C). (4.9) 



qn 

2 



Therefore A = — ^ and the final covariant expression for /self is given by 



/seif=-PEM = i«:((7-^(C',C'))C'- lim -^(7, (4.10) 



which is identical to (2.9) 



4.2 Conclusion 

We have shown that the complete self force may be obtained directly from the 
Lienard-Wiechert stress 3-forms when using the null geometry with the Bhabha 
tube as the domain of integration. This eliminates the need to introduce the extra 



ad hoc term in (2.17). It also proves the reason for the missing term in previous 
calculations is the procedure followed in taking the limits, and not the nature of 
the coordinates used as proposed by Gal'tsov and Spirin jlUj . 

We have seen that a requirement for the term to appear is that the ratio of 
limits A, which describes the way in which the Bhabha tube is collapsed onto the 
worldline, is made finite. This is a natural choice because it demands 5i, 82 and 
Rq to be the same order of magnitude. The specific value A = — | is fixed by the 



orthogonality condition (1.28), however the physical justification for imposing this 



particular geometry on the Bhabha tube is currently unknown. 



Consider figure 2.2 It is easily seen that definition (3.3.2) is equivalent to 



Mro) =iSr\{-^ I SkY (4.11) 

It turns out that the limit ri — )■ T2 may be taken before the other two limits. This 
results in the following lemma. 



67 



Conclusion 



Lemma 4.2.1. 



Pk(to) =Ro^o / i_a_SK 



(4.12) 



where S'^{ti) is the 2-sphere given in Newman-Unti coordinates by 



S^ = Ur,R,e,(p) r = Ti, R = Ro, < < vr, < < 27r| (4.13) 



Proof of 4-S-l- In Newman-Unti coordinates the side Et of the Bhabha tube is 
given by 



Et= |(r,i?,^,(/)) ri <r <r2, R = Rq, 0<9<n, < < 27r|. (4.14) 



It follows from (4.11)that 



lim 



Ti^ro \n — ^2 Jt=ti Js^{t) 



Sk(t, -Ro,6',( 



--Ro^oi lim ( / / Sk{t,Ro,6, 



(4.15) 



Applying theorem (B.4.4) yields result. 



D 



Lemma 4.2.1 shows the important step in the calculation is taking the limits 



Rq ^ and ti = r2 — )■ tq simultaneously. If we use (4.12) instead of (3.3.2) for 



our definition of the self force then we obtain the following in place of (4.6) 



-P(ro) 



.2^2 d 



lim /Ti — To 



•^ °y Ro^o 2Ro °^ n^ro V 2Rq 



Vi, + 0{5l)+0{5l). 



(4.16) 



and the orthogonality condition (1.28) yields the key result 



lim / Ti — To \ lim / 5i 



i?o->ol 



<5i^0 V2/£| 



(4.17) 
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In Dirac's calculation ri = tq and so the Scliott term arises naturally without hav- 
ing to take this limit. However when using Dirac geometry the Lienard-Wiechert 
potential and stress forms have to be expanded in a Taylor series about the re- 



tarded time Tr (see appendix D.2). In this process the relationship (D.20) is used 



which gives a relationship between Rdq and 6r = Tjj — Tr which is similar to (4.17) 
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PART II 
A new approach to the 
reduction of wakefields 
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Chapter 5 



Introduction 



5.1 Collimation and Wakefields in a particle ac- 
celerator 

It is common for accelerators to have bunches of order 10*^ particles or more. For 
example, ALICE, at Daresbury Laboratory, uses bunches with bunch charges of 
20pC to 80pC, which represents approximately L25 x 10^ to 5 x 10^ electrons. As 
the bunch traverses the accelerator some of these particles will be perturbed from 
the ideal orbit or trajectory. This may be due to collective instabilities elsewhere 
in the beam or deflection due to residual gas that could not be removed from the 
vacuum chamber. In addition particles from the wall of the beam pipe can be 
accelerated in a process known as self injection. All of these stray particles will 
form a low density region of charge around the beam which is called the beam halo. 

The presence of a large beam halo is generally undesirable. In colliders the 
halo particles reduce the accuracy of measurements at the interaction region, and 
in medical accelerators they can cause severe consequences as a result of highly 
energetic particles missing the desired target. In order to remove the halo from a 
beam specially designed apparatus called collimating systems are used. 

In general collimating systems incorporate regions where the cross-sectional 
area of the beam pipe is reduced. Collimators are specific sections of the beam 
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pipe which undergo a narrowing in one or both of the transverse dimensions. 
There are many possible configurations depending on the design requirements of 
individual projects. In high energy accelerators the presence of collimators can also 
have an adverse effect on the beam due to so called wakefields. Electromagnetic 
fields due to highly relativistic particle beams can interact with the walls of the 
collimator and induce image charges on the wall. The fields resulting from these 
image charges are known as wakefields and can effect the motion of trailing charges, 
often inducing instabilities and emittance growth. Generally fields caused by large 
scale geometric discontinuities, for example in cavities and collimators, are known 
as geometric wakefields and fields caused by resistivity in the wall are known as 
resistive wakefields. 

5.2 Present approaches to the reduction of wake- 
fields 

There is much interest in methods for reducing the geometric wakefields produced 
by a charged bunch of particles passing through a collimator. The customary 
approach is to reduce the taper angle of the collimator. Early work on the cal- 
culation of wakefields from smoothly tapered structures was pioneered by Yokoya 
[30] . Warnock [3T] and Stupakov [321 133]. More recent investigations by Stupakov, 
Bane and Zagorodnov [311 ESI EH] and Podobedov and Krinsky [33 EB] have also 
looked at the effect of altering the transverse cross section of the collimator. A 
detailed analysis of the numerical and analytic calculation of collimator wakefields, 
including an informative introduction to the topic, may be found in [HH]- All of 
the present methods for minimizing wakefields rely on altering the geometry of the 
collimator. In this thesis we propose a new method where the trajectory of the 
beam is altered. 
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Figure 5.1: Synchrotron radiation 



5.3 The relativistic Lienard-Wiechert field 



A relativistic particle undergoing nonlinear acceleration will generate a radiation 



field primarily in the instantaneous direction of motion (see figure 5.1). This is 
known as synchrotron radiation. For high 7-factors the bulk of the field lies 
inside an angle A0 ~ I/7 where A0 is the angle from the direction of motion. By 
contrast, a relativistic particle with constant velocity will generate no radiation 



field. This is easily seen from the form of Fr in equation (1.112) since when the 



acceleration is zero this term vanishes. It is well known in accelerator physics that 
the Coulomb field Fc generated by a relativistic particle moving with constant 
velocity is flattened towards the plane orthogonal to the direction of motion, and 



is often called a pancake field (see figure 5.2) 



Consider figure |5.3 The magnitude of the Coulombic Fc and radiative Fr 
Lienard-Wiechert fields are plotted as height above the sphere for high 7 . In both 
cases the field is distributed in a narrow spike protruding from the sphere in a 
small angle from the direction of motion. In both cases the bulk of the field lies 
inside an angle A0 ~ I/7 from the direction of motion. 



At first glance the plot of the Coulomb field seems to contradict figure 5.2 
which shows the field flattened in the transverse plane. It is reasonable to ask 
how these two radically different behaviours can be consistent. Consider a particle 
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■> 



■> 



Figure 5.2: The pancake field 




Coulomb field 



Radiation field 



Figure 5.3: The magnitude of the Coulomb and radiative fields for a high 7, given 
as height above the sphere. The bulk of the fields is in the direction of motion. 




Figure 5.4: Showing the communication between a particle and its pancake 
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moving at velocity v along the horizontal line PQ in figure |5.4 Let i? be a point 



in the pancake a distance h from the particle, when the particle is at Q- The last 
point at which the particle could communicate with the point R is at P, a length 
vtg from Q. Here tg is the time it takes for light to travel from P to R and also the 
time for the particle to travel from P to Q. Then ||-P-R|| = ctg and ||-PQ|| = ftg. 
Thus (cts)^ = h'^ + {vt^y. Hence h^ = c\^{l - v^/c^) = c\'^/-f^ so 4 = ih/c 
and ll-PQII = 'yhv/c. Thus a particle needs to have travelled in a straight line for 
a length ||-PQ|| = '~fhv/c in order for a pancake of radius h to develop. Looking at 
the fields which originate at P and arrive at R, they are at an angle approximately 



\\RQ\\ I \PR\ = ^h- This is consistent with figure 5.3 



5.4 Proposal 

We investigate the possibility of reducing wakefields in accelerators by placing 
structures which give rise to geometric wakefields, such as collimators and cavities, 
directly after a bending dipole. We model a beam of charged particles as a one 
dimensional continuum of point charges undergoing the same motion in space, 
but at a different time. In our analysis we envisage a collimator as the source 
of wakefields and we calculate the field strength at the entrance of the collimator 
due to the collective Lienard-Wiechert field of the particles in the beam. We do 
not consider any boundary conditions imposed by the beam pipe or the collimator 
itself. We propose the new method of reduction of wakefields should be used 
parasitically on existing bends so that there is no additional beam disruption due 
to coherent synchrotron radiation wakefields (CSR wakes) or loss or energy due 
to synchrotron radiation (SR). In particular an accelerator which requires the 
following: 

• Short bunches (much shorter bunch length that the aperture of the collima- 
tor). 

• The bending of the bunches, via the use of dipoles. 
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• CoUimation. 

can achieve the colhination for free, i.e. with no additional loss of energy or dis- 
ruption to the bunches from geometric or CSR wakes, by placing the collimator 
just after the bend. 

We have seen that in order for a pancake of radius h to develop the particle 
needs to have travelled in a straight line for a distance 'yhv/c. For highly relativistic 
motion this is approximately 7/1. Our proposed method of reduction of wakefields 
relies on this result. The idea is to bend the beam slightly before it enters the 



collimator (see figure 5.5). Most of the Coulomb field generated by the particle 
before the bend will continue in a straight line. By sufficiently enlarging the beam 
pipe in this direction the wakefield due to this part of the field can be neglected. 
If the distance, Z, of the straight line segment from the terminus of the bend to 
the centre of the collimator is sufficiently small, then the resulting pancake field 
will be too small to reach the sides of the structure. Of course bending the beam 
will generate additional radiation fields, however by judicious choice of geometry 
of the beam these can be minimized. 

Let h denote half the aperture of the collimator and let L represent the spatial 
length of the bunch. The following two scenarios will be considered: 

• Long smooth bunches where L > h and any variation in the density of the 
bunch is over length scales longer than h, 

• Bunches where variation in density is over short length scales less than about 
0.2h. This includes the case of very short bunches where L -C 0.2h. 

These two scenarios are both applicable to present day machines, where the 
bunch length depends upon the specific objectives and engineering considerations 
of individual projects. 

In chapter [6] we show that the coherent electric (magnetic) fields due to a bunch 
modeled as a ID continuum of point particles are given by the convolution of the 
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y 

A 

Field measurement 
point X 




Pathofbea: 



X 



Figure 5.5: Setup for beam trajectory and collimator 

electric (magnetic) field due to a single particle with the charge profile. In chapter 
[7] we carry out a numerical investigation using the mathematical software MAPLE. 
Assuming a Gaussian charge profile we minimize the electric field generated by the 
bunch by calculating the field due to different beam trajectories. Having optimized 
the trajectory we calculate the electric field at a specific point, representing a 
point on the collimator wall, for a selection of different bunch lengths which are 
attainable at present day facilities. Calculating the secondary electromagnetic 
fields generated by the collimator is a boundary value problem, hence calculating 
the full Wakefield kick due to the collimator and a bent beam would require detailed 
knowledge of the geometry and material properties of the beam pipe. This will not 
be undertaken in this thesis. However we will show that the field incident on the 
boundary may be reduced by a factor of 7, and since the wakefields are, to a large 
extent, proportional to the fields at the boundary, the field in the beam pipe will 
automatically be reduced by approximately the same factor. We will find that for 
short bunches, or bunches with large amounts of micro-bunching, it is possible to 
make a significant reduction in wakefields. This is applicable to present day free 
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Table 5.1: Bunch lengths for some modern colliders and FELs 



Collider 


Year of 
Commissioning 


Bunch length [ps] 


SLC, SLAG 

ILC 

CLIC 


1989 
>2015 
>2025 


3 

1 
0.15 


Free Electron Laser 




Min. bunch length [ps] 


FLASH, DESY 
LCLS, SLAG 
XFEL, DESY 


2005 
2009 
2014 


0.05 
0.008 
0.08 



electron lasers, which employ bunch compressors to produce very short bunches, 
for example in LCLS L/c ^ 0.008ps. Assuming a collimator of half aperture 
h = 0.5mm then in this case L = O.OOASh. It turns out that electromagnetic fields 
due to long smooth bunches may not be reduced significantly. In many present day 
colliders the bunches are designed to be long and smooth, however in the future 



short bunch colliders may be desirable (see Table 5.1). 
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The field of a ID continuum of 



point charges 



In this chapter we consider the field generated by a ID continuum of point charges 
on an arbitrary trajectory. The key resuh is that the electric field for the continuum 
is given by the convolution of the electric field for a single point charge with the 
charge profile. This result will be used in the next chapter where we adopt the ID 
continuum as our model for a bunch of particles in an accelerator. 

6.1 The Lienard-Wiechert field in 3- vector nota- 
tion 

Definition 6.1.1. Given a choice of time coordinate t such that ^ is Killing we 

can write A^ = M x Ai, where Al is Euclidean three space. We denote the points 
X G {M\C) and C(r) e M hj 

x={cT,X), and C{t) = {ct,x) (6.1) 
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where T,t eM. and X,x E Ai. The null displacement vector X is given by 



X = {cT — ct,X — x), where t = '~frr{cT, X) 



(6.2) 



where the difference X — a; is a 3-vector at point X G Ai. It follows from the 
definition of r^ that T > t. 

Definition 6.1.2. The spatial displacement between the field point X and the 
emission point x will be denoted by 



r = \\X — x\ 



(6.3) 



where 1 1 . 1 1 is the Euclidean norm. We define the unit 3-vector n G TxAl by 



X — X X — X 



n 



\X -x\ 



n .n = 1 



(6.4) 



where the dot denotes the standard scalar product. 



Lemma 6.1.3. It follows from the null condition that 



r = cT — ct 



(6.5) 



Proof of 6.1.3. 



= g{X, X) =g{{cT -ct,X - x), (cT -ct,X - x)) 

= - {cT - ctf + \\X - xW^ 



Thus 



I cT — ct II = ||X — a;|| = r 



(6.6) 



The result (6.5) follows from noticing T > t. 



D 



It follows trivially from definitions 6.1.1 and 6.1.2 and lemma 6.1.3 that the 
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null 4-vector is given by 



X = r{l,n) 



(6.7) 



Definition 6.1.4. The normalized Newtonian velocity f3 and acceleration a are 
defined by 



Idx I dx df3 1 df3 

c dt cj dr' dt 'J dr 



(6.8) 



Lemma 6.1.5. Thus the A-vectors C,C E Tc(t)M. are given by 



C = c7(l,/3) 
and C' = c7^(a./3)(l,/3) + c7^(0,a) 



(6.9) 



Proof of 6.1. 5. First note that 



^.|(l_^./3)-|._l,l_^.^)-i|,_/5.;3), 



thus since 4:i—l3 ■ /3) = —2a . f3 it follows that 



^^ .,3 



dt 



Ya.f3 



(6.10) 



From (6.1.1) we have C = {ct,x). Thus 



■ dC dC dx. 



Also 



" dC dC d'y .^ ^. o d ,^ _, 



Substituting (6.10) and (6.1.4) yields the result 6.1.5 



D 
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Lemma 6.1.6. The following relations are true 

g{X,C) = rcj{n.f3-l) 
and g{X,C) = rc-f\l3.n-l){a./3) + rc-f^{a.n) (6.11) 



Proof of 6.1. 6. 



^(X,C)=^(r(l,n),c7(l,/3))=rc7(/3.n-l) 



Also 



g{X,C) =(7(r(l,n),7c^(l,/3)) +^(r(l,n),f c(0,a)) 
=rc'^^[{a . (3)(3 .n— {a . (3)) + rc'y'^{a . n). 



D 



Lemma 6.1.7. In chapterUl equations (1.16) and (1.22) we define the electric 
and magnetic 1-forms S,B E FTA^ for a timelike observer curve U . Given a 
coordinate chart {y^ , y^ , y'^ , y^) let U = dyo and let E = Ec + Er where 



Ec = «a„oFc and Er = id^F^- 



(6.12) 



// Ec = Ecady" and Er = ER^dy" for a = 1,2,3 then 



E, 



Ca 



q (n - f3)a 



and E 



Ra 



q {nx {n- f3) X a)a 



Ivreo rc7(l — n . /3)' 



(6.13) 



Proof of 6.1.1. Let Fc = Fc^hdz'' A rf^^ then from ( |1.113[ ) and ( |6Tl| it follows 
that 



1. 



K 



■Fcafe 



-C (XflCft — XfoCa) 

{rc^{n . /3 — 1))3 



(6.14) 



82 



The Lienard-Wiechert field in 3-vector notation 



where k = -r^—. Thus 

47reo 



— Eca — — FcaO 
K K 



— C [XaCo — XoCa) 



(6.15) 



It follows from (6.7) and (6.9) that 



Xn 



Xa = rua, Co = -C7 and Ca = cy^^^ 



thus 



1 c^'^rna - rc^jf3^ 



{n-(3)a 



K " (rc7(n . /3 — 1))3 r'^'y'^il — n . I3Y' 



(6.16) 



Similarly let Fr = Yjiabdz"" A dz\ It follows from (1.112) and (6.11) that 



Ip ^ rc7(n ■ (3 - l)(X,a - X,CJ rc-fu ■ a(X.a - X.Cg) 
K ° {rc'^{n . (3 — l)Y (rc7(n . /3 — 1))^ 

(X„a - Xi,Ca) rc-in . a(X,a - X,Ca) 



(rc7(n . /3 — 1))2 (rc7(n . /3 — 1))' 



Thus 



-Er = ^ Fr , = ^^"^° ~ ^°^"^ '^''^^'' • ^^(-^"^° ~ ^°^") (6 17) 

K "" K "" (rc7(n . /3 — 1))2 (rc7(n . /3 — 1))^ 



It follows from (6.9) that 



Co = -C7^a . /3and Ca = C7^aa + C7^(a . /3)/3a, 



(6.18) 



thus the first term in (6.17) yields 



first term 



-rc7^(a . f3)na + r[c-f'^aa + C7^(a . /3)/3„) 



a„ 



(rc7(n . (3 — 1))^ 

7'(a./3)(/3,-na) 



+ 



rc{n . (3 — ly rc(n./3 — 1)2 



(6.19) 
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and the second term yields 

((— rc7na + rc7/3^)(rc7'^(a . /3)(/3 . n — 1) + rc^'^{a . n))) 



second term = — 



(r'C7(n . f3 — 1))' 



l^iria - /3J(a . /3) _^ (n„ - /3„)(a . n) 



rc(n./3— 1)^ rc(n. /3 — l)^ 

Thus adding the two term yields 



1 {n.f3-l)aa , {ua - Pa){a . n) 



K rc{n . (3 — 1)'^ rc{n . (3 — 1)^ 

The result follows from the rule for triple vector products. D 

Definition 6.1.8. Similarly let B = Be + Br where 



Be = -«a * Fc and Bn = -ig^-k Fr, (6.20) 

c ^ c ^ 



Then if Be = ^cady"" and Br = BRarfy" it can be show that 



Bca = -(n X Ec)a and ^Br^ = -(n x ^r)„ (6.21) 

c c 



6.2 The model of a beam 

We model our bunch of particles as a one dimensional bunch where each particle 
undergoes the same motion in space but at a different time. This bunch is moving 
at a constant speed with relativistic factor 7. Let v label the points in the bunch, 
which will be called body points. The profile of the bunch is given by p(z^)r] 

Definition 6.2.1. Let Xyij) represent the position of body point v at proper time 
r, and for each body point z/ let 

t,(r) = (r + z/)/7. (6.22) 

^ Note that v has the dimension of time. 
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Definition 6.2.2. The retarded time for the body point u corresponding to the 
fields measured at X at laboratory time T is denoted by t{X^T,v). Similarly 
the arrival time at X of the field generated by body point v at proper time r is 
denoted by T(z/, r, X). 
For the z/ = particle we define 

fo(X,T) = f(X,T,0) and fo(r, X) = f (0, r, X). (6.23) 

Lemma 6.2.3. It follows that 

ToiX,f{u,T,X)-uM=T. (6.24) 

and 

f(X,T,z/)=fo(X,T-z//7). (6.25) 

Proof of 6.2.3. The retarded time condition is given by 

cT-ct,{f{X,T,u)) = \\X-x,{r{X,T,u))\\, (6.26) 

and hence 

cT - cf(X, T, z/)/7 - CZ//7 = \\X - x^ (f(X, T, u))\\ . (6.27) 

Thus 



CT{y,T,X) = Ctyij) + ||X - Xy{T)\\ 

= c{T + u)h+\\X-x,{T)\\. (6.28) 
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From (6.27) and (6.28) 



cT = c{r{X, T, v) + z/)/7 + ||X - x,{t{X,T, v)) 
= cf (i/,f(X,T, z/),X). 



(6.29) 



Since T is increasing and the range of f is from — oo to +cxd it follows that T and 



f are inverse to each other, yielding (6.29) and 



f{X,f{u,r,X),iy) = T. 



(6.30) 



Now T(i^, r, X) and f{X, T, u) may be written in terms of To(r, X) and To{X, T). 



From (6.28) 



f(z/,r,X)=fo(r,X) + z//7. 



(6.31) 



From (6.29), (6.30) and (6.23), 



fo(fo(X,T),X)=T 



(6.32) 



and 



fo(X,fo(r,X)) = r. 



(6.33) 



Substituting (6.31) into (6.33) leads to 



fo{X, T(z/, r, X) - Z//7) = r. 



(6.34) 



Substituting r = f{X,T, u) and using (6.29) yields 



f{X, T, z/) = fo(X, T - Z//7). 



(6.35) 



D 
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Statistics for independent identical distributions 

Definition 6.2.4. We assume the z/ for each particle has the identical distributions 

pii'), where 



p{u)du = 1 (6.36) 



That is the probabihty that particle k has displacement u^ is p{i'k)di'- 

Definition 6.2.5. Given a function of Hiyi, . . . , z/tv) of all the random variables 
we define the expectation of H as 

(H) = / dvip{vi) ■■■ dvNp{i^N)H{vi, . . . , z/at) (6.37) 



^Af 



Lemma 6.2.6. For a function which is simply the sum of functions X]fc=i ^('^fe) 
we have 



N 



J2h{iyk))=N{h),p (6.38) 



ifc=i 



where 



{h)^p= / p{v)h{v)dv (6.39) 



is the one particle expectation. 
Proof of 6.2.6. 



^ h{vk) ) = duip{ui) ■■■ dvNp{vN) ^ h{vk) 

= ^ c?z/ip(z/i) ■ ■ ■ / dvNp{vN)h{vk) 
k=l -^ -^ 

= ^ / dvkp{i>k)h{vk) = N dup{u)h{i 

1 1 ^ ^ 
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D 
Lemma 6.2.7. The expectation sum of product of the two functions 

N N 

k=l m=l 

is given by 

{H) = N{hg),^ + (N' - N){h),p{g),^ (6.40) 

This is important since it corresponds to components of the energy, momentum and 
stress of the electromagnetic field. In particular, the energy of the electromagnetic 
field is determined in section\6^ 

Proof of 6.2.7. 

IN N \ N N 

(^) =( ( E ^(^'^)) ( E ^('^™)) ) = E E (hii^M^m)) 

\ k=l m=l I k=l m=l 

N N ^ „ 

= E E / ^'^l/'('^l) ■ ■ ■ / d^Np{vN)h{Vk)g{Vrn) 

= E E / ^'^lP('^l) ■ ■ ■ / dVNp{l'N)h{Vk)g{l^m) 
k=l m=k -^ "^ 

+ ^ ^ / duip{ui) ■■■ duNp{iyN)h{uk)g{iyrn) 

fc=l m^k 

= E / dvkp{i^k)h{vk)g{h'k) + Y^Y^ j dukpiuk) / dlJrnP{l^m)h{Vk)g{l^m) 

k=l-^ k=lm^k-^ "^ 

N 

=N / dvp{v)h{v)g{v) + E E ( / '^^kp{vk)h{vk)j y / du^p{u^)g{ 



^k- 

N 

k=l mj^k 

N 

=N{hg),p + E E (^)ip(^)iP = ^(^^)iP + (^' - ^)(^)ip(^)ip 

fc=l m^k 

Note the structure of the expectation of if, in particular the appearance of N and 

N^ -N. U 
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Lemma 6.2.8. The one particle expectation of a shifted function is given by 

(^(T-7'V))^p = j p^,^{T-T')g{T')dT'. (6.41) 

Proof of 6.2.8. 



((7(T-7-V))^p= I p{u)g{T-^-'u)diy 

lp{liT-T'))g{T')dr 
= JpL.y.iT-T')giT')dr 

where T' = T — 7~^z/, and 

PLab(T) = 7P(7T) (6.42) 

is the charge density as measured in the laboratory frame. D 

6.3 Expectation of electric and magnetic fields 

Definition 6.3.1. For a particle of charge q undergoing arbitrary motion x{t), 
where r is the particle's proper time, the Lienard-Wiechert fields at point X and 
time T are given |10] by 

E{X,T) = e(x - x(r^),/3(r^),a(r^,)) (6.43) 



and 



B{X, T) = B (X - x{TR),(3{Tn),a{TR)'^ . (6.44) 



where E and B are defined in lemma 6.1.7 and definition 6.1.8 respectively. 



For the body point u the Lienard-Wiechert electric and magnetic fields at point 



89 



Expectation of electric and magnetic fields 



X and time T are given by substituting tr = f{X,T, v) into (6.43) 



E{X,T,iy) = E(^X-x{T{X,T,u)),(3{T{X,T,u)),a{f{X,T,u)) 



and likewise for B{X, T, v). 

Let Eo{X,T) be the electric field at point X and time T due to the body 
point 1/ = given by 



Eo{X,T) = e(^X -x{ro{X,T)),f3{fo{X,T)),a{fo{X,T)) 



Using (6.25) it follows 



EiX,T,iy) = EoiX,T-u/^). 



(6.45) 



and 



B{X,T,iy)=Bo{X,T-u/^). 



(6.46) 



Lemma 6.3.2. The total electric and magnetic fields at time T at the point X 
are given by 



N 



ErotiX, T, i/i, . . . , un) = J2 -^(^' ^' 



^k 



k=l 



and 



N 



BTot{X,T,iyi,...,UN) = y^^B{X,T, 



^k 



k=l 



It follows that 



(6.47) 



(6.48) 



{E^,,{X,T,iy,,...,iyN)) = N J p{iy)EoiX,T-u/^)du 



(6.49) 
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and 



BTot{X,T,ui,...,UN)) = N p{u)Bo{X,T-u/j)du. (6.50) 



Proof of 6.3.2. 



The result follows from lemma 6.2.6 and equations (6.45) and (6.46) 



D 



Let total electric field at the point X at time T be given by 



^Tot(X,T) = -(^Tot(X,T,Z/i,...,I/^)) 



(6.51) 



We notice that (6.45) and (6.46) are functions where the dependence on u is simply 



shifted g{T — '-) ^v). Thus by lemma 6.2.8 it follows 



^Tot(X,r)= p{v)E{X,T,v)dv 



p{u)Eo{X,T-uh)du 



jp{j{T-T'))Eo{X,r)dT' 
PL.h{T-r)Eo{X,r)dT', 



where T' = T — 1^/7 , and qphahiT) = q'^pi^T) is the charge density as measured 
in the laboratory frame. Thus the key result is that the total electric field is given 
by the convolution 



^Tot(X,T) = y"pLab(T - T')Eo{X,T')dT' 



(6.52) 



The above can be repeated for the total magnetic field Brp^t{X,T). Clearly 
Eo{X,T') will depend on the energy of the beam 7 and the path of the beam 



XiT . 
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6.4 Expectation of field energy and coherence 

Definition 6.4.1. The energy of the electromagnetic field at time T at the point 
X for the A^ particles is defined as the expectation 

(/.(X,T) = (||^Tot(X,T,z/i,... ,1/^)11' + ||BTot(^,T,z/i,...,z/jv)||'> (6.53) 

Lemma 6.4.2. 

0(X, T) = iV0i,,(X, T) + {N^ - iV)0,oh(X, T) (6.54) 

where the incoherent field is given by 

^UX,T) = {\\EiX,T,u)f +\\BiX,T,u)f)^^ (6.55) 

and the coherent field is given by 

0eoh(X,T) = \\E,UX,T)f+\\B,ts{X,T)f (6.56) 

where the one particle continuous electromagnetic fields are given by 

E,UX,T) = {E{X,T,u)),^ and B,ts{X,T) = {B{X,T,u)),p (6.57) 

I.e. Ects{X,T) and Bcts{X,T) correspond to the electric and magnetic fields due 
to a continuous distributions of charge with distribution given by piv). 



Proof of 6.4.2. Expanding (6.53) we see that this is simply a sum of products 
3 

0(X, T) = ^ (Ei,Tot(^, T, Z/i, . . . , Vn) Ei,Tot(^, T, I/i, . . . , Vn)) 
i 
3 

+ ^ (5i,Tot(^, T, I/i, . . . , Z/7v) 5i,Tot(^, T, Z/i, . . . , Vn)) 
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where -Ej^Tot is the z'th component of -Exot- From (6.40) we have 



i i 

3 3 

+ iV 5^ (B,(X, T, z/f >^p + (iV^ - iV) ^ (S,(X, T, J.)) 



IP 



where Ei{X, T, v) is the i'th component of E{X, T, u) 



D 



We've aheady seen from (6.45) and (6.46) that the electric and magnetic fields 



are simply shifted functions so we can use (6.41 ) to give the coherent and incoherent 
fields in terms of convolutions 

<^i„e(X,T)= /"pLab(T-T')(||^o(X,T')f +||Bo(X,T')f )rfT' (6.58) 



and 



E,ts{X,T) = I pi^^^{T-T')Eo{X,r)dT' 
and Scts(X, T)=f pLab(T - T')Bo(X, r)dT' 



(6.59) 
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Numerical results 



In this chapter we present the results of a numerical investigation carried out with 
the mathematical software MAPLE. The relevant code can be found in appendix 



G In the following we give a brief outline of the calculations involved and state 



the main results. 



7.1 The field at a fixed point X for a single par- 
ticle 



Consider figure 5.5 Half the aperture of the collimator is given by distance h. We 



have seen (figure 5.4) that for high 7, a pancake of radius h can develop only if 
the particle has been travelling in a straight line over a displacement of at least 
jhv/c ~ 7/1. Therefore for our proposal to be effective the distance Z in figure 



5.5| should be less than •yh. Preliminary results show that the optimum value for 
Z is Z < lOh, with the field varying little with lower values, thus in the following 
analysis we fix the field measurement point X = (0, h, lOh) and consider the 
magnitude of the electric field at X due a single particle approaching and passing 
through the collimator. In all calculations we use q = —1.60217 x lO^^^C 

We consider the path constructed from a straight line followed by an arc of 
a circle of radius R followed by another straight line. Observe that this path is 
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unrealistic since it would require large magnets to remove the magnetic leakage. 
Magnetic leakage is defined as the passage of magnetic flux outside the path along 
which it can do useful work. In general in a bending dipole the magnetic leakage 
causes the path of the charge to be slightly smoothed out at the ends of the dipole, 
so that in a real bending magnet the path of the charge would not be precisely the 
arc of a circle. We assume that the smoothing of the path corresponding to real 
dipoles would not significantly change the nature of the result. 

Let O denote the angle of arc. The coordinate system is chosen so that the 
direction of the second straight line is along the z axis and the arc is in the x — z 
plane, finishing at the origin. We refer to this trajectory as the pre-bent trajectory 
in contradistinction with that of a particle approaching from (x, y, z) = (0,0 — oo) 
on a straight line towards the origin, which we refer to as the straight trajectory. 

The pre-bent trajectory is given by x{t) = (a;(r),t/(r), 2;(r)) where 



x{t) = < 



R{cos - 1) + {OR + 7-ur) sin for - oo < r < -RQ/^v 
r( cos(7t;r/i?) - 1 j for - RQ/'jv < r < 
for < r < oo, 



Vv) - i for - oo < r < oo, 



(7.i; 



and z{t) = < 



— i?sin + (0i? + 7t>r) cos for — oo < r < —RQ/'jv 
RsinijvT / R) for - RQ/^v < r < 
'-fVT for < r < oo. 



The straight trajectory is given by 



{x{r),y{r),z{T)) = {0,0, -fvr) for 



oo < r < oo. 



(7.2) 
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Calculating the field at X due to a specific path 

We use a coordinate system {r, r, 6, 0} adapted from the Newman-Unti coordinates 



{T,R,6,(f)}. The coordinate transformation is given by (E.2). Comparison with 



(1.87) yields r = — where R = —g{X,V) is the Newman-Unti radial parameter. 
We require the electric and magnetic fields due to a particle on a given trajectory. 
For fixed field point {X,T) = {Xq,Yq, Zq,Tq) there exist parameters r, 9, and (p 
which satisfy 

cTo = C°(r) + r 

Xo = C^(r) + rsin(^)cos(0) 

yo = C^(r) + rsin(^)sin(0) 

Zo = C^(r) + rcos(^). (7.3) 

Rearranging yields the relations 

r =v/(Xo - Ci)2 + (Fo - C2)2 + (Zo - C^y 



cTo(r) = 


=r + C° 




cos(^) = 


r 




sin(^) = 


v/(Xo-ci)2 + (ro- 


-C2)2 


r 




cos(0) = 


Xo-C^ 




v/(Xo-ci)2 + (yo- 


-C2)2 


sin(0) = 


1^0 -C^' 





(7.4) 



These relations can be substituted into the expressions (E.24) for the radiative 
ER(r, r,6',0) and Coulombic Ec(r, r,6',0) electric fields (or magnetic fields). This 
gives the electric field (magnetic field) as a function of the components 0^,0^,0"^, C^ 
and the coordinates Tq, Xq, Fq, ^o- 

When considering the electric field due to a particle on a specific trajectory 
we need only substitute the correct components for C . For example in order to 
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calculate the electric field for the pre-bent path we consider the three sections 



of the path independently. For each of the three intervals in (7.1) we input the 
trajectory by defining components 

C\r)=x{r) 
C\t) = y{T) 
C^{t = z{t) (7.5) 



where the corresponding values for x{T),y{T) and z{t) are defined in (7.1). See 



appendix G lines 131-147. 



In the Maple code we have written a procedure which will take a selection 



of variable input parameters and output any field as a function of r. See (G.l). 
The variable input parameters are the the components C^,C^,C'^,C^ and a list 
of numerical values for the parameters Xq, Yq, Zq and R,Q, 7 as well as an initial 
value for r. 

Lab time 

The ranges of r for the three different trajectories are obtained by substituting 



the numerical input values for R, B and 7 into the intervals in (7.1). Thus for a 
given set of input variables we are able to plot any desired field component for 
a particular section of the path against r for the range of r appropriate to that 
section. In order to plot the field component against r for the whole path we 
simply display the three plots corresponding to the three sections of the trajectory 
on the same graph. 

The lab time To(r) is a different function of r for each of the three sections. This 



follows from (7.3). For a particular section we may obtain Tq{t) by substituting our 



variable input parameters into (7.3) and thus we may plot any desired field against 



To for that section of path by plotting the field and the time Tq as parametric 



97 



The field of a single particle 



equations in r. To plot the field over the whole range of Tq we simply display all 
three plots on the same graph as before. 

Optimizing the values of R and 

We have control over the variable parameters R and O. In order to establish the 
optimum set of parameters to minimize the field at X we calculate the peak energy 
of the electric field | |-Eo(X, T) 1 1 for a range of T, performing a parameter sweep for 
a selection of values for R and O. We set 7 = 1000, which represents an energy level 
easily obtained in modern accelerators. The procedure used in MAPLE is given in 



G.2| and the results are displayed in figure |7.1[ We have displayed three different 
views of the same graph. The numerical values for the electric field are absent 
because we are interested only in the relative values for the different trajectories. 
The values for 6 and R used in these plots are a selection of the values tested, 
however they are sufficient to show the trend. 

Recall that R determines the curvature of the bend and determines the length 
of the bend. Looking at the first graph we see that in the far corner of the graph, 
where and R are at a minimum, the field is at a maximum. As we approach 
the near region of the graph the magnitude decreases very rapidly with increasing 
0. We interpret this as follows. For a short bend the radial distance from the 
point X to the continuation of the straight section will be small, and thus the 
pancake which developed on the straight section will be strongly encountered at 
point X. For a longer bend this contribution will be greatly reduced due to both 
the increased radial distance of X from the continuation of the straight section 
and the increased longitudinal distance of X from the terminus of the straight 
section. The latter distance is important because once the straight section ends 
and the bend begins the pancake is no longer travelling with the particle and the 
field strength within the pancake is decreasing. In consequence we interpret the 
ridge in the graph where the steep section ends as the cut off where the point X 
no longer encounters a significant field due to the pancake. 
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R 


e 


mill 


500 


1/95 




1000 


1/90 




1500 


1/85 




2000 


1/80 




2500 


1/75 




3000 


1/70 




3500 


1/65 




4000 


1/60 




4500 


1/55 




5000 


1/50 




5500 


1/45 




6000 


1/40 




6500 


1/35 




7000 


1/30 




7500 


1/25 




8000 


1/20 




8500 


1/15 




9000 


1/10 




9500 


1/5 


max 


10000 


1 



Table 7.1: Input values for R and 0. 

In reality we cannot adopt the smallest R and largest G because they are 
impractical in the design considerations of real machines. We chose to restrict the 
trajectory to the values B = 0.13rad and R = 0.5m because the bend is sufficient 
to reduce the field at X while also maintaining a minimal length and intensity in 
order to suppress radiation and CSR wakes. In addition a bend of this size would 
be practical from an engineering perspective. 



The field due to a particle on the optimised pre-bent tra- 
jectory compared with that of a particle on the straight 
trajectory 



Consider the two cases given in figure 7^ in which 7 = 1000, x = 0, y = h and 
z = lOh. In the straight line case the peak field is ~ 75Vm~^ and the majority 
of the field arrives within an interval of 0.015ps. In fact it is easy to show that 
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Figure 7.1: Field strength for different values of R and G. We see clearly that the 
minimum field energy occurs when the R is at its minimum value and G is at its 
maximum value. 
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The coherent field at X due to a bunch 



for a straight line path the peak field increases with 7 and the width decreases 
with 7 leading to the classic pancake. By contrast for the pre-bent case the peak 
field is significantly reduced to only ^ 7.7Vm~^, however the interval over which 
the field arrives is now 0.35ps for the right hand peak, and O.lps for the left hand 
peak. The reason for these two peaks is that the left hand peak is the coulomb 
field due to the first straight line segment, whereas the second peak is due to the 
radiation from the circular part of the beam path. The discontinuity is a result of 
the discontinuity in acceleration for this trajectory. Repeating the calculation with 
higher 7-factors does not significantly change the height or shape of the second 
peak. 



Figure |7.3| shows the cartesian components of the electric field for the particle 
on the pre-bent trajectory. We see that the field is largely in the x and y directions, 
with the peak field in the y direction. By contrast for a particle on the straight 
trajectory the y component dominates with the z component negligible and the x 
component zero. This means that for a straight trajectory the field is primarily 
directed in the transverse direction as expected. The nonzero x component in the 
pre-bent case is due to the incident angle of the initial straight section as well as 
the radiation caused by the bend. 

7.2 The coherent field at X due to a bunch 

Consider a bunch modelled as a one dimensional continuum of point particles with 
a low density halo. The fields generated in the halo will not be considered. The 
one dimensional continuum is a good model for beams with transverse dimension 
significantly smaller than the bunch length. The assumption is made that the 
majority of the bunch charge is contained in the one-dimensional core and only 
the halo is removed by collimation. Within this model, each particle in the core 
undergoes the same motion in space but at a different time and is moving at a 
constant speed with relativistic factor 7. 



Using (6.52) we calculate the field due to a Gaussian particle distribution pLab 
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^0 [ps] 
Pre-bent path 

Figure 7.2: The electric field strength ||^o(^,T)|| at X = (0,/i, lO/i), with h = 
0.5mm, due to a body point following a straight path along the z-axis and a body 



point following the pre-bent path given in (7.1) with = O.lSrad and R = 0.5m. 
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Figure 7.3: Electric field components {Eo)x, (-E'o)j/ and {Eo)z at point X = 
(0, h, lO/i), with h = 0.5mm, due to a body point following a pre-bent path with 
e = O.lSrad and R = 0.5m. 
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Table 7.2: Peak field strength for different sized bunches with h=0.5mni. 



Bunch Length 


Peak Etoi 


(X,T)|| [Vm-i] 


L[h] 


(L/c)[ps 


straight 


pre-bent 


1.8 X 10^ 


3.00 X 10° 


1.97 X 10-1 


1.97 X 10-1 


6.00 X 10-1 


1.00 X 10° 


5.91 X 10-1 


5.89 X 10-1 


9.00 X 10-2 


1.50 X 10-1 


3.93 X 10° 


3.48 X 10° 


4.80 X 10-2 


8.00 X 10-2 


7.33 X 10° 


5.27 X 10° 


3.00 X 10-2 


5.00 X 10-2 


1.16 X lOi 


6.36 X 10° 


4.80 X 10-3 


8.00 X 10-3 


5.12 X lOi 


7.53 X 10° 



for the two cases in figure 7.2 The code for the convolution can be found in 



G.4[ We input the time Tq = t at which the convolution should be centered, the 
bunch length (FWHM of the Gaussian) as upper and lower values of Tq, and the 
number of points A^ over which the samples should be taken. The procedure can 
be summed up in the following steps which are followed in a do loop for j = 1..N. 

• Define pLab := {t, a,b)- > l/{a^{27i)) exp ((-(t - 6)2)/(2a2)) 

• Solve To{tj) = tj for Tj, where tj = t — a — {b — a)/N){j + 1/2) and a and b 
are the upper and lower bounds on the range of Tq respectively. 



Substitute r = Tj into the electric field E(To(r)) to give the field strength at 
time To = tj 



Multiply E(tj) by Phahitj), where pLab is a specific Gaussian defined by in- 
putting FWHM. 

N 

Sum the result over j, sum = \^E(tj)pLab(^j) 

N 

Normalize by dividing sum by / ^PLab(^j) 



Table 7.2 displays the results for a selection of bunch lengths which are attain- 
able in some present day machines. 
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7.2.1 Long bunches 

If the bunch is long and smooth, i.e. longer than the collimator aperture, so that 
there is no significant change in pLab over the width of Eo{X, T'), then ^o(^) T') 
may be crudely regarded as a ^-function and E^oti^, T) is given by 

ETot{X,T) ^ pi^^^iT) J Eo{X,r)dT'. (7.6) 

Integration of Eo{X,T') for the straight and pre-bent trajectories reveals that 

||^Tot(X,T)||^ / ^11 ff (7.7) 

ZTTCoC ||A|| 

This value of -Exot is independent of R and G for all paths where R is large 
compared to L. To see why this is the case consider our one dimensional beam 
of particles as a continuous fiow of charge, similar to a line charge in a wire but 
without the background ions. The fields due to this fiow may be calculated using 
the Biot-Savart law. Since h <^ R the field is dominated by the nearby current 
and hence no variation of R, <d or Z will alter the fields. We find that calculations 



using the Biot-Savart Law agree very closely with equation (7.7), thus providing 
verification of our code. 

7.2.2 Short bunches 



If the beam has bunches of length L < 0.05h then it follows from (6.52) and 



figure 7.2 that a considerable reduction in fields is possible. If pLab has full width 
at half maximum L/c = O.OOSps with corresponding bunch length L = O.OOASh, 
then the peak value for the total electric field in the straight line case is given by 
~ 51.2Vm^^. By contrast, in the pre-bent case the peak value for the total electric 
field is ~ 7.5Vm~^, giving an approximate factor of 7 reduction in field. This is 
approaching the maximal factor of 10 improvement one can achieve with 7 = 1000, 
which occurs when the bunch length is small enough that the convolution gives 
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the peak values for the fields in figure 7.2 With higher energies and shorter bunch 



lengths the radiation peak remains unchanged, whereas the electric field for the 
straight path grows linearly with 7. Thus even greater improvements can be made. 

7.3 Conclusion 

In our analysis we have chosen a specific point X = (0, 0.5mm, 5mm) and mini- 
mized the field at this point. We have shown that the magnitude of the electric 
field due to a single particle can be considerably reduced by altering the path of the 
beam, however the duration for which the field is non-zero is increased. We have 
used this result to show that the coherent field for a short bunch can be reduced 
significantly by bending the beam, with reductions of up to 85% feasible for some 
present day FELs and future colliders. No reduction in the coherent field can be 
made for long smooth bunches. We assume the coherent field will dominate the 



incoherent field because of the A^ Vs A^^ behaviour given in (6.54), however the 
incoherent fields are always present. 

If the field point X is instead displaced in the positive x direction, then a 
significant increase in field strength is observed. This increase results from both 
the Coulomb field from the straight section of the path before the arc and the 
radiation from the circular part of the path. 



Consider figure TA The beam has been pre-bent from the left before passing 
through the origin of the graph, hence while on the bend the direction of motion is 
in the positive x direction. The magnitude of the field is shown as a contour plot. 
We see the magnitude increasing as we pass from the origin along the x-axis and 
then decreasing again after a very dense region. This pattern is what we expect 
to find from the synchrotron radiation emitted from the bend. The darkest parts 
of the graph are where the majority of the synchrotron radiation passed through 
the x-y plane. There are four discrete spots; two very dark spots at approximately 
X = 1.4mm and two slightly less intense spots to the left of these at x = 1.0mm. 
The two dark spots are approximately ten times the magnitude of the other two. 
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Figure 7.4: Contour plot for the maximum fields ||-Eo(X,T)|| in the {x,y) plane 
transverse to beam at z = Z. The plot represents a 4mmx4mm region around 
the beam. The beam has been bent from the left before passing though the origin 
of the graph. Twenty timesteps were taken and the lack of smoothness is due to 
numerical errors. Most of the field is between OVm~^ and lOOVm"^ in magnitude 
however the two black spots represent regions where the field is approximately 
1000 times greater. 

All the remaining field shown in the plot is negligible in comparison to these four 
peaks. It will be necessary to alter the shape of the collimator to avoid the high 
field regions interacting with the material in the collimator. This need not affect 



the efficacy of the collimator to remove the halo, for example see figure 7.5 



One criticism of our work is the fact that we have used the Lienard-Wiechert 
field, which is strictly accurate only for a particle in free space. In the accelerator 
community the following formula is often employed to describe the field of a bunch 
traversing a circular beam pipe 



Er 



2A 
r 



(7i 
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Figure 7.5: Modified collimator in the plane transverse to the path of the beam. 

where A is the longitudinal charge distribution and r is the radial distance from 
the axis. This formula assumes o"^ > r/7, where az is the rms bunch length. For 
typical machines this is the normal regime, for example with r = 1cm, 7 = 1000 
then az > lOfim. In our investigation we have shown that the most substantial 
reductions in wakefield are expected for very small bunch lengths, therefore even 
without the difficulty introduced by the bend this formula would be inappropriate. 
The correct form for the field in a curved beam pipe is a boundary value problem 
depending on the intricate geometry of the beampipe, see for example 



In this investigation we have adopted the rigid beam approximation so that 
the charge profile remains constant throughout the whole trajectory. This approx- 
imation is generally adopted in calculating the field generated by a relativistic 
beam traveling in a straight line, however for a pre-bent trajectory there will be 
CSR wakes and energy loss due to radiation which in general will disrupt the 
charge profile. In order to determine whether there will be an advantage in bend- 
ing the beam before collimation it will be necessary to calculate the net effect of 
bending plus reduced collimation wakes compared with collimation wakes on a 
straight beam trajectory. It is probable that the adverse effects of implementing 
an additional bend will be too severe for there to be any advantage in this new 
approach. However all accelerators, even linacs, already have to bend the beam 
using dipoles in certain places. Therefore it seems natural to place a collimator 
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directly after a bending magnet in order to prevent additional adverse effects. The 
optimum design of the beam path, beam tube and collimator shape, for particular 
machines will require a combination of analytic, numerical and experimental re- 
search. Clearly long tapers will reduce the advantage gained by bending the beam 
since it will give time for the pancake to form. However it may be advantageous 
to use a short taper. 
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Appendix A 



Dimensional Analysis 



The SI base units are given by 

L = metre, m 
T = second, s 
M = kilogram, kg 
A = Ampere, A (A.l) 

It is more convenient to use the dimension of charge Q instead of A, with derived 
unit the Coulomb C = sA. We use square brackets to denote the dimensions of 
the enclosed object. The constants eo,/io where c~^ = eo/io have dimensions 

[H ~ (A.2) 

The electric current 3-form has the dimension of charge [J'] = Q. Here the (3) 

(3) 

denotes the degree of the differential form. Dimensions of the electric and magnetic 



111 



Dimensions in chapter\l\ and Part II 



fields are given by 



(0) 


ML 

M 
QT' 


\£]- "" 
\B]- ^ 


and 
and 


(1) 

\B] 
(1) 


ML^ 
ML 


(0) 


^ J LQT^ 


QT 



(A.3) 



and dimensions of tlie electromagnetic 1-from potential A and 2-forms J-" and -kj^ 
are given by 

[A] = [J^] = [^^] = ^ (A.4) 



Also 



[5k] =[Pk] = ^ and [Pk] = [force] = ^. (A.5) 

(3) (0) J- (0) 1^ 



Base quantities have dimensions 



d 



x"] = [dx"] =L, and t^ = y (A.6) 



dx°- 



1 



L 



such that [(7] = L^ and [(? ^] = ;p-. 

A.l Dimensions in chapter [T] and Part II 

We choose proper time r to have the dimension of time T such that 



[C%t)] = L, [C%r)] = ^, [C%r)] = ^ (A.7) 
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It follows 



[X]=[x«-C»(r)][A 
d 



[V] =[C'^iT)] 
[A] =[C%t)] 



d 



1 

T' 
1 



dx"! T2' 



2^2 



[A] = [g{-,A)] 



(Ai 



and 



[^(V^, V^)] 



"2^2' 



[^(^,V^)] 



L2 



[^(^,^)] 



L2 

2^2 



(A.9) 



A. 2 Dimensions in Part I 



We choose proper time r to have the dimension of time L such that 



[C«(r)] = L, [C^ir)] = 1, [C-{r)] 



L 



(A.IO) 



It follows 



d 



[X]=[x«-C»(r)][— 



dx° 



[V] =[C"^(r)] 
[A] =[C"^(r)] 



d 

dx"- 
d 

dx"' 



1 

L' 

1 

Z2' 



1, [X] = [(?(-, X)]=L2, 

[V] = [9{-,V)] = L, 
[A] = [g{-,A)] = 1 



(A.ll) 



and 



[^(\/,\/)]=l, [g{X,V)]=L, [g{X,A)] = l 



(A.12) 
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Appendix B 
Differential Geometry 

In this appendix we present a brief introduction to the geometric constructs and 
notations encountered in the thesis. This introduction is by no means complete, 
and for a deeper understanding of the subject we refer the reader to the vast 
collection of introductory books on the subject, of which [IJIHS] are good examples. 

B.l Tensor Fields and differential forms 

Vector fields 

Definition B.1.1. Let M be an arbitrary differential manifold of dimension m, 
and X an arbitrary point in M. The space of smooth real valued c°° functions 

over M is denoted J-'(M), 

/ e J^(M) implies / : M ^ M, x^ f{x) (B.l) 

Definition B.l. 2. A (contravariant) vector at a point V|a. is a map from J^(M) 
to M, 

V\,:TiM)^R, f^VUf), (B.2) 



114 



Tensor Fields and differential forms 



which satisfies 

VUf + 9) = VUf) + VUg), 

vuxf) = \vuf), 

and VUfg) = VUf)g + fVUg), (B.3) 

where f,gE ^(M) and A is an arbitrary scalar. 

Definition B.1.3. For every x G M the tangent space to M at point x, written 
T^^M, is the m dimensional vector space whose elements are the vectors at x, i.e. 
^U £ T^^M. The tangent bundle is the 2m dimensional manifold given by the 
set theoretic union of the tangent spaces T^M for all x G M, 

TM = y T^.M. (B.4) 

We call M the base space. 

Definition B.1.4. Given the projection map 

7r:TM^M, V\^ ^ x, (B.5) 

a section of the tangent bundle is a continuous map 

V -.M^ TM, X H^ V\^ 
such that T^{V\^) = x for all x eM. (B.6) 

The map V identifies a vector at a point for each point in the base spacq^ therefore 
when acting on a smooth function it is the map 

V : ^(M) ^ M, / ^ V{f) (B.7) 

^We have used the notation V\x to denote a vector at a point as weU as a vector field evaluated 
at a point 
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with the properties (B.3). A smooth section of the tangent bundle is a vector 



field . The space of vector fields over M is written FTM. 

Definition B.1.5. Given a local coordinate basis y"' on M there exists an induced 
JL 

9y" 



local basis J^ on TM. In terms of this basis a vector field V E TM is given by 



V = V—- (B.8) 



where the Einstein summation convention is used for a = l..m and V^ are smooth 
functions on M. In terms of a different local coordinate basis z" 

V = ¥"— = ¥" — — . (B.9) 

Qya dy'^dz^ ^ ' 

where V"" are the components of V in the y"" coordinate basis and V^dz^ /dy"" are 
the components of V in the z'^ coordinate basis. 

Differential l-forms 

Definition B.1.6. The space dual to the tangent space T^M is called the cotan- 
gent space at x and is denoted by T*M. Elements CU- ^ T*M are covariant vectors 
or covectors at x and satisfy 

CU:T,M^M, V^CUV), (B.IO) 

with the properties 

CWvi) = fCUvi)- (B.ii) 

Definition B.1.7. The cotangent bundle is the 2m dimensional manifold T*M 
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defined as the set theoretic union of the cotangent spaces T*M for all x G M, 

T*M = [j T*.M. (B.12) 

xeM 

Definition B.1.8. Given the projection map 

VTc : T*M ^ M, CU^x, (B.13) 

a section of the cotangent bundle is a continuous map 



such that TTciClx) = X for all x G M. (B.14) 

The map C, identifies a covector at a point for each point in the base space. When 
acting on a vector field it is the map 

C : TTM ^ J^(M), V^CiV) (B.15) 

with the properties 

c{v+w) = av)+aw), 

afV) = K{V). (B.16) 

A smooth section of the tangent bundle is called a covector field or (differen- 
tial) l-form . The space of 1-forms is written rT*M. 

Lemma B.1.9. The duality o/T^-M and T*M demands 

CUVl) = VUCl) (B.17) 



where V\x{C\x) satisfies the reversal of (B.ll) with respect to vectors and covector s 



Definition B.1.10. Given a local coordinate basis y'^ on M there exists an induced 
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local basis dy"' on T*M. In terms of this basis a differential l-form ( G T*M is 
given by 

C = Cady^ (B.18) 

where (a sue smooth functions on M. In terms of a different local coordinate basis 



C = Cady^ = Ca^dzK (B.I9) 



where C" are the components of C in the y"" coordinate basis and Cag^F are the 
components of ( in the z"- coordinate basis. 



Tensor fields 

Definition B.1.11. The degree of an arbitrary tensor will be represented as an 
ordered list s of or more entries. Each entry is either the symbol F (for l-form) 
or V (for vector) e.g. s = [V, F, F, V]. The space of tensors of degree s over M is 
denoted 1^^ M with sections F 1^'^ M. Let the tangent space and the cotangent 
space be denoted 

TM = (gjf^^M and T*M = (g)f^^M (B.20) 

then arbitrary degree tensors are constructed using the tensor product, 

® : (g)"M X (g)*M ^ (g)'"'*^M, (T,S)^T®S (B.21) 

where s and t are ordered lists and [s, t] is simply the concatenation of the two 
lists. For example given vector field V G FTM and 1-forms a, (3 E rT*M we may 
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define a degree [V, F, F] tensor field by 

V^a^Per(S)^^'^'^^M. (B.22) 

The tensor product satisfies 



T®(S®R)=(T®S)®R 
(Ti + T2) ® S =Ti ® S + T2 ® S (B.23) 

and 

/(T®S)=(/T)®S = T®(/S) (B.24) 

for T,Ti,T2 G (gj'M, S G (8)*M, R G (8)"M and / G J^(M). The dual space of 
(3)'^M is denoted (3)''M, where s is the list obtained by interchanging the symbols 
F and V in s. The total contraction of elements in (^''M with elements in (^*M 
is written 

(g)"M X (g)'M -^ -F(M), (T, R) ^ T : R (B.25) 

where T G (^*M and R G (^''M. It is defined inductively via 

V ■.a = a:V = a{V) where a G (g)'^lM and V G (g)'^', (B.26) 
and extended to arbitrary tensors by 

(S®T) : (R®U) = (S:R)(T:U) (B.27) 
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The metric 

Definition B.1.12. Of special importance is tlie symmetric, non-degenerate de- 
gree [F, F] tensor field g E ^ ' M with the properties 

giV,W) = giW,V), 
and g{V, W) = for all V ^O^W = 0. (B.28) 

This tensor is called the metric. It provides an isomorphism between the covari- 
ant and contravariant vector fields. The metric dual V of vector field V is the 
differential 1-form given by 

V = giV,-), (B.29) 

where (— ) denotes an empty argument. There exists a symmetric, non-degenerate 
degree [V, V] tensor field g^^ E (^ ' M which satisfies 

g-\V,W)=g{V,W) and V = g-\V,-). (B.30) 

Given the local coordinate basis y^ on M the metric g is given by 

g = gabdy'' dy'' (B.31) 

where the functions gab are determined by 



*'=K|^-<|') (B'^^' 



Differential p-forms 

Definition B.1.13. An important subspace of (3) M is the space of totally anti- 
symmetric degree [Fi, ..., Fp] tensors denoted A^M. The space of smooth sections 
of A^M is denoted FA^M and elements \E' G FA^M are called (differential) p- 
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forms. For example for z/, cj G r(^' ' M the totally antisymmetric part of the 
degree [F, F] tensor field z/ (g5 cj is the difference 



-(z/Ow-wOz/) = ^ (B.33) 



since reversing the positions of u and u yields 



-{u <S) u - u <S) u) = -"^ , (B.34) 



threfore \E' G FA^M is a differential 2-form. The space of differential 0-forms 
FA°M is defined as the space of smooth functions over M, i.e. J-'(M) = FA°M 
and the space of differential 1-forms FT*M = F(3) M is now also written as as 
FA^M. 

Higher degree forms are obtained using the exterior or wedge product. 
The wedge product of a G FA^M and /3 G FA'^M is the map 

AiFA^MxFA^M -^ FA^+'^M, a,/3 ^ a A /3, (B.35) 

where the (p + g)-form a A /3 is the totally antisymmetric part of the tensor field 
a ^ (3. The wedge product satisfies 

a A (/3 A7) =(a A/3) A 7, 
(ai + a2)Af3 =ai A /3 + ^2 A /3, (B.36) 



and 



f{aAf3)={fa)Af3 = aA{ff3), (B.37) 



and 



a A (3 = {-1}^" (3 A a (B.38) 
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for a, ai, a2 e A^M, /3 G A«M, 7 G A'^M and / G J'(M). It follows any arbitrary 
p-form can be reduced to a linear superposition of wedge products of p differential 
1-forms. Given local coordinate basis y"" on M, 



a 



G TA^M, a = aa,a,..a,dy''' A rfy"^ A .. A dy''" . (B.39) 



B.2 Differential operators 

Definition B.2.1. For the following definitions it is useful to define the map 

7] : A^M -^ A^'M, a ^ r]{a), 

where r]{a) = (-l)^a for all a G A^M. (B.40) 

Exterior derivative 

Definition B.2. 2. The exterior derivative of a differential form is the map 

d : TA^M -^ TA^+^M, a ^ da, (B.41) 

where 

df{V) = V{f), 
d{a A l3) = da A l3 + r]{a) A djS, 

d{da) = (B.42) 

for all / G rA^M, V G TTM, a G rA^M and /3 G TA^M 

Lemma B.2. 3. For 1-form v G TA^M and vector fields V,W e TTM the follow- 
ing relation holds 

dviV, W) =V{v{W)) - W{v{y)) - v{\V, W]), (B.43) 
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where [,] denotes the Lie Bracket. 

Proof of B.2.3. It is sufficient to prove for the 1 — form v = fdg, where f,g& 
J^(M). In this case du = df A dg and thus 

diy{V,W) =df Adg{V,W), 

= l{df{V)dg{W)-dg{V)df{W)), 

= l{Vif)Wig)-Vig)Wif)). (B.44) 

Here du^V, W) is the action of the 2— form du on the ordered pair of vector fields 
{V,W). Also 

V(iy{W)) - W{u{V)) - u{[V, W]) 

=V{fdg{W)) - W{fdg{V)) - fdg{[V, W]) 
=VifiWig))) - WifVig)) - f[V, W]ig) 
=Vif)Wig) + fViWig)) - Wif)Vig) 

-fW{Vig))-f[V,W]ig) 
=Vif)Wig)-Wif)Vig) 
=2diy{V, W) (B.45) 



D 



Interior contraction 

Definition B.2.4. The interior contraction of a p-form with respect to a vector 
field is defined by 

i : TTM X rA^M -> rA^-^M, V,a^ iya (B.46) 
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where 

iv{a A (3) = iyo. A /3 + 77(a) A ivP 

iyiya = (B.47) 

for all V e TTM, u e TA^M, a G rA^M and ^ G TA^M. 

Hodge Dual 

Definition B.2.5. Introduce a gf-orthonormal frame e" such that 

g = Qabdy'' ® dy^ = r/,fee" ® e'' (B.48) 

where 77^6 = ±1 for a = 6 and rjab = for a 7^ 6. Then 

^1gA"M, *l = e° Ae^ a.. Ae'""\ (B.49) 

is the volume form on M. 

Definition B.2.6. The Hodge dual is the map 

^ : A^M -^ A^-PM, a^M (B.50) 

where 

^(1) = .^1 

■kP = iy-k 1, 

■k{a A u) = i-p -ka, 

i,(fa) = fi.a. (B.51) 



for all u G rA^M, a G rA^M and f3 G TA^M. Properties ( |B.39[ ) and ( |B.2.6[ ) 
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define the action of tlie Hodge dual on any arbitrary p-forni. 



Lemma B.2.7. Given two vector fields V,W & FTM the following relation is true 



VA^W = g{V, W) ^ 1 



(B.52) 



Proof of B.2.7. 



V Ai^W =V Aiw^l 



dzb 



(B.53) 



Using tlie alternating Leibniz rule (B.47) yields 



= ig^b^dz"- A *1) = iQ^.bdz"' A ^1 — dz"" A ig^t t^ 1 
and thus dz" Aig^t i< I = 5^ i<l. (B.54) 



Substituting (B.54) into (B.53) yields 



VAi^W =VbW'' ^ 1 
=g{V,W)^l 



(B.55) 



Lemma B.2.8. Given v E FA^M and a E FA^M the following is true 



D 



u A-ka = — -kU-pa 



-rvV-] 



(B.56) 



where a^ = ri{a). 



Proof of B.2.8. By lemma B.2.7 it is clearly true for deg(a) = 1. Assume true for 
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deg{a) = p, then for u G FA-'^M 



•kip(a Aid) = -k (i-pa Au) + •k(a'^ A i-pu) 
= -k {ipa Aoj) + g(u, u) i< a^ 
= -k (ipa Au) + iz,{u A -ka^) + u Ai^i^ a^ 



Evaluating the first two terms yields 



Thus 



k ip(a Au) = — vAiz,ka'^ = vA -kia A iS) 
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^{ipa Aoj) + i^{u A i^a^) ='k {u A ipa^) — i^ii^kiipa)) 

= -k {oj A ipa"^) — k:(ipa A u) 
= -k{uj A ipo^) — *(a; A iya"^) 
=0 (B.58) 



(B.59) 



D 



Lemma B.2.9. Given 1-forms v,u,a E FA^M the following is true 

V A ^{uj A a) =g{y-, a)^uj — g(u, u) k^a (B.60) 



Differential operators 



Proof of B.2.9. By lemma B.2.8 we have 



v A -ki^oj A a) = — -kiiyiuj A a) 



- — -*^{iuOO A a — a; A i-pa) 
= -k{g{i^,a)u) --k{g{u,u)a) 
-g{u, 5) * a; — g{v, cD) * a 



D 



Lemma B.2.10. Given one forms a,P,'y,i',uj G A^M the following is true 

2,:^ * (a A /3) A 7 A z/ = {l3{^)u{u) - u{^)l3{u)) ^ a 

+ (cj(7)a(z/) - a{^)u{u)) -k j3 

+ (a(7)/3(?)-/3(7)«(^))*c^ (B.6i; 



Proof of B.2.10. 



i[3-k{aAf3)A'yAh'=-k{aAf3Auj)A'yAu 

= --i A-k{aA(3 Auj) Av (B.62) 



Using (B.2.S) yields 



7 A k{a A 13 Alu) = - k[i7y{-a A /3 A w)) 
= k [i:y{a A l3 Au)) 
= k {a{^)l3 Au - l3{^)a Au + u{^)a A /3) 



(B.63) 



Thus substituting (B.6S) into (B.62) yields 



i^k{aA(3) Aj Au = -u A *(a(7)/3 Au- (3{'y)a Au + u{'y)a A (3) (B.64) 
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Now using (B.2.8) again yields 



iiD * (a A /3) A 7 A z/ = ^ f is {a{-f)/3 A w) - i^ {P{l)a Au) +i^ (^(7)" A /3) 

+ (w(7)a(z/) - a{^)u{u)) -k (3 



D 

Lemma B.2.11. For two forms a, /3 G FA^M of the same degree the following is 
true 



a A-kfi = (3 /\ -ka. 



(B.65) 



Proof of B.2.11. By (B.2.7) it is clearly true for deg(a;) = deg(/3) = 1. Assume true 



for deg(a) = deg(/3) = p, then for for a E TAp+^M, (3 e rA^M and u E TA^M we 
have 



a A ^{P A z/) =a Aiu^ (3 

=iy{a^ A 7^/3) - ijjo;'' A ^{3 
= — P A ki-j^a^ 
=13 Av Ai^a 



(B.66) 



Lie Derivative 



Definition B.2.12. The Lie derivative is the the map 



D 



L : FTM X r(g)WM -^ F^^M, V, T f^ r^T, (B.67) 
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which is additive hnear in both arguments 



Cv{T + S) = CvT + CvS, 
Cv+wT = CvT + CwT, (B.68) 

has the properties 

Cvf = V{f), 

CvW=[V,W], (B.69) 

and obeys the Leibniz rule for tensor products, wedge products and contractions 

£y(T®S) =£v'T®S + T®£yS, (B.70) 

Cv{aAP)=CvaAP + aACvP, (B.71) 

Cv{a{W)) = Cva{W) + a{CvW). (B.72) 

Lemma B. 2. 13. Cartan's formula 

Ly = div + iyd (B.73) 



Proof of B. 2. 13. Trivial for O-forms. First prove for 1-form v G A^M. From (B.72) 



and (B.43) 



Cvy{W) =Cv{y{W)) - u{CvW) 

=V{p{W))-p{\V,W]) 

=2dv{y,w) + w{v{y)) 

=ivdv{W) + d{v{y)){W) 

=ivdv{W) + diyviW) 

={ivdu + divy){W) (B.74) 
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Now assume true for p-form a G FA^'M, show true for {p + l)-forni a A u. 

{iyd + div){ci A z/) =iv{da A z/ + (— l)^a A dv) + (i(zya A z/ + (— l)^a; A iyz^) 

=iv{da A z/) + (-1)^zv/(q; A (iz/) + (i(iyQ; A z/) + (-l)Prf(a; A iyz/) 
=ivda A z/ + ((-1)^+^ + (-l)^)(ia A iyi^ + diya A z/ 

+ ((-1)^"^ + {-lY)iva Adu+ (-l)^^(a A iydu + a A dzy'^) 
= {ivd + (izi/)ct A z^ + a A (iyd + diy)^ 
=Cya A z/ + a A £yz/ 
=Cy{a A z/) 



Thus by induction true for all p-forms. 



D 



Lemma B.2.14. The Lie derivative commutes with the exterior derivative 



Cyd = dCi 



(B.75) 



Proof of B.2.14- follows trivially from B.73 



D 



Levi-Civita Connection 



An affine connection is a map 



V : TTM X r(g)''lM -^ T^'^'m, ( V, T) h^ VyT 



(B.76) 



which is additive linear in both arguments, 



Vv-(T + S) = VvT + VyS, 

Vy+wT = VyT + VvkT, 



(B.77) 
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and satisfies 

Vy/ = V{f), 

VfvT = fVvT, 

Vv(T ® S) = VyT ® S + T ® VyS, 

Vv(a A/3) = Vya A/3 + a A Vy/3. (B.78) 

The Levi-Civita connection on M is the unique torsion free metric compatible 
affine connection. 

B.3 Pushforwards, pullbacks and curves 
Pushforward map 

Definition B.3.1. Given differentiable manifolds M and N and the smooth map 
: M — > N, X I — > 4>{x), the pushforward of a vector at a point V\x G T^-M 
with respect to is the map 

0, : T,M ^ T<^(,)N, V|, ^ 0,V^U, (B.79) 

where 

<P*vuf) = vufo^), 

and MVUW\,) = (P.V\,(P,W\, (B.80) 

for l^|a;,W^U ^ TjjM and / G FA^N. The pushforward of a vector at a point is 



naturally extended using (B.6) to obtain the pushforward of a vector field 



!), : TTM -^ TTN, V ^ 0,1/. (B.81) 
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Let M, N and O be differential manifolds and let 0=,, : T^jM — > T^(2:)N and 
ij^ : T^N — > T<^(a.)0, then 

Lemma B.3.2. The composition of the pushforwards is the pushforward of the 
composition 



i)^ o (j)^ = [if) o (p). 



(B.82) 



Proof of B.3.2. 






u 



Lemma B.3.3. Let {x^.,x'^ , ...^x™") he a coordinate basis o/M*^ and {y^,y'^, ...,y^) 
a coordinate basis o/M^. Given a diffeomorphism where 



)M , wdN 



X = a;"(a;) = {x^{x), ..,x"'(x)) 



x)=y'{<P{x\x),..,x"^{x))), 



then 



d 



' dx"- 



d 



X dx"- dy 



a flo,b^ 



b „.b 



where 4> = y o 



(B.83) 



Proof of B.3.3. 



d 



' dx"- 



(/) 



d 



dx" 
d 



dx" 
d 

dx" 



{fo<P) 

(/o/°0) 
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where y^{(j){x)) = 0^ by definition. Evaluating using the chain rule yields 



^ ■(/(/)) ''''f 



Qj.a Q^a Qyb 

(90^ d 



(/) 



U 

Fullback map 

Definition B.3.4. For the smooth maps : M — )> N and / : N — )> R the pullback 
of / with respect to is given by the composition: 

07 = ^ o / (B.84) 

Lemma B.3.5. For the smooth maps : M — )■ N, -0 : N — )> O, and / : N — )> M 
the pullback of the composition is the composition of the pullhacks 

{(f) o ipy = (f)* o ip* (B.85) 

Proof of B.3.5. 

= /o0o^ 

= (0O^)*/ 

D 

Definition B.3.6. The pullback of a differential p-form with respect to the smooth 
map : M — > N, x i — > 4>{x) is given by: 

0* : TA^N — > TA^M, a ^^ (p*a (B.86) 
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where 

(f)*{a + /3) = (j)*a + (j)*/3 

0*(aA/3) = 0*aA0*/3 (B.87) 

for all a e TA^N and (3 e TA^N. 

Definition B.3.7. The pullback of a 1-form acting on a vector field is the l-form 
acting on the pushforward of the vector field 

<t>*df{V)=df{(P,V), (B.88) 

for all / G TA^N and V E TTM. 

Lemma B.3.8. the pullback commutes with the exterior derivative 

d(j)*a = (j)*da (B.89) 

for all a E TA^M. 

Proof of B.3.8. First show for / E TA^M 

<p*df{y) = df{<p,v) 
= <P*vif) 

= \/(/o0) 

= V{<j>*f) 
= d<l>*f{V) 



134 



Pushforwards, puUbacks and curves 



Now show for u = gdf e TA^M where gje ThPM 

(P*d{gdf)=^*{dgAdf) 

= {(l>*dg) A ifdf) 
= d(j)*g A 0*rf/ 
= d{<P*g<P*df) 
= d(P*{gdf) 

proof for a general 1-form u = v-idx'^ follows by linearity. 
Now, assuming true for a G FA^M 

(f)*d{i' A a) = (f)* [du A a — u A da] 

= (j)*{dv) A (j)*a - (jfv A (f)*da 
= dcffv A (f)*a - 0V A (i0*a 
= d[(/)*uA(t)*a] 
= d(t)*{v A a) 

Hence by induction the relation holds for all p-forms. D 

Lemma B.3.9. The internal contraction of a pullback with respect to a vector 
field V , is equal to the pullback of the internal contraction with respect to the 
pushforward of V 

iy(f)*a = 0*(i0.^a) (B.90) 

Proof of B.3.9. (by induction) 
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Trivial for O-forms. First prove for l-form u = gdf G FA^N 

ivi(l)*9df) = iv{,<P*g<P*df) 
= (t)*giv(j)*df 
= <P*g{(t^*df).V 

Noticing that (f)*df{V) = df{(f)^:V), and remembering that the pullback of a number 
doesn't change the number, ie 

rf/(0.\/) = 0*(rf/.0,\/) 
= 4>*{i,t>,vdf) 

We can therefore write 

'i'v{4>*gdf) = (t)*g(j)*{i<t>,vdf) 

= <P*{.gH,vdf) 

= <t>*h.v{9df)] 

proof for a general l-form u = v-idx'^ follows by linearity. 
Now, assuming the relation holds for aeFA^N, we have 

iy(t)*{u A a) = iv{(p*i^ A (f)*a) 

= iv(f)*i^ A (f)*a — 0*z/ A (f)*a 

= 'f{i<t>,vi^) A (jfa - (jfv A 0*(^</,.\/") 

= 0*(v.v(«^ Aa)) 

Thus we have proved by induction that the relation must hold for all (p+ 1) forms, 
and therefore for any general form /3 G FA'^N. D 
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Curves 

Definition B.3.10. A smooth parameterized curve C{s) on a manifold M is a 
smootli map from an open interval / C M to M, 



C : J ^ M, s ^ C(s) 



(B.91) 



If y"" are local coordinates on M then we use the notation 



ya^Cis)=C^is), 



(B.92) 



thus if C(so) = a; is any point on the image of C then 



y'^ix) = C^iso). 



(B.93) 



The tangent vector to C at x is 



^ U £ -tx-l-Vl, G la; — G* 



d_ 
ds 



so 



(B.94) 



For any / G J^lM) 



G* 



9_ 
9s. 



d 



)L.(^' = ai(^"^) 



So 



9 a/ac« 



so 



(B.95) 



hence 



• d 



9C" 



ds 



d 

^odr' 



(B.96) 



There is an induced vector field C G FTM where (71^: is the tangent vector at 
X for all X G C{s). 
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B.4 Integration of p-forms 

Definition B.4.1. Let a be a diffeomorphism from the submanifold S C M of 
dimension n to tlie differentiable manifold M of dimension m. 

a-.E^M (B.97) 

If y" are local coordinates on M at cr{x) then a* acting on the local basis of 1-forms 
dy"' is given by 

a*dy'' = diy" o a) (B.98) 

and for any / G J-'(M) 

a*ifdy^) = ifoa)diy^oa) (B.99) 

If we define a local coordinate system for S at a; G S by 

z'' = y''oa (B.lOO) 

then 

(T*{fdy° A dy^ A .. A dy"') = (/ o a)dz^ A dz^ A .. A dz"" (B.lOl) 

Definition B.4. 2. If m-form a G FA^M has compact support then so does the 
n-form cr*Q; G rA"S, and 



a = / a*a. (B.102) 

M is 



Theorem B.4. 3. // S is an oriented differential manifold of dimension n, with 
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boundary dT, of dimension (n — 1) then 



da= a, (B.103) 

for all a G rA"~^S with compact support. This theorem is often called the gener- 
alized Stokes' theorem. 

Theorem B.4.4. Let t be a choice of coordinate on a manifold M such that ^ is 
Killing and let t foliate M into surfaces S^ . Then for a G FA^M 



, , a= Caa, (B.104) 



and thus 



a = dt iQ^a (B.105) 

M(ti,to) Jt=to JT,t 

where M(ti,to) is a submanifold ofM. with range oft between to and ti. 
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C.l Definitions 

The space of C°° functions with compact support is called the space of test func- 
tions. We extend this notion to the space of test p-forms. 

Definition C.1.1. Let M be a differential manifold of dimension m. The space 
of test p- forms on M is denoted FqA^M, 

TqA^'M = ye VM'M.\ if has compact support}. (C.l) 

Definition C.l. 2. The space of p-form distributions T^iA^'M. is the vector space 
dual to the space of test {m — j9)-forms FqA^^^M, 

TdA^M X ToA^'-^M -^ R, (^, v?) ^ ^M G M, (C.2) 

which satisfies 

^\^ + ij] = \^^] + ^j], (C.3) 

for \eR,<^,ije ToA'^-^'M and ^ G TdA^M. 

Definition C.1.3. The subspace of F/jA^M comprising piecewise continuous 
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p-forms is the space of regular distributions. The action of a regular p-form 
distribution ip^ on an {m — p)-test form ip is given by the integral 



^^[V^] = / v^A^ (C.4) 

'm 



for any ip G FoA™ ^M and where ip G FA^M is piecewise continuous. We say that 
ip^ is the j9-form distribution associated with the p-form ip. 

Definition C.1.4. The exterior derivative of a p-form distribution is defined 
as: 

d : FbA^M -^ Fz)AP+^M, ^ ^ d^ (C.5) 

and satisfies 

d^[^] = -^[rfv?"] (C.6) 

For any ^ G FqA^-^p+i^M 

Lemma C.1.5. // M has no boundary then for any regular distribution ip^ G 
FflA^M 

d^^M = (rfV)''[v^] (C.7) 
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Proof of C. 1.5. 



M 



ipAdip- / d{ip'' Alp) 

M Jm 

ipAdip- / {(p'' Alp) 

M JdM. 

ip Adip 

M 



{di^rM 



u 



C.2 Criteria for regular distributions in N-U co- 
ordinates 

Theorem C.2.1. Let the 1-form a G rA"^(A^\C) he represented in Newman-Unti 
coordinates by 

a = aidz'^, where z^ = r, z^ = R, z^ = 6, z^ = 0, (C.8) 

and where the functions a-i = ai{T,R,6,(f)) are polynomials in R and are singular 
on the worldline. Let the most divergent terms in the polynomial functions ai he 
denoted 



a,^^^. (as) 



where a[{T,6,(f)) are hounded and Pi are positive constants. The distrihution a^ G 
r^iA-'^A^, where 



c^^M = / V^Aa is finite for all ipeToA^M, (C.IO) 

'm 
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is well defined providing the four constants (3i satisfy 



l3o < 3, /3i < 2, (32 < 2, /33 < 3. 



(C.ll) 



Proof of C.2.1. An arbitrary test 3-form ^p G T^h^Ai is given in Minkowski co- 
ordinates hy (p = (pijkdy^ A dy^ A dy'^. Applying a coordinate transformation such 
that !f = ipijkdz^ A dz^ A dz^ where {z*} are NU coordinates yields the following 
form for the coefficients (pijk-. 



<^123 —R Yi23, 



<^012 — -RY012, 

and ^023 =R^^l2^ + ^'Yj^ss- 



(C.12) 



Here the functions ^\jk depend on the test functions ipijk-, sines and cosines of 

9 and 0, and the functions Ci and Cj. The key result is that they are bounded 

functions of r, 9 and cj). 

We are interested in the boundedness of a^ [if] therefore it is sufficient to show 

that / ip /\ aidz"^ is bounded for all (p G ToA^Ai. In component form we have 
Jm 



ip A (iidz^ 



M 



^-(pusOiQ + V^023ai - V^oi3a2 + Vouaajdz 



0123 



M 



where dz 



0123 _ j^O A ^~1 



(iz° A dz^ A dz"^ A dz^. 



Substituting the relations C.12 yields 



ip A cxidz^ 



M 



)2a^2 



R Y;^23'^0'^-2 



0123 



M 



+ 



+ 



RYl^^a2dz 



0123 



M 



023 + R'Yl23)&ldz'''' 



I (R'Yl 

Jm 

f RYl^^a^dz^^"^^ 
Jm 



(C.13) 
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Substituting C.9 and separating with respect to /^-dependence yields 



ip A djdz* 



M 



< 



max 



r=oo i'9=TT i'S=2-K 



a',Y\,,dz'^' 



^ dz' 



^* 



+ 



+ 



+ 



r=-oo Je=0 J(j>=0 

'T=oo pO^Tv p(^^2n 

max I / / / a[{Yl23 + Yl23)dz" 



max 



max 



r=-oo Je=0 J(j}=0 

T = 00 /•0=TT l'4>='2-K 



dz 



Rf"^ 



a'A,dz'^' 



r=-oo Je=0 J(j>=0 

T = 00 /•0=TT l'4>='2-K 



^'^^ 



^dz^ 



r=-oo J6»=0 J</)=0 



«'3Yjl2rf^°^^ 



i?'^^ 



^dz^ 



(C.14) 



We now consider the integrals w.r.t z^ = R. The standard integral 



R^ 



dR 



1 + 7-/3 



(C.15) 



J 



where e G M"*", is bounded in the limit e ^^ providing /3 < 1 + 7. Comparison 

D 



with the powers in C.14 yields the conditions C.ll 



Theorem C.2.2. Let the 2-form a G rA^(A^\C) he represented in Newman-Unti 
coordinates by 



a = ttijdz^ A dz\ where z^ = r, z^ = R, z^ = 9, z^ 



(C.16) 



and where the functions Uij = aij{T, R, 6, (p) are polynomials in R and are singular 
on the worldline. Let the most divergent terms in the functions aij be denoted 



Cxij 






(C.17) 



where a[j{T,6,(j)) are bounded, and Pij are positive constants. The distribution 
a^ e Vr,k^M, where 



a 



D\ 



if /\ a is finite for all </) G FqA^A^, 



(C.18) 



M 
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is well defined providing the six constants /3ij satisfy 



/3oi < 1, /3i2 < 2, /3i3 < 2, 
/3o2 < 2, /3o3 < 2, /323 < 3. 



(C.19) 



Proof of C.2.2. An arbitrary test 2-form G TqA^M is given by 



V? = (fijdy' A rfy^ 



(C.20) 



Applying a coordinate transformation such that ip = 'pijdz'^ A dz^ where {z*} are 
NU coordinates yields the following form for the coefficients ipij, 



(pi2 =i?Yi2, 

^03 =i?Y^3 + R^Yl,, 
'Poi =Yoi, 
and (^23 =-R^Y23. 



(C.21) 



Here as before the functions Y- ■ are bounded functions of r, ^ and 0. We are 
interested in the boundedness of tt'^[v^] therefore it is sufficient to show that 
if A Uijdz^ A dz^ is bounded for all if G TqK^M.. Hence 



M 



ip A aijdz^Adz^ 



M 



(C.22) 



(-'^13«02 + '^12tt03 + '^03"12 " '^02ttl3 + '^23"01 + 0OlOi23)dz 



0123 



M 
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Substituting the relations C.21 yields 



ip A Ciijdz^ A dz^ 



M 



HY\.,aQ2dz 



0123 



M 



EY\2aQ'idz 



0123 



M 



+ 



+ 



+ 



(i^Yjg + R^Yl^)ai2dz 



0123 



M 



{RYl^ + R^Yl2)ai^dz 



0123 



M 



R'Yl,ao,dz'''' 



M 



Y^Q^a23dz 



0123 



M 



(C.23) 



Substituting C.17 and separating with respect to /^-dependence yields 



Q.ijdz'' A dz^ A if 



M 



< 



max 



r=oo f6=TT f(h=2n 



a',2Yl,dz'^' 



+ 


max 


+ 


max 


+ 


max 


+ 


max 


+ 


max 



r=-oo Je=0 J(t)=0 

r=oo /'9=iT /'(i>=2n 



R 

^ 



c^os^Uz''' 



r=-oo Je=0 J(j>=0 
r=oo /'9=iT /'(h=2Tr 



R^°^ 



■dz^ 



^ dz' 



a[M, + Yl,)dz''' 



r=-oo Je=0 J(t)=0 
r=oo p9=it p4>=2Tr 



^^ 



^ rfz^ 



«'l3(YS2+Y^2)rf^°'' 



r=-oo J6»=0 J(j>=0 
r=oo p9=it p4>=2tt 



^''- 



^ rf^^ 



«'oiYLrf^°^^ 



r=-oo Je=0 J(t)=0 
T=oo p9=it r<t>=2it 



e p2 



r=-oo Je=0 J(t)=0 



«23YSl^^°^^ 



i?''' 

1 



R^' 



-dz^ 



(C.24) 



Once again comparing the integrals with respect to z^ = R with the standard 

D 



result C.15 yields the relations C.19 
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Dirac Geometry 



D.l Definitions 

Definition D.1.1. Consider the region A^ = N\C wliere A^ C A^ is a local 
neighborhood of the worldline. For every field point x & N there is a unique point 
td{x) at which the worldline crosses the plane of simultaneity according to an 
observer comoving with the charge at x. 

C:R^M, t^C{t) (D.l) 

Td: M^^, x^ td{x) (D.2) 

Definition D.l. 2. Dirac geometry uses a spacelike displacement vector Y = x — 
C{td{x)), which satisfies 

g{Y,C{TD)) = 0, Rl=g{Y,Y), (D.3) 

to associate a spacetime point with a point on the worldline. We observe that 
Re > is the magnitude of Y. 

Definition D.1.3. We use the notation Co = C{t£){x)). The vector fields 
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Vd, Ad, Ad G FTA^ are defined as 

VD\. = C'{rD{x)) — , ADU = C\rD{x))— and AdU = C\Mx)) — , 
oy^ oy^ oyJ 

(D.4) 
Lemma D.1.4. The exterior derivative of the Dirac time td is given by 

Proof of D.1.4. 

It follows from definition ID.1.21 that 

0=dg{Y,VD), 
=dg{x, Vd) - dg{CD, Vd), 
=Vd + {g{x, Ad) + 1- g{AD, CD))dTD, 
=VD + {g{Y,AD) + l)dTD. (D.6) 



D 



Lemma D.1.5. 



dRD = -i^ (D.7) 



R 



D 



Proof of D.1.5. 
Let 



xgTTTW, x|^ = a;"7r^ and Cd eVTM, Cd\x = C''d^—, 



— and CbGLTM, C^|, = C^ 



(D.8) 
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It follows from definition ID. 1.21 that 

dRo =dy/g{Y,Y) 
1 



dg{Y,Y) (D.9) 



2,/^,(YX) 
dg(Y, Y) =dg{x - C^, x - Co) 

=dgi^, x) + dgiCo, Cd) - 2dg{^, C^) (D.IO) 

(i5f(x,x) = d{gabx'^x^), 

= gab{dx'')x^ + gabx'^idx''), 

= Xadx'' + Xadx"', (D.H) 



Now dx"" = dy"" since x = {y ,y\y ,y ), therefore 



(i5f(x,x) = 2xady"-, 

= 2x. (D.12) 



Also 



dg{Cn,Cn) = d{gabC^oC'o), 

= idC^D)9abC'r> + {dC%)gabCl 

= ^CoaV^dTD, 

= 2g{Cn,Vn)dTD, (D.13) 



and 



dg{CD,^) = d{gahx''C\j), 

= gab{dx'')C'j, + gabX^d{C'^), 

= CDadx'' + XadiC'r)), (D-14) 
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where d{C1)) 



diC% 



„ dr = V?)dT, therefore 

OT 



dg{CD,^) = Cd + XaVgdTD, 

= Cd + g{y^,VD)dTD. 



(D.15) 



Thus 



dR 



D 



2Rd 

1 



(d^(x, x) + dg{CD, Cd) - 2dg{^, Cd)) 



R 



D 



{Y + g{Y,VD)drD) 



(D.16) 



The definition D.3 yields 



dR 



Y 



D 



R 



D 



(D.17) 



Lemma D.1.6. 



D 



dgiY^Ao) 

dgiY^Ao) 

dg{AD,AD) 



--A 



D 



--A 



D 



Vd9{Y,Ad) 
g{Y,AD) + l 

VD{g{Y,AD)+g{An.AD)) 



g{Y,An) + l 
2g{AD,AD)VD 



g{Y,An) + l 



(D.18) 



Definition D.1.7. We define the normahzed vector field 



Y 

nD = ——, where g{nD,nD) = l and ginoyVo) = 0. (D.19) 

Rd 
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D.2 The Lienard-Wiechert potential expressed 
in Dirac Geometry 

Dirac geometry is not a natural choice to use to describe electromagnetic phe- 
nomena because all retarded (and advanced) quantities are given only as Taylor 
expansions around the Dirac time td. The retarded stress form must be calculated 
as such an expansion. Below we give the advanced and retarded Lienard-Wiechert 
potentials. 

Lemma D.2.1. The difference Sr = td ~ Tr is given in terms of Rd by 

6r =Rd - Igin, C)Rl + {^g{nD, Aof + \g{nD, Ad) - ^9{Ad, Ad))RI + 0{R%) 

(D.20) 

and the difference Sa = Ta — td is given by 

6a =Rd - Igin, C)Rl + {^giriD, Ad)' - ^g{nD, Ad) - ^g{AD, Ad))RI + 0{R%) 

(D.21) 

Proof of D.2.1. 

r2 ^3 r4 

drr) = Cd- VD5r + Ad^ - Ad^ + Ad^, + 0{5l) (D.22) 

2 o 24 

and thus the null vector X is given by 

5' : 51 , 5^ 

51 : 51 ■: 5' 



X = x- C{Tr) =x-Cd + VD5r " Ad- + Ad^ - Ad^. + 0{5l), 



Y + VD5r - Ad^-^ + AD^f - Ad^^ + 0{5',). (D.23) 

2 b 24 
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Substituting (D.23) into the lightcone condition (1.61) gives 



g{X, X) =g{Y, Y) + 2g{Y, Vd)<5, - (1 + g{Y, AdM + o^(^, Ad)S: 



-{g{Y, Ao) + g{An, Ar>M + 0{6',). 



(D.24) 



Definition (D.1.2) and (D.1.7) yield 



R 



D 



g{X, X) =RJ, - (1 + RDgiriD, AnM + ^9{nD. Ad)S: 



12 



{RoginD, Ad) + g{AD, AD))5t + 0{5l). (D.25) 



We may solve this equation to obtain 5r and 5a in terms of Rr,. 
Let 5r = cliRd, then equating coefficients of order Rj^ yields 



al = 1. 



(D.26) 



We choose Sr > 0. Knowing that Rd > it follows that ai = +1. Now let 
6r = Rd + ct2-Ri) then then equating coefficients of order i?|, yields 



=2a2 + g{nD,C) 
gijiD^C) 



02 



(D.27) 



Let 5r = Rd — 2' -^g + '^s-^i); then then equating coefficients of order i?^ yields 



3 1-1 

^3 = ^9{nD, AdY + -giriD, Ad) - T^jdiAD, Ad) 
6 b 24 



(D.28) 



Thus to third order 5r is given by (D.20). 



A similar calculation may be performed in to obtain an expression for 6a = 
Ta — td in terms of Rd- In this case all quantities on the left hand side are 
evaluated at the advanced time r^, so that instead of solving the retarded null 
condition g{X,X) = we must solve the advanced null condition g{W,W) = 0. 
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Since r^ — td is positive this means that all terms with odd powers of 6 will have 
opposite sign to those in the retarded calculations. The resulting expression for 6a 



is given by (D.21). 

D 

Lemma D.2.2. In terms of the Dirac time td and the Dirac radius Rd the re- 
tarded Lienard- Wiechert potential is given by 

K = ~-^ + {Ad + -9{nD, Ad)Vd) 

+ {vB[-g{Ao,AD) - ^g{no,ADf - -(^(n^,^^)) - ^A^ - ^g{no,AD)Aj^Rr, 
+ 0{Rl), (D.29) 

and the advanced Lienard- Wiechert potential is given by 

A, =-^ + {Ad - -g{nD,AD)VD) 
no I 

+ (- ^D(g^(^D,^D) - -g{no,Aof + ^g{nD.Ao)) + -Ad - ^g{nD.AD)AT^Ro 
+ 0{R\,) (D.30) 

Proof of D.2.2. 

We evaluate the retarded Lienard-Wiechert potential as a series in R^. 

V = Vo~ Aob, + Ao^ - Aoijo)^ + 0{t^) (D.31) 



Substituting (D.20) yields 



V =Vd - AdRn + 2 (^^ + ^^^(^^' ^d))RI 



+ (^Ad {-g{nD, Ad^ + -g{nn, Ad) - ^9{Ad, Ad)) - -Ad - -g{nD, Ad)Ad^ RI 
+ OiR%) (D.32) 
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Also 



g{X,V) = - {g{Y,AD) + l)6r + g{Y,AD)f 

+ {g{An, An) - g{Y, An)) | + 0{5t) 



(D.33) 



Again substituting (D.20) yields 



g{X,V) = - Rn -lg{nD,An)Rl 



+ y-g{nD,ADf + -g{nD,AD) 



^giriD, Ad) + ^g^Ao, An)) R'n + 0{RV) 

(D.34) 



Dividing (D.32) by (D.34) gives (D.29). Evaluating the advanced potential 



A 



adv X 



C{ra) 



9{W,C{Ta)) 



(D.35) 



using the same procedure leads to (D.30) 



D 



The retarded and advanced Lienard-Wiechert fields Fj-et and Fadv are obtained 
by taking the exterior derivative of Aret and Aadv respectively. In 1938 Dirac [T7] 
showed that the difference between the retarded and advanced fields is finite on 
the worldline and given by 



^(F.et-Fadv) = ^(^((7, C')C'-C) 



(D.36) 



It is easily seen that taking the sum of expansions 



Fret — ^(Fadv + Ffct) + ^(Fadv " Ffct) 



(D.37) 



is equivalent to expanding Fj-et only. This point was emphasized by Infeld and 
Wallace [H], and later by Havas [45j . 
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Adapted N-U coordinates 



(r, r, 0, (P) 



For the numerical investigation presented in chapter [7] we use a coordinate system 
(r, r, 6*, (p) adapted from the Newman-Unti coordinates. This change in coordinates 
was initially motivated by our interest in the ultra-relativistic Lienard-Wiechert 
fields. The N-U coordinate system breaks down in the ultra-relativistic limit since 
R = —g{X, V) = when V is null. In the new coordinate system the radial 
parameter is given by 



-^ = -g(X,dyo) (E.l) 



which remains non-zero in the ultra-relativistic limit. If {y^, y^, y"^, y^) is the global 
Lorentzian coordinate chart then the coordinate transformation is given by 

yO = C\t) + r 

yi = C^(r) + rsin(e)cos(0) 

y2 = c'2(^) + rsin(^)sin((/.) 

y^ = C^{T) + rcos{9). (E.2) 

Lemma E.0.3. In terms of the new coordinates the vector fields X, V^ G rT(A^\C) 
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are given by 



X = r| (E.3) 



Proof of E. 0.3. 



X = x-C{r) 

d d d d 

= r-— 7 + rsin(6') cos(0)-— r + rsin(^) sin(0)7— r + rcos^-- 

oy^ oy^ oy^ ay 

d_ 

dr 
d dy^ d dy^ d dy^ d dy^ d 



3 



dr dr dy^ dr dy^ dr dy'^ dr dy^ 
d 



= C{t) 



Lemma E.0.4. The Minkowski metric g G (^ M zs given by 



[F,F] 
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D 



g = - c^dr (g)dr + r'^dO ® dO + r^ sin^ i 

+ a[dT ® dr + dr ® dr] + raeldr ® dO + dO ® dr] + ra(f,[dT rf0 

+ d(f)^ dr] (E.4) 



APPENDIX E. ADAPTED N-U COORDINATES 



and the inverse metric g ^ G ' M zs given by 

1/d d d d \ ae / d d d d 

a\dT dr dr drJ ar\dr 86 36 dr 
a^ / d d d d 

Q;rsin^(^) VSr del) d4> dr 



1 



( 



d d 



r^sin 6'V(90 dcj) 



(E.5) 



Where a is defined by (1.88) and ag and a^ are the derivatives of a with respect 



to 6 and respectively. Let z^ = r, z^ = r, z"^ = 6, z^ = (p, then the matrices 
G' = G'^fj = g{dza^dzb) and G'^^ = G"^^^ = g^^{dz"-, dz^) are given by 



G' = 9idz-,d^b) 



' —(? a rae ra^ * 

a 

me r^ 

\ ra^ r^ sin^ 6 J 



G 



/-I 



0-0 

a 
1 c^ sin^ {6) + aj sin^ {6) + al ae 

a 




ae 







ar 



ax axsin^{6i) 
1 

1 



arsin^(6') 







r2 sin^ 6 ) 



Proof of E. 0.4. 



g = -dy^ ® dy^ + dy^ ® dy^ + dy"^ ® dy'^ + dy'^ ® dy^ 



(E.6) 



rfyO = C^{T)dT + dr 

dy^ = C^{T)dT + sm{6) cos{(f))dr + rcos(6') cos{(f))d6 — rsin(6') sin(0)(i0 
dy'^ = C'^rdr + sm{6) sin(0)(ir + rcos(^) sin(0)(i^ + rsin(^) cos(0)(i0 
dy^ = C^rdr + cos{6)dr - rsm{6)d6 
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Thus 



g = -c^dr ®dT + r^d9 ® dd + r^ sin^ i 
+ (-6*° + C^ sin 6* cos + (7^ sin 6^ sin + C^ cos 9) [dr dr + dr (g) dr] 
+ (C^r cos 6* cos + C^r cos 6* sin - C^r sin 6) [dr ® d9 + dO ® dr] 
+ (C^r sin 6^ cos - C^r sin 6 sin 0) [dr ® dcf) + dcj) ® dr] 

= -c^dr (g)dT + r'^dO ® dO + r'^ sin^ 0)d(j) ® dcf) 

+ a[dT ® dr + dr ® dr] + rae[dT ® dO + d9 ® dr] + ra^ldr ® d(j) + d(j)® dr] 



g ^ follows from the matrix (£^.0.4) D 

Corollary E.0.5. 



dr = —dr 
a 



c^ sin^ (6) + al sin^ [9] + al 1 



sin^(6')a2 ' ' a ^ ar " arsin'^{9) 



d'^ = _,_ 2//1N -9 ^<- + T^T - -de - —TT—^TTT^^'t^ 



dO = —df) dr 



r^ ar 



r2sm''6' arsm^(6') 

Proof of E.0.5. follows from definition oi g^^. D 

Lemma E.0.6. The 1-forms X,V E FA-'^M are given by 

X = radr (E.8) 

V = —e^dr + adr + raedO + rarpd(f) (E.9) 
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Proof of E. 0.6. 



.?-4,-) 



r(-C" + C sin ^ cos + C sin ^ sin + C cos e)dT 

radr 
d . 



^'=*<Sr- 



Lemma E.0.7. 



Proof of E. 0.7. 



Lemma E.0.8. 



-c^dr + adf + raerf^ + ra 



^1 = A/|det(^)|rfr Adr Ade A 



- a"^ r^ sin^ 9\dT A dr A de A ^ 
ax^ sin 6'(ir Adx Add AdS 
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D 



^1 = -ar^ sin 6'rfr A rfr A ^61 A rf0 (E.IO) 



D 



A = adr + rdgdO + rd^dcp (E-H) 
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Proof of E. 0.8. 



^_^^«...a 



dr 



■dy" 



-C\T)dy^ + C\T)dy^ + C^{T)dy'' + C\T)dy'' 
-C^{t) \c^{T)dT + dr 



C^{T)dT + sin(^) cos(0)(ir + rcos(6') cos{4))d9 — rsin(6') sin(0)(i0 
C'^{T)dT + sin(^) sm(0)(ir + r cos(^) sin(0)(i^ + r sin(^) cos(0)(i0 
C'^(r)rfr + cos(^)rfr - r sin(^)rf^ 



+ C\t 

+ C\t 

+ C\t 

g{A, V)dT + adr + rcie(i6' + rd 

diir + rdefi^ + rd. 



Lemma E.0.9. 



a 



A = -a 



a 



(?d 



a' 



1 / . aa0\ 

+ -iag de 

r V a / 



Proof of E.0.9. 






a^ 



rsin'(^) 



a^ 



a / 



a J 



dr 



U 



(E.12) 



A = g'\A,-) 

= g{A, V)g'^{dT, -) + dg'^{dr, -) + rdeg'^iydO, -) + rd^g'^{d(j), -) 



a \ sm (6')q;^ a ar arsm (fc*) 



sin2(^)a2 

+ rdel—dg -dA + rdJ—— 

Vr ar / Vr^sm 



d^ 



sin\ef'^ arsin2(^)'^^ 



a 
a 



dr 



g{A,V) ^ cd\^ ^ fdaj ae^eA 



a 



a^ 



1 



+ 



1 / . ««e\o , 

+ -{ae )d0+ , ^ 

r V a J rsm (^) 



a^ a 

da^ 



sin (^) V a 



a J 



dr 



a 



dA 
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Lemma E.0.10. 



Proof of E. 0.10. 






dr dr 
r(-C'^ + C^ sin 9 cos 6 + C^ sin 9sm6 + C^ cos i 



ra 



For g{X, A) the only relevant term in the metric is adr ® dr, thus 



giX,A) = r-gid^,dr) 
a 

= xa 



Lemma E.0.11. 



^2 



q [otde — aa0)dT A d6 + {ad,j> — da^)dT A 
47reo a^ 



Proof of E.0.11. (53) follows directly from (20) (47) (51) 

q g{X,V)X AA-g{X,A)X AV 



F 



R 



47reo giX,Vy 
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D 



^(X,\/) = ra (E.13) 

g{X,A) = m (E.14) 



D 



A = ^(—dr+-dr+^d9+^d(P] (E.15) 

47reo Var r a a 



Fr = 5 (E.16) 



^c = —c^(%dT Adr+ ^dr Ad9+ ^dr A d(^\ (E.17) 

47ren Vr ra^ ra^ / 
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Using the relations: 

X AV = ( radr ) A f — c^dr + adr + raedO + ra;</ 

= ra^dr Adr + r'^aaedr A d9 + r^aa^rfr A dcj) (E.18) 

X A A = i radr ) A i g{A, V)dT + adr + rdgdO + rd^dcf) j 

= raddr A dr + r'^adgdr A dO + r'^ad^dr A dcp (E.19) 

along with (51)(52)gives 

q ar{raddT A dr + r'^adgdr AdO + r'^adfpdr A d(j)) 



F 



R 



47reo (ar)3 

q rdira^dr A dr + r'^aa^dr A dO + r^aa^dr A dcj)) 
47reo (a;r)3 

q 1 



47reo «^ 



ade — dag)dT AdO + {ad^ — da^)dT A 



c — 


q &X AV 
47reo^(X,V)3 
q (?{ro?dT A (ir + r'^aaedr AdO + r'^aa^dr A dcj)) 


= 


47reo (ar)3 

/ c^f \rfrArfr+ "V^Arf^+ "t^^A^) 



Lemma E.0.12. 



D 



^ Fh = ^ f T~r^' ^^ A dO - MO) (aed - ade) ^^ ^ ^ A 
47reo V a sin(^) a^ / 

.Fc = ^^^^rf^Arf0 (E.21) 
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Proof of E. 0.12. 

*(X AV) = ivix * 1 



ng^igXcur'^ sin Odr A dr A d9 A d(j)) 

r^asm{e)d9 Adr (E.22) 



-k{X AA)= lAix * 1 



2 . //ix/ . .\, ,, r^iaaA, — aA,a) ^ 

r sm.[d)[aa0 — aeajar A a0 ^-7^ — dr A dO 

sin(6') 

r^dsin(^)rf^Arf0 (E.23) 



therefore 



^Fr 



q -c^i.{XAV) 
47reo ^(X,\/)3 

g c^ sin ^ 
47reo a^ 

q g{X,V)^{X AA)-g{X,A)^{X AV) 



-dO A dcj) 



4vreo ^(X,\/)3 

g /a<^d-ad<^ sin(6')(a0d - ad^) 

-aT A ac7 r dT A 



47reo V a^ sin(6') a^ 



D 



Lemma E.0.13. The couloumhic and radiative terms of the 1-forms E and B take 
the form 

Ec =-, c ^rf0 + -iri^dr H rda* H 5-i^ar 

47reo Vro;-^ a^r^ ro;'^ a'^r"' 

"fl , '^^ , 

-I ^dr -\ ^ dr 

~ q [aa^-ad^ ade - dae 
Er =- ^ d0 de 



/aeiade - dae) _ a;^(d;Q;0 - Q;d<^)\ \ 
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and 



Be = ^cf^^^rfc/) + C'M<^)-C^cosi<^)^0 



\- 



47reo V ra 



ra" 



-csin(6') ( — -dS — 
V fa-^ 



raisin (6') 



-d9 



B 



R 



4vreL 



cAnen 



sm{9) I -^"^3 — (i(/. + "T. oTr^ cig 



ra3sin^(^) 



a3sin2(^) 



raisin (6') 



dr 



('. 



--, sin(6') d0 H , ■ 2//1N "^ 

c47reo \ a^ a^ 8111^(6') 



+ a 



raisin (6') 



(ir 



(E.25) 



Proof of E. 0.1 3. 



_d_ 



dr d dr d dO d dch d 



dy^ dr dy^ dr dy^ 86 dy^ dcj) 

c ( d -Q d ae d a^ d 

a\ OT dr r 39 rsm{6)d(p 



(E.26) 



Results follow from definitions (1.16) and (1.22). 



D 
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MAPLE Input for Part I 



In this thesis the we use the mathematical software MAPLE to implement the com- 
putations which support the results presented in parts I and II. In principle there 
are other programming tools which could have be used, such as MATHEMATICA 
and MATLAB, each of which has its own advantages and disadvantages. In gen- 
eral, MATHEMATICA and MAPLE are more suited to symbolic computation, 
whereas MATLAB is more suited to numerical computation. 

In part I of the thesis we require heavy use of symbolic computation. In 
particular we utilize the tools of differential geometry to manipulate tensors and 
differential forms. These tools were readily available to us in MANIFOLDS package 
[16] written by Robin Tucker and Charles Wang for use with MAPLE. There 
are similar packages available for use with other software, such as RICCI for use 
with MATHEMATICA, and Tensor Toolbox for use with MATLAB, however the 
availability of the MANIFOLDS package and supporting documentation was an 
important factor in deciding to use MAPLE instead of other possible programming 
tools. In addition, the procedural language of MAPLE was appealing to the author 
based on his experience with C^^ and FORTRAN programming languages. 

The calculations carried out for part II of the thesis are more numerical by 
nature, however rather than adopting a programming tool more suited to numerical 
calculations we decided it would be more economical to build on the code already 
written in MAPLE. 
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The following script was written in MAPLE 15 and can be run with the pack- 
ages Plots, Linear Algebra and the additional package Manifolds |46j with tools for 
differential geometry. 

1 # set up coordinate system 

2 Manif oldsetup(M, [taUjRjthetajphi] , [e,E,0] , 

3 map(x->simplify(x, symbolic) , 

4 [e[0]=d(tau), 

5 e[l]=d(R), 

6 e[2]=d(theta) , 

7 e[3]=d(plii)])): 



8 Constants ( [epsilon, q, ep, bO, bl, b2, b3, a3, RO] ) ; 

9 Manfdomain(M, [a, ad, ath, aph, athd, aphd, adphph, adthth] , [tau, 

theta, phi] ) : 

1 Manfdomain(M, [C0,C1 ,C2,C3,Cd0,Cdl,Cd2,Cd3,Cdd0,Cddl,Cdd2,Cdd3] , [tau]) 

2 : 

3 g := (-l+2*R*ad/a)*d(tau) &X d(tau) 

4 - (d(tau) &X d(R)+ d(R) &X d(tau)) 

5 + (R"2/a~2)*(d(theta) &X d(theta)) 

6 + (R"2/a~2)*sin(theta)"2 *d(phi) &X d(phi) : 

7 Mancovmetric(M,g) : 

8 Manvol(M) := -(R~2/a"2)*sin(theta)* '&" ' (e [0] , e[l], e [2] , e[3]) : 

9 Basisl := {d(tau) ,d(R) ,d(theta) ,d(phi)} : 

20 Basis2 := {d(tau)&~d(R) , d(tau)&"d(theta) , d(tau)&~d(phi) , 

21 d(R)&"d(theta) , d(R)&"d(phi) , d(theta)&~d(phi)} : 

22 Basis3 := {d(tau)&~d(R)&"d(theta) , d(tau)&~d(R)&"d(phi) , 

23 d(tau)&"d(theta)&~d(phi) , d(R)&"d(theta)&"d(phi)} : 

24 Basis4 := {e(0) &"e(l) &"e(2), e(l) &~e(2) &"e(3),e(2) &~e(3) 

25 &"e(0),e(3) &"e(0) &~e(l)}: 

26 a_sublist :={diff (a,tau)=ad,diff (a,theta)=ath,diff (a,phi)=aph, 

27 dif f (ath,tau)=athd,dif f (aph,tau)=aphd, 

28 dif f (ath, phi) =athph, dif f (aph, theta) =aphth, dif f (ad, theta)=athd, 

29 dif f (ad, phi)=aphd, dif f (aph, phi)=aphph, 

30 dif f (ath, theta)=athth, diff(adph, phi)=adphph, diff(adth, 

31 theta) =adthth} : 

32 Cd_sublist := {diff (CO,tau)=CdO,diff (CI ,tau)=Cdl, 

33 diff (C2, tau) =Cd2, diff (C3, tau) =Cd3} : 

34 Cd_inv_sublist := {CdO=diff (CO,tau) ,Cdl=diff (CI ,tau) , 

35 Cd2=dif f (C2 , tau) , Cd3=dif f (C3 , tau) } : 

36 Cdd_sublist := {dif f (CO,tau,tau)=CddO,diff (Cl,tau,tau)=Cddl, 

37 diff (C2 , tau , tau) =Cdd2 , diff (C3 , tau , tau) =Cdd3 , 

38 diff (CdO , tau) =CddO , diff (Cdl , tau) =Cddl , 

39 diff (Cd2, tau) =Cdd2, diff (Cd3, tau) =Cdd3} : 

40 Cddd.sublist := { diff (CddO, tau) =CdddO, diff (Cddl, tau) =Cdddl, 
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42 
43 



dif f (Cdd2 , tau) =Cddd2 , dif f (Cdd3 , tau) =Cddd3} : 

aa : = 
-Cd0+Cdl*cos(phi)*sin(theta)+Cd2*sin(phi)*sin(theta)+Cd3*cos(theta) 



45 aath := dif f (aa,theta) : 

46 aaph := diff(aa,phi) : 

47 aad := subs(Cdd_sublist ,dif f (subs(Cd_inv_sublist,aa) ,tau)) : 

48 aathd := subs (Cdd_sublist, dif f (subs (Cd_inv_sublist, aath) , tau)) : 

49 aaphd := subs (Cdd_sublist, dif f (subs (Cd_inv_sublist, aaph) , tau)) : 

50 aathth:=subs(Cdd_sublist , dif f (subs (Cd_inv_sublist, aath) ,theta)) : 

51 aaphph:=subs(Cdd_sublist , dif f (subs (Cd_inv_sublist, aaph) ,phi)) : 

52 aadphph:=subs(Cdd_sublist,diff (dif f (subs (Cd_inv_sublist, aad) ,phi) , 

53 phi)) : 

54 aadthth:=subs(Cdd_sublist,diff (dif f (subs (Cd_inv_sublist, aad) ,theta) , 

55 theta)) : 

56 aa_sublist :={a=aa, ath=aath, ad=aad, aph=aaph, athd=aathd, 

57 aphd=aaphd , athth=aathth , 

58 aphph=aaphph , adphph=aadphph, adthth=aadthth} : 



63 
64 
65 



69 
70 



72 
73 



80 
81 
82 
83 
84 



xO 
xl 
x2 
x3 



= CO-(R/a): 

= Cl-(R/a)*sin(theta)*cos(phi) 
= C2-(R/a)*sin(theta)*sin(phi) 
= C3- (R/a)*cos (theta) : 



J := Matrix (4, 4) : 

for i from to 3 do J[i+1, 1] := diff(x 



I i, tau) : 



J[i+1, 2] 
J[i+1, 3] 
J[i+1, 4] 



= diff (x II i, R) : 

= diff (x Mi, theta) : 

= diff (x II i, phi) end do: 



subs(a_sublist , J): 

DetJ := simplify(Determinant(J)) : 
detJ := (R"2/a~2)*sin(theta) : 



71 AdJ := simplify (eval (subs ( a_sublist, Adjoint (J)))) 



df_tau_0 :=AdJ[l, l]/detJ 
df_tau_l := AdJ[l, 2]/detJ 
df_tau_2 := AdJ[l, 3]/detJ 
df_tau_3 := AdJ[l, 4]/detJ 
df_R_0 := AdJ [2, l]/detJ 
df_R_l := AdJ [2, 2]/detJ 
df_R_2 := AdJ [2, 3]/detJ 
df_R_3 := AdJ [2, 4]/detJ 
df_theta_0 := AdJ [3, l]/detJ 
df_theta_l := AdJ [3, 2]/detJ 
df_theta_2 := AdJ [3, 3]/detJ 
df_theta_3 := AdJ [3, 4]/detJ 
df_phi_0 := AdJ [4, l]/detJ 
df_phi_l := AdJ [4, 2]/detJ 
df_phi_2 := AdJ [4, 3]/detJ 
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B7 df_phi_3 := AdJ[4, 4]/detJ : 

88 PD_0:=df_tau_0*PD(tau)+df_R_0*PD(R) +df _theta_0*PD(theta) 

89 +df_phi_0*PD(phi) : 

90 PD_l:=df_tau_l*PD(tau)+df_R_l*PD(R) +df _theta_l*PD(theta) 

91 +df_phi_l*PD(phi) : 

92 PD_2:=df_tau_2*PD(tau)+df_R_2*PD(R) +df _theta_2*PD(theta) 

93 +df_phi_2*PD(phi) : 

94 PD_3:=df_tau_3*PD(tau)+df_R_3*PD(R) +df_theta_3*PD(theta) 

95 +df_phi_3*PD(phi) : 

96 VX := R*PD(R) ; 

97 dualX := F2C(& (VX)) ; 

98 VV := PD(tau)+VX*(ad/a) ; 

99 dualV := collect(F2C(& (VV)), Basisl) ; 

100 dualA :=collect (expand(R*(ad"2/a"2)*d(tau) + (-ad/a)*d(R) 

101 +R* ( (ad*ath) /a"2-athd/a) *d(theta) + R* ( (ad*aph) /a"2-aphd/a) *d(phi) ) , 

102 Basisl) ; 

103 VA :=collect(expand(F2C( & (dualA))), Basis6) ; 

104 ALW := collect(expand(F2C(dualV/(-R))) , Basisl) ; 

105 FLW := collect (expand(subs(a_sublist, d(ALW))), Basis2) ; 

106 starFLW:=collect(F2C(&i (&star(FLW))) , Basis2) ; 

107 stress :=proc (kill) ; 

108 collect (subs (Cd_sublist,F2C(((ep/2)*((PD_ I I kill &i FLW) &"(&star 

109 FLW)-(PD_| Ikill &i(&star FLW))&"FLW)))), Basis3) ; 
end proc: 



stress_0 
stress_l 
stress_2 
stress_3 



=stress(0) 
=stress(l) 
=stress(2) 
=stress(3) 



expansion_cdot_sublist :={epsilon=l, Cd0=l+(b0*tau"2/2)+0(tau"3) , 
Cdl=(bl*tau"2/2)+0(tau"3) , Cd2=(b2*tau"2/2)+0(tau~3) , 
Cd3=a3*tau+(b3*tau"2/2)+0(tau"3) , Cdd0=b0*tau+0(tau"2) , 
Cddl=bl*tau+0(tau"2) , Cdd2=b2*tau+0(tau"2) , 
Cdd3=a3+b3*tau+0 (tau"2) } ; 

120 S_k_cdot :=proc(sublist, kill) 

121 local spl; 

122 spl:=stress_| Ikill; 

123 subs(sublist, subs(aa_sublist , 

124 collect (expand (subs (Cdl*cos (phi) *sin(thet a) 

125 +Cd2*sin(phi)*sin(theta)+Cd3*cos(theta)=a+CdO, 

126 -Cdl*cos(phi)*sin(theta)-Cd2*sin(phi)*sin(theta) 

127 -Cd3*cos (theta) =-a-CdO , -Cdl*cos (phi) *cos (theta) 

128 -Cd2*sin(phi)*cos(theta)+Cd3*sin(theta)=-ath,Cdl*cos(phi)*cos(theta) 

129 +Cd2*sin(phi)*cos(theta)-Cd3*sin(theta)=ath, 

130 -Cdl*sin(phi)+Cd2*cos(phi)=aph/sin(theta) ,Cdl*sin(phi) 

131 -Cd2*cos(phi)=-aph/sin(theta) , Cd_sublist ,spl)) , Basis3))); 

132 end proc: 
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137 
138 
139 
140 



intgrd_0:= series(coef f (S_k_cdot(expansion_cdot_sublist , 0), 

'&"'(d(tau), d(theta), d(phi))), tau=0) : 

intgrd_l:= series(coef f (S_k_cdot(expansion_cdot_sublist , 1), 

'&"'(d(tau), d(theta), d(phi))), tau=0) : 

intgrd_2:= series(coef f (S_k_cdot(expansion_cdot_sublist , 2), 

'&"'(d(tau), d(theta), d(phi))), tau=0) : 

intgrd_3:= series(coef f (S_k_cdot(expansion_cdot_sublist , 3), 

'&"'(d(tau), d(theta), d(phi))), tau=0) : 

141 get_integrands := procO; 

142 print (t, intgrd_0) ; 

143 print (x, intgrd_l) ; 

144 print (y, intgrd_2) ; 

145 print (z, intgrd_3) ; 

146 end proc: 

147 get_integrals := procO; print (t, 

148 f actor(simplify(int(int(int(intgrd_0, phi=0. .2*Pi) , theta=0. 

149 tau) ) ) ) ; 

150 print(x, f actor(simplify(int(int(int(intgrd_l, phi=0. .2*Pi) , 

151 theta=0. .Pi) , tau)))); 

152 printCy, f actor(simplify(int(int(int(intgrd_2, phi=0. .2*Pi) , 

153 theta=0. .Pi) , tau)))); 

154 print(z, f actor(simplify(int(int(int(intgrd_3, phi=0. .2*Pi) , 

155 theta=0. .Pi) , tau)))); 

156 end proc: 

157 get_integrals() ; 



Pi), 
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Comments 

1-7 Set up the Newman-Unti coordinate system {t,R,6,(J)) = (tau, R, theta, phi) 
8-12 The global variables are defined. For i = 0..3 we use notation C* = Ci, 
C* = Cdi, C* = Cddi. Also a = a, d = ad, oq = a.th.,a^ = aph, d^ = aphd etc. 
The constants a, 6* defining the comoving frame are given by a and bi respectively. 
Also q and ep are constants. 



13-17 The metric (1.92) is input. This associates the manifold M with Minkowski 

space A4. The function Mancovmetric(M, g) identifies g as the metric on M. The 

Manifolds package will automatically give the inverse metric and the vector and 

covector bases on TM and T*A^. Note that there is no factor of c^ in the metric 

because we use dimensions such that g{C, C) = —1. 

18 Manvol (M) defines the volume 4-form. Notice the negative orientation. 

19-25 Define coordinate bases to simplify output 

26-31 These lines define the relationships between a and its derivatives. 

32-41 These lines define the relationships between the components of C, C,C and 

C ■ 

42-58 The here we define the parameters aa, aad, aath, aaph, aaphd. . . which 

assign the coordinate representations to the variables a, ad, ath, aph, aphd... 

59-62 The coordinate transformation from Newman-Unti (tau, R, theta, phi) to 

Minkowski coordinates (xO, xl, x2, x3). 

63-70 We determine the Jacobian J and its determinant. 

71-87 We calculate the partial derivatives of the Newman-Unti coordinates with 

respect to the Minkowski coordinates. 

88-95 These lines define the Minkowski basis vectors PD_t=Tj|Tj,PD_x=g|Y) ^^-J—Im^ 

PD_z=^ in terms of Newman-Unti coordinates. 

96-103 Defines the vectors X = VX, V = VV, and A = Yk and their duals using 



(1.90) and (1.91) and (1.99). 



104-106 The Lienard-Wiechert potential A = ALW is defined using (1.106). The 



2-form field F = FLW may is calculated by taking the exterior derivative. This 
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is included in the Manifolds package. The Hodge dual is also used to define 
^F = starFLW 

107-114 These lines define the four stress 3-forms S;^ =stress_i for i = 0, 1, 2, 3. 
115-119 Defines the expansion around the momentarily comoving frame 
120-32 A procedure for substituting the expansion into either of the stress 3-forms 
and simplifying the resulting expression. 

133-145 These lines provide the procedure get .integrands for obtaining the inte- 
grands. 

147-156 These lines provide the procedure get_integrals for carry out the inte- 
gration. 



157 This calls the procedure get_integrals. The result is stated in (4.5) 
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MAPLE Input for Part II 



For the numerical investigation in Part II we use MAPLE to perform many different 
calculations, integrals and plots for a wide range of input parameters. As a result 
I have many different files with variations on the code. With hindsight I would 
have liked to have kept all the code in one file, beautifully annotated and ready 
to reproduce any calculation. However coding in MAPLE is a skill which I have 
learnt throughout my PhD and the code I have written hasn't always been the most 
simple or the most elegant. In this section I present some of the most important 
code which has been used to obtain the results stated in chapter [7| Once again 
we use the packages Plots, Linear Algebra and Manifolds\ 

G.l Part 1 - Setup 



1 # set up coordinate system 

2 Manif oldsetup(M, [tau,r,theta,phi] , [e,E,0] , 

3 map(x->simplify(x, symbolic) , 

4 [e[0]=d(tau), 

5 e[l]=d(r), 

6 e[2]=d(theta) , 

7 e[3]=d(phi)])): 



Constants (epsilon, Lp, Rp, thetap, v, XO, YO, ZO, q_e, ep, mu, c) 
Manf domain (M,gAV) : 

Manfdomain(M, [a,ath,aph] , [tau,theta,phi] ) : 
11 Manf domain (M, [ad, athd,aphd, athph, aphth] , [tau,theta,phi] ) : 
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12 ManfdomainCM, [C0,C1 ,C2,C3,CdO,Cdl ,Cd2,Cd3,CddO,Cddl,Cdd2,Cdd3] , [tau] ) 

13 : 

14 ManfdomainCM, [rhat, cthhat, sthhat, cphhat, sphhat, TO], [tau]): 

15 g := -c"2*d(tau) & X d(tau) 

16 + a*(d(tau) & X d(r)+ d(r) & X d(tau)) 

17 + r*ath *(d(tau) & X d(theta)+ d(theta) & X d(tau)) 

18 + r*aph*(d(tau) & X d(phi)+ d(phi) & X d(tau)) 

19 + r"2*(d(theta) & X d(theta)) 

20 + r"2*sin(theta)"2 *d(phi) & X d(phi) : 

21 Mancovmetric(M,g) : 

22 G:=Manconmetric(M) : 

23 Cd.sublist := {diff (CO,tau)=CdO,diff (CI ,tau)=Cdl, 

24 diff (C2, tau) =Cd2, diff (C3, tau) =Cd3} : 

25 Cd_inv_sublist := {CdO=diff (CO,tau) ,Cdl=dif f (CI ,tau) , 

26 Cd2=dif f (C2 , tau) , Cd3=dif f (C3 , tau) } : 

27 Cdd.sublist := {dif f (CO,tau,tau)=CddO,dif f (Cl,tau,tau)=Cddl, 

28 diff (C2, tau, tau) =Cdd2, diff (C3, tau, tau) =Cdd3} : 

29 aa : = 

30 -Cd0+Cdl*cos(phi)*sin(theta)+Cd2*sin(phi)*sin(theta)+Cd3*cos(theta) 

31 : 

32 aath := dif f (aa,theta) : 

33 aaph := diff (aa, phi) : 

34 aad := subs(Cdd_sublist , dif f (subs (Cd_inv_sublist, aa) , tau)) : 

35 aathd := subs (Cdd_sublist, diff (subs (Cd_inv_sublist, aath) , tau)) : 

36 aaphd := subs(Cdd_sublist,diff (subs (Cd_inv_sublist, aaph) , tau)) : 

37 Basisl := {d(tau) ,d(r) ,d(theta) ,d(phi)} : 

38 Basis2 := {d(tau)& "d(r), d(tau)& ~d(theta) , d(tau)& ~d(phi) , d(r)& 

39 "d(theta), d(r)& ~d(phi) , d(theta)& "d(phi)} : 

40 Basis3 := {d(tau)& "d(r)& ~d(theta) , d(tau)& "d(r)& "d(phi), 

41 d(tau)& "d(theta)& *d(phi) , d(r)& "d(theta)& "d(phi)} : 

42 VX := r*PD(r) : 

43 VV := PD(tau) : 

44 dualX := & (VX) : 

45 dualV : = & (VV) : 

46 dualA := subs(gAV *d(tau) + ad*d(r) + r*athd*d(theta) + 

47 r*aphd*d(phi) ) : 

48 VA := & (dualA) : 

49 ALW : = q_e* (dualV/ (r*a) ) : 

50 FLW : = 

51 collect (subs({dif f (a,tau)=ad,diff (a,theta)=ath,dif f (a,phi)=aph, 

52 dif f (ath,tau)=athd,dif f (aph,tau)=aphd, 

53 diff (ath , phi ) =athph , dif f (aph , theta) =athph} , 

54 simplify (F2C(d(ALW)))),Basis2) : 

55 FLWc := subs (ad=0,athd=0,aphd=0, FLW) : 

56 FLWr := collect (simplify(FLW - FLWc) ,Basis2) : 
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67 
68 



xO 
xl 
x2 
x3 



= (CO+r)/c: 

= Cl+r*sin(theta)*cos(phi) 
= C2+r*sin(theta)*sin(phi) 
= C3+r*cos(theta) : 



J := Matrix (4, 4) : 

J[l, 1] :=CdO/c: 

for i from 1 to 3 do J[i+1, 1] 

J[i+1, 2] 

J[i+1, 3] 

J[i+1, 4] 



= Cd II i : 

= diffCx II i, r): 

= diff (x Mi, theta) : 

= diff (x II i, phi) end do: 
J[l,2] :=diff (xO, r) : J[l,3] :=diff (xO, theta) : J[l ,4] :=dif f (xO, phi) 
J: 



DetJ := simplify(Determinant(J)) : 
detJ := -(l/c)*a*r"2*sin(theta) : 
ManvoKM) :=-(l/c) a*r"2*sin(theta)*'&" ' (e [0] , e[l], e [2] , e[3]) 

AdJ := simplify(Adjoint(J)) : 



73 


df_tau_1 


: : = 


AdJ[l, l]/detJ 




74 


df _tau_} 


c := AdJ[l, 2]/detJ 




75 


df _tau_] 


[ := AdJ[l, 3]/detJ 




76 


df _tau_5 


2 := AdJ[l, 4]/detJ 




77 


#df_r_t 


:= AdJ [2, l]/detJ : 


78 


df_r_t 


= ((CdO+a)*c)/a : 


79 


df_r_x 


= AdJ [2, 2]/detJ : 


80 


df-r_y 


= AdJ [2, 3]/detJ : 


81 


df_r_z 


= AdJ [2, 4]/detJ : 


82 


#df_the 


ta_t := AdJ [3, l]/detJ 


83 


df.thet 


a_t 


= (ath*c)/(r*a) : 


84 


df_thet 


a_x 


= AdJ [3, 2]/detJ : 


85 


df_thet 


a-y 


= AdJ [3, 3]/detJ : 


86 


df.thet 


a_z 


= AdJ [3, 4]/detJ : 


87 


df_phi_1 


: := AdJ[4, l]/detJ 




88 


df_phi_} 


c := AdJ [4, 2]/detJ 




89 


df_phi_] 


r := AdJ [4, 3]/detJ 




90 


df_phi_2 


^ : = 


AdJ [4, 4]/detJ 





91 PD_t:=df_tau_t*PD(tau)+df_r_t*PD(r) +df_theta_t*PD (theta) 

92 +df_phi_t*PD(phi) : 

93 PD_x:=df_tau_x*PD(tau)+df_r_x*PD(r) +df_theta_x*PD (theta) 

94 +df _phi_x*PD(phi) : 

95 PD_y:=df_tau_y*PD(tau)+df_r_y*PD(r) +df_theta_y*PD (theta) 

96 +df _phi_y*PD(phi) : 

97 PD_z:=df_tau_z*PD(tau)+df_r_z*PD(r) +df_theta_z*PD (theta) 

98 +df_phi_z*PD(phi) : 

99 PDt_Fc:=PD_t &i FLWc : 

100 PDt_starFc:=collect(subs(aph=aaph,Cd_sublist ,F2C(PD_t &i 

101 (&star(FLWc)))) , Basisl,simplify) : 

102 Elec_c : = (l/c)*PD_t &i FLWc : 
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103 Mag_c : = (l/(c*c))*collect(subs(aph=aaph,Cd_sublist,F2C(PD_t &i 

104 (&star(FLWc)))) , Basisl, simplify) : 

105 Elec_r := (l/c)*collect(PD_t &i FLWr, Basisl) : 

106 Mag_r := (l/(c*c))*collect(PD_t &i F2C(&star(FLWr)) , Basisl) : 

107 Elec_cx := simplify (PD_x &i Elec_c) 

108 Elec_cy := simplify (PD_y &i Elec_c) 

109 Elec_cz := simplify (PD_z &i Elec_c) 
no Elec_rx := simplify (PD_x &i Elec_r) 

111 Elec_ry := simplify (PD_y &i Elec_r) 

112 Elec_rz := simplify (PD_z &i Elec_r) 

113 Mag_cx := simplify (PD_x &i Mag_c) 

114 Mag_cy := simplify (PD_y &i Mag_c) 

115 Mag_cz := simplify (PD_z &i Mag_c) 

116 Mag_rx := simplify (PD_x &i Mag_r) 

117 Mag_ry := simplify (PD_y &i Mag_r) 

118 Mag_rz := simplify (PD_z &i Mag_r) 

119 Energy _res : = (l/2)*(ep*((Elec_cx+Elec_rx) "2+(Elec_cy+Elec_ry) "2 

120 + (Elec_cz+Elec_rz) "2) + (1/mu) * ( (Mag_cx+Mag_rx) "2+ (Mag_cy+Mag_ry) "2 

121 + (Mag_cz+Mag_rz) ~2) ) : 

122 hat_ sublist := {TO=(sqrt((XO-Cl) "2 + (Y0-C2)"2 + (Z0-C3)"2) +CO)/c, 

123 rhat=sqrt((X0-Cl)"2 + (Y0-C2)"2 + (Z0-C3)"2), 

124 cthhat=(Z0-C3)/(sqrt((X0-Cl)~2 + (Y0-C2)"2 + (Z0-C3)~2)), 

125 sthhat=sqrt((X0-Cl)"2+(Y0-C2)"2)/(sqrt((X0-Cl)"2 + (Y0-C2)"2 + 

126 (Z0-C3)"2)), 

127 cphhat=(X0-Cl)/(sqrt((X0-Cl)~2+(Y0-C2)"2)), 

128 sphhat= (Y0-C2) / (sqrt ( (XO-Cl) "2+ (Y0-C2) "2) ) } : 

129 prehat_ subslist :=cos(theta)=cthhat,sin(theta)=sthhat , 

130 cos(phi)=cphhat ,sin(phi)=sphhat,r=rhat : 

131 Curve3_ def := { 

132 COa=epsilon*gamma*tau, 

133 Cla=epsilon*gamma*v*tau, 

134 C2a=0 , 

135 C3a=0 } : 

136 Curve2_ def := { 

137 COa=epsilon*gamma*tau, 

138 Cla=Lp-epsilon*(Rp*sin((Lp/(epsilon*Rp))-(gamma*v*tau)/Rp)) , 

139 C2a=epsilon*Rp* ( 1-cos ( (Lp/ (epsilon*Rp) ) - (gamma* v*tau) /Rp) ) , 

140 C3a=0} : 

141 Curve 1_ def :={ 

142 COa=epsilon*gamma*tau, 

143 Cla=epsilon*(gamma*v*cos(thetap)*tau+Lp 

144 -Rp*sin(thetap)+cos(thetap)*(thetap*Rp-Lp/epsilon)) , 

145 C2a=epsilon*(-gamma*v*sin(thetap)*tau + 

146 Rp* (1-cos (thetap))-sin(thetap)*(thetap*Rp-Lp/epsilon)) , 

147 C3a=0} : 

148 Curve3_ sublist := eval (subs (Diff=diff,eval (subs (Curve3_ def, 
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149 I 

150 CO=COa,CdO=Diff (COa,tau) ,CddO=Dif f (COa,tau,tau) , 

151 Cl=Cla,Cdl=Diff (Cla,tau) ,Cddl=Dif f (Cla,tau,tau) , 

152 C2=C2a,Cd2=Diff (C2a,tau) ,Cdd2=Dif f (C2a,tau,tau) , 

153 C3=C3a,Cd3=Diff (C3a,tau) ,Cdd3=Dif f (C3a,tau,tau)} 

154 )))) : 

155 Curve2_ sublist := eval (subs (Diff=diff ,eval (subs (Curve2_ def, 

156 I 

157 CO=COa,CdO=Diff (COa,tau) ,CddO=Dif f (COa,tau,tau) , 

158 Cl=Cla,Cdl=Diff (Cla,tau) ,Cddl=Dif f (Cla,tau,tau) , 

159 C2=C2a,Cd2=Diff (C2a,tau) ,Cdd2=Dif f (C2a,tau,tau) , 

160 C3=C3a,Cd3=Diff (C3a,tau) ,Cdd3=Dif f (C3a,tau,tau)} 

161 )))) : 

162 Curvel_ sublist := eval (subs (Diff=diff , eval (subs (Curvel_ def, 

163 I 

164 CO=COa,CdO=Diff (COa,tau) ,CddO=Dif f (COa,tau,tau) , 

165 Cl=Cla,Cdl=Diff (Cla,tau) ,Cddl=Dif f (Cla,tau,tau) , 

166 C2=C2a,Cd2=Diff (C2a,tau) ,Cdd2=Dif f (C2a,tau,tau) , 

167 C3=C3a,Cd3=Diff (C3a,tau) ,Cdd3=Dif f (C3a,tau,tau)} 

168 )))) : 



169 get_range3 := proc(Values_sublist) 

170 local Taub ; 

171 Taub := subs(Values_sublist,XO/(epsilon>Kgainma*v)) ; 

172 0. .Taub ; 

173 end proc : 

174 get_range2 := proc(Values_sublist) 

175 local Taua ; 

176 Taua := subs(Values_sublist,-Rp*thetap/(gainma*v)) ; 

177 Taua . . ; 

178 end proc : 

179 get_rangel := proc(Values_sublist) 
ISO local Taua ; 

181 Taua := subs(Values_sublist,-Rp*thetap/(gainma*v)) ; 

182 subs (Values.sublist, StartTau) . .Taua ; 

183 end proc : 



184 Get_Fields := proc(Cnum,Values_sublist) 

185 local Curve_sublist ; 

186 Curve_sublist := Curve I I Cnum I I .sublist ; 

187 I 

188 Elec_cx_res = 

189 subs (Values_sublist , subs (Curve_sublist , 

190 subs (hat_sublist , subs (prehat_subslist , 

191 subs(a=aa, ad=aad, ath=aath, athd=aathd, aph=aaph, aphd=aaphd, 

192 Elec_cx))))) , Elec_cy_res = 

193 subs (Values_sublist , subs (Curve_sublist , 

194 subs (hat_sublist , subs (prehat_subslist , 
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195 subs(a=aa, ad=aad, ath=aath, athd=aathdj 

196 Elec_cy))))) , 

197 Elec_cz_res = 

198 subs (Values_sublist , subs (Curve_sublist , 

199 subs (hat.sublist , subs (prehat_subslist , 

200 subs(a=aa, ad=aad, ath=aath, athd=aathdj 

201 Elec_cz))))) , 

202 Elec_rx_res = 

203 subs (Values_sublist , subs (Curve_sublist , 

204 subs (hat.sublist , subs (prehat_subslist , 

205 subs(a=aa, ad=aad, ath=aath, athd=aathdj 

206 Elec_rx) ) ) ) ) , 

207 Elec_ry_res = 

208 subs (Values_sublist , subs (Curve_sublist , 

209 subs (hat.sublist , subs (prehat_subslist , 

210 subs(a=aa, ad=aad, ath=aath, athd=aathdj 

211 Elec_ry))))) , 

212 Elec_rz_res = 

213 subs (Values_sublist , subs (Curve_sublist , 

214 subs (hat.sublist , subs (prehat_subslist , 

215 subs(a=aa, ad=aad, ath=aath, athd=aathdj 

216 Elec_rz))))) , 

217 Mag_cx_res = 

218 subs (Values_sublist , subs (Curve_sublist , 

219 subs (hat_sublist , subs (prehat_subslist , 

220 subs(a=aa, ad=aad, ath=aath, athd=aathdj 

221 Mag_cx) ) ) ) ) , 

222 Mag_cy_res = 

223 subs (Values_sublist , subs (Curve_sublist , 

224 subs (hat_sublist , subs (prehat_subslist , 

225 subs(a=aa, ad=aad, ath=aath, athd=aathd, 

226 Mag_cy) ) ) ) ) , 

227 Mag_cz_res = 

228 subs (Values_sublist , subs (Curve_sublist , 

229 subs (hat_sublist , subs (prehat_subslist , 

230 subs(a=aa, ad=aad, ath=aath, athd=aathd, 

231 Mag_cz))))) , 

232 Mag_rx_res = 

233 subs (Values_sublist , subs (Curve_sublist , 

234 subs (hat_sublist , subs (prehat_subslist , 

235 subs(a=aa, ad=aad, ath=aath, athd=aathd, 

236 Mag_rx) ) ) ) ) , 

237 Mag_ry_res = 

238 subs (Values_sublist , subs (Curve_sublist , 

239 subs (hat_sublist , subs (prehat_subslist , 

240 subs(a=aa, ad=aad, ath=aath, athd=aathd, 

241 Mag_ry) ) ) ) ) , 

242 Mag_rz_res = 



aph= aaph , aphd= aaphd , 



aph=aaph, aphd=aaphd, 



aph=aaph, aphd=aaphd, 



aph=aaph, aphd=aaphd, 



aph=aaph, aphd=aaphd, 



aph=aaph, aphd=aaphd, 



aph= aaph , aphd= aaphd , 



aph= aaph , aphd= aaphd , 



aph= aaph , aphd= aaphd , 



aph= aaph , aphd= aaphd , 
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243 subs (Values_sublist , subs (Curve_sublist , 

244 subs (hat_sublist , subs (prehat_subslist , 

245 subs(a=aa, ad=aad, ath=aath, athd=aathd, aph=aaph, aphd=aaphd, 

246 Mag_rz))))) , 

247 TO_res = 

248 subs (Values_sublist , subs (Curve_sublist , 

249 subs (hat_sublist , subs (prehat_subslist , 

250 subs(a=aa, ad=aad, ath=aath, athd=aathd, aph=aaph, aphd=aaphd, 

251 TO))))) , 

252 Cl_res = 

253 subs (Values_sublist , subs (Curve_sublist , 

254 subs (hat_sublist , subs (prehat_subslist , 

255 subs(a=aa, ad=aad, ath=aath, athd=aathd, aph=aaph, aphd=aaphd, 

256 CI))))) , 

257 aa_res = 

258 subs (Values_sublist , subs (Curve_sublist , 

259 subs (hat_sublist , subs (prehat_subslist , 

260 subs(a=aa, ad=aad, ath=aath, athd=aathd, aph=aaph, aphd=aaphd, 

261 a))))) , 

262 Energy _res = 

263 subs (Values_sublist , subs (Curve_sublist , 

264 subs (hat_sublist , subs (prehat_subslist , 

265 subs(a=aa, ad=aad, ath=aath, athd=aathd, aph=aaph, aphd=aaphd, 

266 sqrt ( (Elec_cx+Elec_rx) ~2+ (Elec_cy+Elec_ry) "2+ (Elec_cz+Elec_rz) "2) ) ) ) ) ) 

267 } ; 

268 end proc: 



Comments 

1-98 This section of the code is almost identical to that of Part I with a few notable 
exceptions: 



We use the coordinate system given by (E.2) where (r, r,6',(/)) = (tau, r, 
theta, phi). 



The metric is now given by (E.4) 



• We We define Fc =FLWc by setting all components of acceleration to zero in 
FLW. We define Fr =FLWr as the difference FLW-FLWc. 
99-106 Calculate Ec =Elec_c, Er =Elec_r, Be =Mag_c and Br =Mag_r. This is 



easily done using (6.12). 



107-118 Calculate the components of these vectors in the xl , x2 and x3 directions 
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by taking the internal contractions with respect to PD_x, PD_y and PD_z respectively. 
119-121 Calculate the total energy of the electric field | |E(r, r, ^, 0) | p =Energy_res. 



122-130 Input the substitutions given by (7.4) 



131-147 Define the three sections of the pre-bent path by inputting the components 



(7.1) according to (7.5). The labels for the axis in the code are different to those 



given in figure 5.5| due to the way I initially set up the trajectory. The axes x, y, z 
in the figure correspond to y, z, x in the code, and correspondingly the point 
X = {Xq,Yq, Zq) is given by (YO, ZO, XO). The coordinate system is aligned so 



that instead of being located at the terminus of the bend as in the figure (5.5) 



the origin is located at the end of the small straight line section. As a result the 

parameter Lp is defined as the negative of the distance Z. We use notation R =Rp 

and 6=thetap. 

148-168 Define the three corresponding list of substitutions which will associate a 

field with a particular trajectory. 

169-183 Calculate the ranges of r =tau for each of the three sections of the path. 

The values Taua and Taub are the tau values at the start and the end of the bend 

respectively. The value StartTau is the initial value for tau. 

184-186 The procedure get .fields uses the substitutions in the previous section 

to output the listed fields as functions of r =tau for a given section of path and a 

given set of input parameters. The inputs are the number Cnum=l, 2 or 3 which 

tells maple which of the three sections of the path 131-147 we are considering, and 

a list of numerical inputs of the following format 

Values_sublistO := 

subs (gain=1000 , {X0=0 , Y0=0 , Z0=1 , epsilon=l , v=sqrt ( l-l/gam"2) , Lp=0 , 

Rp=1000, thetap=0.1, gamma=gajii, StartTau=-20}) ; 

247-251 Notice the lab time To(r) =TO_res and the total energy of the electric 
field \\S{t, r, 0,0) IP =Energy_res are also obtained as functions of tau. 



G.2 Part 3 - minimize peak field 



269 get_list3 := proc(Values_sublist) 

270 local Taub ; 

271 Taub := subs(Values_sublist,Lp/(epsilon>Kgainma*v)) ; 

272 $ (round (Taub) . .0) ; 

273 end proc : 

274 get_list2 := proc(Values_sublist) 

275 local Taua, Taub ; 

276 Taub := subs(Values_sublist,Lp/(epsilon*gainma*v)) ; 

277 Taua : = 

278 subs(Values_sublist,Lp/(epsilon*gamma*v)-Rp*thetap/(gainma*v)) ; 

279 $ (round (Taua) . .round (Taub) ); 

280 end proc : 

281 get.listl := proc(Values_sublist) 
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2B2 local Taua, Taub ; 

2B3 Taub := subs(Values_sublist ,Lp/(epsilon*gainma*v)) ; 

2B4 Taua : = 

285 subs (Values_sublist , Lp/ (epsilon*gamma*v) -Rp*thetap/ (gainma*v) ) ; 

286 $(round(subs(Values_sublist,StartTau)) . .round(Taua) ) 

287 end proc : 

288 Max_ field := proc(Values_ sublist) 

289 local Fieldl,Field2, Fields, taurngl,taurng2,taurng3, FUNCTl , FUNCT2, 

290 FUNCT3,Taua, Taub,VALl, VAL2, VAL3, il, 12, iS, Fl, FFl, F2, FF2; 

291 Taub := evalf (subs (Values, sublist ,Lp/(epsilon*gaiimia*v))) ; 

292 Taua := evalf (subs (Values. 

293 sublist ,Lp/(epsilon*ganma*v)-Rp*thetap/(gamma*v))) ; 

294 



301 
302 
303 
304 
305 



Fieldl 

295 Field2 

296 Fields 

297 taurngl 

298 taurng2 

299 taurngS 



=evalf (subs(Get_ Fields (l,Values_ sublist) ,Energy_ res)) 
=evalf (subs(Get_ Fields (2, Values, sublist) , Energy, res)) 
=evalf (subs(Get_ Fields (S, Values, sublist) , Energy, res)) 

= get. listKValues. sublist) 

= get_ list2(Values_ sublist) 

= get_ listS(Values_ sublist) 
VALl :=(abs (subs (Values, sublist, StartTau))-(abs (round (Taua)))) : 

for il from 1 to VALl do: 

for 12 from 1 to 100 do: 

FUNCT1:= max(subs(tau=taurngl[il] , Fieldl), 

subs (tau=taurngl[il] -12/100, Fieldl)) ; 

end do : 

306 end do : 

307 ARR:=Array(l. .19) : 

308 #for il from 1 to (abs (round (Taua))) do: 

309 for IS from 1 to 19 do: 

310 ARR [IS] : = (evalf (subs (tau=-iS* (0.05) , Field2))); 

311 FUNCT2:=max(ARR); 

312 end do : 

313 #end do : 

314 max (FUNCTl, FUNCT2) ; 

315 end proc : 

316 Values.sublistl : = 

317 subs (gain=1000 , X0=0 , Y0=0 , Z0=1 , epsilon=l , v=sqrt (l-l/gam~2) , 

318 gamma=gam , Lp=0 , Rp=Rpp , thetap=thetapp, StartTau=-20) ; 

319 thetap.range:= 1/95, 1/90, 1/85, 1/80, 1/75, 1/70, 1/65, 1/60, 1/55, 

320 1/50, 1/45, 1/40, 1/S5, 1/SO, 1/25, 1/20, 1/15, 1/10, 1/5, 1; 

321 Rp.range : =500 , 1000, 1500, 2000, 2500, SOOO, 3500, 4000, 4500, 5000, 

322 5500, 6000, 6500, 7000, 7500, 8000, 8500, 9000, 9500, 10000; 

323 B:=Matrix(20, 20); 

324 for 11 from 1 to 20 do: 

325 for 12 from 1 to 20 do: 

326 subs.LthetaR : =subs(thetapp=thetap .range [il] ,Rpp=Rp .range [12] , 

327 Values_sublistl) ; 

328 B[il, 12] := Max_field(subs_LthetaR) ; 

329 end do ; 
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330 end do ; 



Comraents 



269-287 The procedures get_list I iCnum will round the values of Taua and Taub 
to the nearest integer and output the range of tau as a sequence of integers. 
288-315 The procedure Max_ field will compare the peak value of Energy, res 
for the initial straight line and the bend for a number of values of tau. The field 
given by the second straight line is negligible. The peak field for the straight 
segment is given by the local variable FUNCTl and the peak field for the bend is 
given by FUNCT2. 

316-330 These lines of code will create a 20 x 20 matrix B whose elements are 
the peak fields corresponding to the given values of thetap and Rp. These values 



correspond to those given in table 7.1 and the resulting matrix was used to plot 



figure 7.1 using the MAPLE function matrixplot. 



G.3 Part 4 - Plots 



331 Plot_Field_tau := proc(Values_sublist, Field) 

332 local Fieldl,Field2, Fields, taurngl,taurng2,taurng3; 

333 Fieldl :=subs(Get_Fields(l,Values_sublist) , Field) 

334 Field2:=subs(Get_Fields(2,Values_sublist) , Field) 

335 Fields :=subs(Get_Fields (3, Values_sublist) , Field) 

336 taurngl := get_rangel(Values_sublist) 

337 taurng2 := get_range2(Values_sublist) 

338 taurngS := get_rangeS(Values_sublist) 

339 display ( 

340 plot (Fieldl ,tau=taurngl,color=BLACK,_rest) , 

341 plot (Field2,tau=taurng2,color=RED,_rest , numpoints=1000) , 

342 plot (Fields , tau=taurngS , color=BLUE , _rest ) 

343 ) : 

344 end proc: 

345 Plot_Field_TO := proc (Values_sublist, Field) 

346 local Fieldl, Field2, Fields, taurngl, taurng2,taurng3; 

347 taurngl := get_rangel(Values_sublist) ; 

348 taurng2 := get_range2(Values_sublist) ; 
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349 taurngS := get_range3(Values_sublist) ; 

350 Fieldl :=subs(Get_Fields(l,Values_sublist) , [TO_res, Field, tau=taurngl] ) 

351 ; 

352 Field2:=subs(Get_Fields(2,Values_sublist) , [TO_res, Field, tau=taurng2] ) 

353 ; 

354 Fields :=subs(Get_Fields (3, Values_sublist) , [TO_res, Field, tau=taurng3] ) 

355 ; 

356 display ( 

357 plot (Fieldl, color=BLACK,_rest) , 

358 plot (Field2 , color=RED , _rest) , 

359 plot (Fields , color=BLUE , _rest ) 



360 
361 



386 



) 



end proc : 

362 Values_sublistl : = 

363 subs (gam= 1000, 

364 {X0=0 . 005 , Y0=0 , Z0=0 . 0005 , epsilon=l , v=sqrt ( l-l/gain"2) , gaiimia=gam , Lp=0 , thetap=0 . 13 

365 StartTau=-100 

366 , c=3*10"(8), q_e=-l. 80951262*10" (-8)}); 

367 Values_sublist2 : = 

368 subs (gain= 1000, {X0=0.005,Y0=0, Z0=0.0005, 

369 epsilon=l ,v=(sqrt(l-l/gajn"2)) ,gamma=gain,Lp=0,thetap=0,Rp=0.5, 

370 StartTau=-100, c=3*10"(8), 

371 q_e=(-1.80951262*10"(-8))}); 

372 EEx:=subs(Get_Fields(l,Values_sublist2) , Elec_cx_res+Elec_rx_res) 

373 EEy :=subs(Get_Fields(l,Values_sublist2) , Elec_cy_res+Elec_ry_res) 

374 EEz:=subs(Get_Fields(l ,Values_sublist2) , Elec_cz_res+Elec_rz_res) 

375 TT:=subs(Get_Fields(l,Values_sublist2) , (TO_res)) : 

376 EExl :=subs (Get .Fields (l,Values_subli St 1) ,Elec_cx_res+Elec_rx_res) 

377 EEyl :=subs (Get .Fields (l,Values_subli St 1) ,Elec_cy_res+Elec_ry_res) 

378 EEzl :=subs(Get_Fields(l,Values_sublistl) ,Elec_cz_res+Elec_rz_res) 

379 TTl:=subs(Get_Fields(l,Values_sublistl) , (TO_res)) : 



380 EEx2 

381 EEy2 

382 EEz2 



=subs (Get_Fields (2 , Values_sublist 1) , Elec_cx_res+Elec_rx_res) 
=subs(Get_Fields(2,Values_sublistl) ,Elec_cy_res+Elec_ry_res) 
=subs (Get .Fields (2 , Values_sublist 1) , Elec_cz_res+Elec_rz_res) 



383 TT2:=subs(Get_Fields(2,Values_sublistl) , (TO_res)) : 



384 EEx3 

385 EEy3 



EEz3 



=subs (Get .Fields (3 , Values_sublist 1) , Elec_cx_res+Elec_rx_res) 
=subs(Get_Fields(3,Values_sublistl) ,Elec_cy_res+Elec_ry_res) 
=subs(Get_Fields(3,Values_sublistl) ,Elec_cz_res+Elec_rz_res) 



TT3:=subs(Get_Fields(3,Values_sublistl) , (TO_res)) : 



388 partlx:=plot([10"(12)*TTl, abs(EExl), 

389 tau=get_raiigel(Values_sublistl)] , color=black, numpoints=10000) : 

390 part2x:=plot([10"(12)*TT2, abs(EEx2), 

391 tau=get_range2(Values_sublistl)] , color=red,resolution=600, 

392 numpoints=50000) : 

393 part3x:=plot([10"(12)*TT3, abs(EEx3), 

394 tau=get_rajige3(Values_sublistl)] , color=blue, numpoints=10000) : 
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400 
401 
402 
403 



395 partly :=plot( [10" (12) *TT1, abs(EEyl), 

396 tau=get_raiigel(Values_sublistl)] , color=black, numpoints=10000) 

397 part2y:=plot([10~(12)*TT2, abs(EEy2), 

398 tau=get_range2(Values_sublistl)] , color=red,resolution=600, 

399 numpoints=50000) : 
part3y:=plot([10"(12)*TT3, abs(EEy3) , 

tau=get_raiige3(Values_sublistl)] , color=blue, numpoints=10000) : 
partlz:=plot([10"(12)*TTl, abs(EEzl) , 
tau=get_rangel(Values_sublistl)] , color=black, numpoints=10000) 

404 part2z:=plot([10"(12)*TT2, abs(EEz2), 

405 tau=get_range2(Values_sublistl)] , color=red,resolution=600, 

406 numpoints=50000) : 

407 part3z:=plot([10"(12)*TT3, abs(EEz3), 

408 tau=get_range3(Values_sublistl)] , color=blue, numpoints=10000) : 
Resize(display(partlx, part2x, part3x, axes=boxed, 
view=[15.8. .16.8, 0..8], axesfont= [TIMES, ROMAN, 20], 
thickness=3)) ; 

Resize (display (partly, part2y, part3y, axes=boxed, 
view=[15.8. .16.8, 0..8], axesfont= [TIMES, ROMAN, 20], 

414 thickness=3)) ; 

415 Resize(display(partlz, part2z, part3z, axes=boxed, 

416 view=[15.8. .16.8, 0..8], axesfont= [TIMES, ROMAN, 20], 

417 thickness=3)) ; 
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Coraments 

331-334 This procedure will plot any field in the list Get .fields (or combination 
thereof) against tau for a given set of inputs. We can plot the field due to the 
straight trajectory by setting thetap= 0. 
345-361 This procedure will plot any field in the list as a function of TO. 



362-417 This will make the plots given in figure 7.3 



G.4 Part 5 - Convolution 



418 rho_ box:= (t,a, b) -> l/a*(Heaviside(t+a/2+b)-Heaviside(t-a/2+b)) ; 

419 plot(rho_ box(t, 0.0005, 0) ,t=-0.01 . .0.01,title="box 

420 distribution" ,colour=brown,axes=boxed) ; 

421 rho_ Gauss := (t, a, b) -> l/(a*sqrt(2*Pi))*exp((-(t-b) "2)/(2*a"2)) ; 

422 plot(rho_ Gauss(t, 0.5, 0) ,t=-l . . l,title="Gaussian 

423 distribution" ,colour=brown,axes=boxed,numpoints=10000) ; 
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Part 5 - Convolution 



424 conv:=proc(PEAK, a, b, N,comp ) 

425 local t, i, E_seq, tau_seq,rho_seq,EO_seq, conv,suml ; 

426 global convxl, convyl, convzl, convx2, convy2, convz2 ; 

427 if PEAK=1 then 

428 for i from to (N-1) do 

429 t: =16. 6667; 

430 #solve(a+((b-a)/N)*(i+l/2)=TT,tau) ; 

431 #print(" ",tau_l_||i 

432 :=solve(t-a-((b-a)/N)*(i+l/2)=10"(12)*TT,tau) ; 

433 #print(tau_l_| I i) ; 

434 EEE_0_||i :=evalf (subs(tau=tau_l_| I i, EEllcomp)); 

435 EEE_lJ|i :=EEE_0_| |i*evalf (rho_Gauss(t-a-((b-a)/N)*(i+l/2), b-a, 

436 t ) ) ; 

437 rho_l_| |i:=rho_Gauss(t-a-((b-a)/N)*(i+l/2) , b-a, t) ; 

438 end do : 

439 suml:=add(EEE_l_| |i, i=O..N-l); 

440 convl IcompI |PEAK:=suinl/add(evalf (rho_Gauss(t-a-( (b-a)/N)*(i+l/2) , 

441 b-a, t)) , i=0. .N-1) ; 

442 print (convl I comp I I PEAK) ; 

443 elif PEAK=2 then 

444 for i from to (N-1) do 

445 t : =16 . 685 ; 

446 tau_l_| |i :=fsolve(t-a-((b-a)/N)*(i+l/2)=10"(12)*TT| |PEAK,tau); 

447 #print(tau_l_| I i) ; 

448 EEE_0_ Mi: =evalf (subs (tau=tau_l_ I I i , EE I I comp I I PEAK) ) ; 

449 EEE_l_||i :=EEE_0_| |i*evalf (rho_Gauss(t-a-((b-a)/N)*(i+l/2) , b-a, 

450 t ) ) ; 

451 rho_l_| |i:=rho_Gauss(t-a-((b-a)/N)*(i+l/2) , b-a,t); 

452 end do : 

453 suml:=add(EEE_l_| |i, 1=0.. N-1); 

454 convl IcompI |PEAK:=suml/add(evalf (rho_Gauss(t-a-((b-a)/N)*(i+l/2) , 

455 b-a, t)) , i=0. .N-1) ; 

456 print (convl I comp I I PEAK) ; 

457 end if : 

458 end proc: 



Comments 

347-352 Defines the charge profile p{iy). We can use either a box profile or a 
Gaussian profile. 



424-458 Procedure for calculating the convolution (6.52). The convolution has to 
be evaluated for the pre-bent path and for the straight path for a selection of 
different bunch lengths. We adopt a Gaussian form for pLab and define the bunch 
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Part 5 - Convolution 



length as the full width at half maximum (FWHM). The results are given in table 
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